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FOREWORD 

This study was undertaken by Northrop-Huntsville , for the NASA/MSFC 
Huntsville, Alabama, under Contract No. NAS8-25182. The program was under 
the direction of MSFC's Flight Data Statistics Office with Dr. F. R. Krause 
as the project monitor. 

The authors wish to express their acknowledgement to Dr. Krause for his 
many helpful discussions m the completion of this work, and to Dr. R. R. 
Legault, University of Michigan, for kindly providing us the multispectral 
data presented in this report. Thanks are also due to Messers W. Campbell and 
R. Brecht, NASA co-op students, for helping with data processing. 
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Section I 

INTRODUCTION 


As more sophisticated space vehicles and orbiting space stations are de- 
veloped, onboard data processing and conditioning becomes an imperative system 
design consideration. The outputs from hundreds or thousands of various sensors 
in advanced mission payloads are fast exceeding the capacity of available and 
near-real-time communication systems. In particular, the airborne sensors for 
earth resources and advanced surveillance have larger amounts of data output 
as requirements increase for greater area coverage and resolution. 

Remote sensing of the earth's resources and environment from airborne and 
satellite-borne sensors is becoming one of the most significant and immedi- 
ate applications of the developed space technology to the whole world. The 
entire spectrum of remote sensing technology will cover sensor design, ground 
truth study, atmospheric physics, data handling, and data analysis. This study 
is concerned with data analysis only. 

The purpose of data analysis in remote sensing problems is to assess or 
identify some desired characteristics of the sensor's targets through analysis 
of multispectral (or multimage) data collected by the sensors. These character- 
istics may be surface temperature, humidity of the ground, salinity of the 
ocean, or classifications among different crops and soils. Different applica- 
tions mainly require different sensors designed to operate at different ranges 
of the electromagnetic spectrum, while the methods of data analysis are quite 
similar in principle. Therefore, it is not necessary to specify any particular 
applications in order to discuss the data analysis techniques . 

Since these techniques are to be used for identifying desired character- 
istics from some given data measurement, they are really within the region of 
pattern recognition problems. These problems have a long history of develop- 
ments especially in the field of experimental psychology, neurophysiology of 
the visual system, and automatic machine character recognition. 
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Handling and processing of the vast amount of data, which are usually ran- 
dom processes, will require compression of the data to smaller quantities of 
more meaningful parameters. These parameters should characterize the physical 
phenomena recorded by the various sensors and usually include the mean value, 
first few orders of central moment, correlation functions, spectral density 
functions, probability density function, and associated error functions with 
respect to each evaluated parameter. To handle many sensor outputs simulta- 
aneously would require multi-channel data processing. Further, to perform 
the data processing in a continuous and near-real-time basis, on-line sequential 
techniques are necessary. 

Statistical methods for the analysis of aircraft survey data are being 
developed by the University of Michigan, Purdue University, and the University 
of Kansas. All these methods require human intervention to select homogeneous 
training areas on the ground. Spectroscopic reference signatures are then 
computed for each training set and all other resolution elements are classified 
by matching them with all the reference signatures. The University of Michigan 
uses the correlation computation for the match, i.e., by computing the corre- 
lation coefficients between the sample spectrum and the various reference spectra 
(ref. 1). Purdue University uses more complex Bayesian statistical decision 
methods for classification (ref. 2). The University of Kansas uses even more 
complex clustering techniques which employ both measurement space and spatial 
clustering (ref. 3). 

Our classification differs from the above techniques in using methods of 
change analysis for the classification of resolution elements. Human inter- 
ventions for the selection of training areas and reference spectra axe not 
needed. The strip map of an electromechanical scanner is subdivided into 
inventory areas whose resolution elements are homogeneous. 

Statistical methods for data flow compression and signature enhancement 
are being developed to improve the utilization of remote earth observations 
from aircrafts and spacecrafts. The objectives of our present research 
program are to extend the automatic data flow compression and signature 
enhancement that has previously been developed for earth based sensors to 
flight systems. 
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Section II presents analytical consideration and -mathematical algorithms 
for automatic analysis of multispectral data. Section III describes in detail 
the computer programs developed based on mathematical algorithms presented 
above. The application of these computer programs to some actual remote 
sensing data by the University of Michigan multispectral scanner from a air- 
craft over some agricultural fields is presented in Section IV. Finally, 
the summary and conclusions are given in Section V. 
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Section i! 

ANALYTICAL CONSIDERATIONS AND ALGORITHMS 

2.1 INTRODUCTION 

Analytical considerations and the corresponding classification algorithms 
for multispectral data sequences are presented in this section. In order to 
have a better understanding of the overall problem involved, a survey of pattern 
recognition problems will be given first. This will be followed by (a) statis- 
tical analysis of feature characteristics, (b) statistical test for inventory 
boundary, (c) boundary enhancement by area smoothing and differentiation, and 
(d) similarity criteria betxjeen different populations. An overall functional 
block diagram of the data analvsis is shown in Figure 2-1. _ . 


2.2 PATTERN RECOGNITION PROBLEMS - A BRIEF SURVEY 

Pattern recognition techniques are methods for making decisions based on 
statistical testing. The purpose is to identify a specified pattern from a 
set of given input measurement or to classify a given input to one of several 
prescribed categories . The input measurement is usually composed of a set of 
N numbers of ordered entries (or components) (i = 1,2, . . ,N) , which may be 
considered as a feature vector X in N - dimensional feature space. The pattern 
recognition technique is to be used to classify X to one of M categories Cj , 
j = 1, 2, . . .M, based on some criterion of similarity between the input 
feature and the given patterns. 

The first important step in any pattern recognition problem is to know the 
characteristics of the feature vector from a sufficiently large ensemble of 
presentative samples of every possible category. C^. The necessity of this step 
is obvious, since, without knowing the characteristic or invariants associated 
with every pattern category, there is simply no rational basis for identifying 
them. Some variations in the characteristics of the feature vector always 
exist and some uncontrollable factors further introduce errors in the measure- 
ment of the feature vector. Hence, the logical way to analyze the feature 
characteristics is by statistical methods. This will require evaluating the 
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mean, rms, probability density of each component, covariance matrix among any 
two pair of components, etc,, from the set of representative samples (usually, 
called a training set). From these parameters and quantities evaluated, some 
statistical models may be formulated. These models, expressed in functional 
form, may replace those often lengthy and volumous numerical values of the 
feature characteristics. Some particular statistical models such as Gaussian 
distribution permit simplification in the latter steps involved in the pattern 
recognition. 

The statistical theory states that if every component of a given feature 
vector is equally important, then the probability of success for recognizing 
a certain pattern will increase with the number of the components employed. If 
every component is not equally important, then the probability of success will 
depend not only on the number of components employed, but also on their relative 
importance. Further, a complete statistical description of an ensemble of 
samples require very large numbers of components to be measured from each 
sample. It is both impractical and undesirable to do so. This consideration 
brings out the second important step in pattern recognition, namely, the 
problem of feature selection. A rationale for feature selection may be stated 
as follows: "Given the desired accuracy, we want to find the optimum group 

(or groups) of a predetermined number of components from the feature space so 
that the probability of error or mis re cognition be a minimum." Actually, the 
optimum feature selection depends on the error criterion chosen, feature 
characteristics, and algorithms for classification. 

Assume one has been given the optimum set of measurement vectors, the 
next step in pattern recognition is to devise some rule by which to classify 
any given feature measurement into one of several pattern categories. So this 
is the classification problem. One of the simplest classification methods is 
to compare the relative distance between the point (in feature space) to be 
classified to every representative mean point of available pattern categories. 
The category with the shortest distance to the point under consideration will 
be assigned as the recognized pattern. Such a simple classification method is, 
however, quite limited. For example, there may be a situation where making 
a wrong classification on some pattern categories is more costly than others. 


2-3 



NORTHROP 


TR-860 


HUNTSVILLE 

Under such a situation, we want to weigh differently on different pattern 
categories according to their relative importance. Many more sophisticated 
classification methods have been developed based on testing hypotheses in 
statistical inference theory. One of the most powerful methods is the maximum 
likelihood ratio test. 

In summary, a pattern recognition system consists of three interrelated 
aspects: feature characteristics analysis, feature selection, and classification 

algorithms. Each of these aspects will be described in more detail in the 
following sections including corresponding pertinent bibliographies . A 
hierarchy of the pattern recognition program is shown m Figure 2-2. 

2.2.1 Feature Characteristics Analysis 

Because the feature measurement taken in natural situations always involve 
some influence by the environment, and because some variations also exist in 
most pattern categories, the logical way to assess the characteristics of the 
feature vector will be the statistical analysis. More specifically, since the 
feature vector is of multiple components, a multivariate analysis will be needed. 
There are many different statistical methods to describe the feature character- 
istics in the feature space such as joint probability density, correlation 
function, spectral density, and characteristic function (which is the Fourier 
transform of the probability density) , just to name a few. In general practice 
and for simplicity, however, it usually takes some combination of these different 
types of statistical descriptions in their respective lowest order. In other 
words, one will usually evaluate at least the following quantities from the 
representative set of available samples: 

• Mean vector, which is the ordered mean values of the feature components 

denoted by p = (u , p , ...p ). 

12 N 

• Correlation matrix, which is a matrix of mean products of any two pair 
of feature components , denoted by R „ . 

• Covariance matrix, which is a matrix of mean products of any two pair 
of feature components subtracted by their respective mean values, 

denoted by C. . . 

il 

• Univariate histogram, which is the probability density of a single 
feature component. 
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From the above evaluated parameters and functions some analytical models 
or approximation of the feature space of the pattern categories may be induced. 

When the multivariate Gaussian distribution is made, the ensuing analysis 
on feature selection and classification algorithms can be greatly simplified. 
This is also the reason why so many approaches developed for pattern recogni- 
tion system employ this particular distribution (or model) . 

2.2.2 Feature Selection 

Not every component of the feature vector (or the set of measurement) are 
equally important for classification. Obviously, those components peculiar 
to pattern categories are of primary importance. It is also desirable to use 
as few components as possible to achieve the prefixed accuracy of successful 
recognition. These two facts suggest the need for an objective study of select- 
ing the feature components. The feature selection is closely related to the 
performance of pattern classification. 

Generally speaking, the more feature components used, the higher the 
probability of successful recognition. So, the problem of feature selection 
might be more rigorously stated as follows: out of a given ph ysica lly access- 

ible or available set of feature components, to select an optimum subset of a 
prescribed number of feature components such that the probability of successful 
recognition is maximized. By this definition, the feature selection problem 
must be analyzed with respect to a particular descision process for pattern 
recognition. However, this logical approach for feature selection is very 
complicated. Therefore, several existing methods for the feature selection 
adopt different approaches based on various measures of effectiveness in the 
choice of components. 

One approach Is to employ a measure of distance between any two (arbitrary) 
probability density of pattern categories (ref. 4). This measure is called 
"divergence", which is defined as the difference between two expected values 
of the likelihood ratios . The second approach uses the idea of entropy or 
average information (ref. 5) as criterion for feature selection. This approach 
calculates a statistic obtained over a large sample set to yield a measure of 
goodness for the feature component, which can be considered as the mutual 
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information of the feature and the pattern classes. The third approach employs 
Karhunen - Loeve expansion which does not require the a priori knowledge of the 
probability density of the feature vector. The essence of the expansion is to 
project the original feature vector onto a less-dimensional vector space spanned 
by a set of deterministic orthonormal coordinate functions. The expansion has 
the following optimum properties: (1) it minimizes the mean-square error for 

the truncated series of -the expansion and (2) it minimizes the entropy function 
defined over the variance of the random coefficient in the expansion (ref. 5). 

2.2.3 Cl ossification Algorithms 

The classification scheme in pattern recognition is to divide the feature 
space into some prescribed number of subspaces, each of which will be designated 
for one of the pattern categories. Different classification schemes differ 
mainly in the criteria used to establish these subspaces. The separating 
boundaries between these subspaces are called decision surfaces which can be 
defined by a set of functions. These functions are called "discriminant func- 
tions" and are so chosen that one and only one of them will give the larger 
value compared to the others for all samples in a particular pattern class 
(i.e. , a subspace). That is, for a given feature vector to be classified, 
a finite number of discriminant functions are to be evaluated with this feature 
vector. The function which yields the largest value will assign the input 
feature vector to the pattern category for which this function is designated. 

This scheme may be expressed as follows : Any given feature vector X is classi- 

fied as a pattern class Ch (i=l, 2...,M), if 

g ± (X) > g_. (X) for all i ^ j, 

where g. (i=l, are the discriminant functions for class C . 

1 1 

The classification methods may be divided into two main groups; namely, 
the parametric method and the nonparametric method. The parametric method 
is more appropriate for classification when each pattern category can be 
characterized by a set of parameters obtained from the feature characteristics 
analysis of the training set. The nonparametric method is more appropriate 
when no simple realistic assumption can be made on the characteristics of the 
pattern classes. These two methods are discussed in more detail in the. following 
paragraphs . 
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Tlie parametric method may be further divided into two distinct groups; 
deterministic and statistical methods. The discrimination functions in the 
deterministic methods can be either linear or nonlinear in the input feature 
vector. The linear discriminant functions separate the feature space by hyper- 
planes, while the nonlinear discriminant functions divide the feature space 
into curvilinear hypersurfaces. The simplest linear deterministic method 
is the so-called correlation method (also called template matching) . This 
method evaluates the scalar product of the input feature vector with a set 
of reference feature vectors, each of which represents a particular pattern 
class. The feature vector will then be assigned to the particular pattern 
class that yields the largest scalar product. Another linear deterministic 
method is the minimum-distance method which uses the shortest distance between 
the input feature vector to the set of reference pattern vectors as the 
criterion. 


The three common statistical classification methods are Neyman-Pearson 
criterion, Bayes criterion, and minimax test. The Neyman-Pearson criterion 
is the simplest of all and demands the least a priori information. The Bayes 
criterion is to minimize the average cost of misrecognitxon, when the priori 
probability of occurences of pattern classes are known, and when_the average 
cost is a linear function of the absolute error probabilities. The minimax 
test is a Bayes test except that the former minimizes the maximum of the 
conditional cost of misrecognition. All three of these tests are optimum under 
the respective criteria and also use the liklihood ratio of the conditional 
probability densities of the pattern classes . These methods hold for any 
arbitrary distribution functions; however, great simplifications may be made 
if the distribution is Gaussian. 


Finally, we shall describe the nonparametnc methods for pattern classi- 
fication. By these methods, no parameter from the training set is evaluated 
a priori, rather the discriminant function (in the linear case) will assume 
the form of 


g. (X) = W 


N 

i-1 


V 


X. 

i 
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.th 


(j=l, 2 , . . . ,M and i«l, 2, ... ,N) , 

where W. (W , W , . . . ,W ) are the weighting functions for the j pattern 
3 3 -t ± Z J IN 

class. These weighting functions are to be determined by a systematic error- 
correction scheme from the training set. The error-correction scheme may be 
written as 


W = W. + cX 
1 3 

This states that each new weighting function is simply the sum of the old one 
plus a fixed fraction c of the input feature vector. It has been proven 
theoretically that within the finite number of error corrections, a suitable 
weighting function can be obtained, provided that the training set is separable 
by hyperp lanes . 


In summary, the advantage of the nonpar ametric method is that no assump- 
tion is needed on the feature space. However, they are applicable to those 
cases where pattern classes may be separable by hyperplanes , while the 
parametric method is not restricted by the separability requirement. 

2.2.4 Pertinent Bibliography For Pattern Recognition 

Some pertinent bibliography on the pattern recognition technique together 
with remote sensing of the earth's resources and environment will be given in 
Appendix A. This bibliography is classified into five categories. 

• Remote sensing of the earth's resources and environment 

• Pattern recognition in general 

• Feature characteristics analysis 

• Feature selection 

« Classification algorithms. 

This bibliography is not intended to be extensive, rather it is considered 
to be more informative for the reader. 


2.3 STATISTICAL TEST FOR INVENTORY BOUNDARIES 
2.3.1 Introduction 

Any statistical evaluation of multispectral observation requires a popu- 
lation of resolution elements that are statistically alike. Such a population 
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occupies a certain area of the strip map, which shall be called an inventory 
area. In many cases such inventory areas are continuous in space. Typical 
examples of such continuous inventory areas are clouds, snow fields, agri- 
cultural fields, forests, etc. The boundary of the inventory area shall be 
called inventory boundary which thus separate populations of homogeneous 
resolution elements. However, all statistical properties will change whenever 
an inventory boundary is crossed. Conversely, the change of statistical 
properties along a scan line indicates that an inventory boundary has been 
crossed . 

The above considerations provide a basis for the development of sequential 
and completely automatic methods for the selection of the inventory boundary. 
Our previous tests for environmental variations (refs. 6 and 7) recognize the 
change of statistical parameters with record length, furthermore, the position 
of an inventory boundary indicates a change of reflectance or thermal emission 
which is not tied to any particular application. Statistical methods which 
utilize inventory areas are thus applicable to crop surveys, yield surveys, 
and land use investigations. Automatic detection of inventory boundaries 
could thus support statistical data flow compression and signature enhancement 
methods that may be sued for many disciplinary applications. 

2.3.2 Formulation of the Algorithm 

The computer selection of inventory boundaries follows from a direct 
application of accumulative error curves (refs. 6 and 7). Consider a sequence 
of resolution areas, i = 1, 2,..., m, along a scan. Each of these resolution 
areas is assumed to be statistically independent. Any average ( ) ^ from a 

single resolution area then represents a statistical experiment. If all 
resolution areas belong to the same inventory area, then the accumulative 
average over these areas 


( ) 


m 


1 

m 


m 


l < ), 

i=l 


( 2 - 1 ) 


should become more accurate as m increases. The statistical error of such 
an accumulative average should thus become smaller with increasing m. 
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The test for inventory boundaries is now derived by monitoring the decrease 
of the actual error with accumulative number m, by comparing this actual error 
curve with the most probable error curve. 


The actual (or sampled) error of an accumulative mean, A ( ) m , can be 
computed from the deviation of the individual experiment ( )^ relative to 
the accumulative average over all experiments . 


A ( ) 


m 


— 1 y 

m(m-l) . ^ 
L—l 


( >4 - ( > m 

i m 


( 2 - 2 ) 


This actual error can also be used to construct a random variable; i,e, Chi- 
square variable with (m-1) degree of freedom. That is 


m(m-l) ( ) 


(m-1) = 


m 


(2-3) 


The value of a is constant as long as the resolution areas, i = 1, 2 ... m, 

are statistically independent and belong to the same population. Under these 

2 

considerations, the probability distribution of the variable x is also known, 

2 

The most probable x (m-1) value is given by 


X (m-1) 


= m-3 for m > 4 


(2-4) 


’most prob. 

This most probable value indicates the existance of a most probable error 

.2 


<5 2 = 
m 


m(m-l) 


(m-1) 


most prob. 
2 


m-3 

m(m-l) 


(2-5) 


The value of the unknoxna variance a can now be estimated by fitting the shape 

of the most probable error curve, 6 , to the actual error curve A ( ) . This 

tn m 

curve fit is chosen such that the relative mean square variation between these 
two curves become a minimum. Let M denote the last accumulative number. The 
above curve fit may thus be expressed by 


_3 

3o 


M 

l 

m=4 


A ( ) - 6 

m m 

6 

m 


( 2 - 6 ) 
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Substituting S m from the previous equation and differentiating yields 


2 

a 


M 

V A 2 

L , tn-3 ' 

l=4 


Y /m(m- 

m^ V m ~- 


1) 



(2-7) 


The actual error A ) should oscillate around the most probable error curve. 

The range of oscillations is given by the confidence factor of the Chi-square 
2 

X distribution. 


X l-p ( m ~ 1 ) 
/m-1 


a 

/in 


< A ( ) < 

— m — 


x p (m-1) 
/S=T 


a 


( 2 - 8 ) 


The values of 


X p (m-1)/ /m-1 and x^ (m-1)/ /m-1 are listed in Table 2-1. 


The above derivation is only for single channel data. It is natural to 
extend the test for inventory boundary to a combination of several channels. 
This extension is achieved by replacing the random variable (defined by 
equation (2-6)) for a single channel by the generalized random variable of 
combining L number of channels . 


X 2 (m-1) 

L 


m(m-l) 

L 


/ 


I 

1=1 



(2-9) 


Selected 

Channels 

Inventory boundaries for a combination of L channels are then derived by 
testing the following confidence interval. 


Xl-p (m - 1 > 
v4n-l 


< < 


L 

l 

1 


A 2 (— ) 


m, A 


a x 2 / m 


1 1/2 

r 1 


Selected channels 



( 2 - 10 ) 


The derivation of this extended test is exactly analogous to the derivation 
for the single channel. 
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Table 2-1. CONFIDENCE FACTORS FOR INVENTORY BOUNDARIES 





p ( m-l )// ii>l 




M 

p = 0. 005 

p = 0.01 

p = 0.10 

p = 0. 90 

p = 0. 99 

p = 0. 995 

3 

1.985 

1.821 

1.17 

.01 

0.032 

0.000 

4 

1.880 

1.752 

1.24 

.08 

0.082 

0.058 

5 

1.789 

1.681 

1.25 

.12 

0.170 

0.134 

6 

1.726 

1 .631 

1.25 

.46 

0.244 

0.203 

7 

1.668 

1.586 

1.24 

.52 

0.304 

0.262 

8 

1.626 

1.549 

1.23 

.56 

0.353 

0.311 

9 

1.593 

1 .521 

1.22 

.60 

0.394 

0.352 

10 

1.563 

1.494 

1.22 

.62 

0.428 

0.386 

11 

1.536 

1.473 

1.21 

.65 

0.457 

0.416 

13 

1.494 

1.435 

1.20 

.68 

0.504 

0.465 

15 

1.459 

1.407 

1.19 

.71 

0.542 

0.505 

17 

1.432 

1.383 

1.18 

.73 

0.572 

0.536 

19 

1.408 

1.362 

1.17 

.75 

0.597 

0.563 

21 

1.389 

1 .345 

1.17 

.77 

0.618 

0.585 

26 

1.351 

1 . 311 , 

1.15 

.79 

0.660 

0.629 

31 

1.320 

1.286 

1.13 

.81 

0.690 

0.661 

62 

1.253 

1.214 

1.10 

.87 

0.768 

0.755 

m ->■ “ 

1 

1 

1 

1 

1 

1 
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The above formulation is discussed for a one-dimensional data sequence. 
Actually, it is applicable equally well for 2-dimensional data sequences such 
as data from a strip map. The above algorithm may be employed to detect the 
inventory boundaries in two directions, simultaneously. 

2.4 BOUNDARY ENHANCEMENT ALGORITHM 

Another method for inventory boundary detection - boundary enhancement 
method - will now be presented. Basically, the idea of the boundary enhancement 
is to precondition the two-dimensional array of original data by converting the 
difference between adjacent samples into the corresponding absolute value or 
by taking the square of the difference. Since, at the boundary between any 
two homogeneous regions, one expects larger differences in the data difference, 
the method of the boundary enhancement will accentuate the existing boundary. 

Let the 2— dimensional multispectral data sequence be denoted by x(i,j,k), 
where i and j represent the scan line and the sample number along a scan line, 
respectively, while k stands ^for the channel number for each data sample. 

Further let 

i = 1, 2, . . . , M 
j = 1, 2, . . . , J 
k = 1, 2 , . . . , K 

The 2-dimensional moving average, which serves as a smoothing operator, is 
defined by 

N N 

y(i,j,k) = l I x(i+i', j+j ’ , k) 

(2N+1) i'=-N j'=-N 

where a(2N+l) x (2N+1) array of given data samples have been averaged. This 
operation is illustrated in Figure 2-3. 

Two types of the boundary enhancement have been used. The first type is 
based on taking the square of the difference of adjacent samples after the 
2-dimensional averages. Specifically, it is defined as 
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Figure 2-3. ILLUSTRATION OF 2-DIMENSIONAL MOVING AVERAGE OVER TARGET SCENES 
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Z 1 (i,j,k) = [y(i,j ,k) - y(i,j-d,k)] 2 
+ [y(i»j,k) - y(i,j+d,k) ] 2 
+ '[y(i»j»k) - y(i-d,j,k)] 2 

+ [y (i>3 jk) - y(i+d,j,k)] 2 

The second type is the same as the first type except replacing the squaring 
operation by taking the absolute value of difference. That is 

Z 2 (i,j,k) = |y(i,j,k) - y(i,j-d,k) j 

+ |y(i>j»k) - y(i,j+d,k)[ 

+ |y(i,j»k) - y(i-d,j,k)| 

+ |y(i»j,k) -y(i+d,j,k)| 

where d denotes the number of samples skipped in calculating the finite 
difference and it can be any positive integer. This operation is illustrated 
in Figure 2-4. 

The new sample sequence Z^(i,j,k) or Z 2 (i,j,k) is enhanced data in the 
sense that their values have been magnified if they happen to locate close to 
inventory boundaries as illustrated in Figure 2-5. 


In order to find a systematic and automatic selection of the inventory 

boundaries based on these enhanced data, the mean value Z 1 „(k), mean 

—2 — ^ 
square value Z (k) , and the standard deviation S (k) must be calculated 
1 ,/ 

for each channel over the entire 2-dimensional sample sequence as follows. 



1 

MJ 


M 

l 

i=l 


J 

I Z (i,j,k) 
3=1 1,2 


! 
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Figure 2-4. ILLUSTRATION OF ENHANCING DATA AROUND INVENTORY BOUNDARIES 
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The threshold value, T, for determining the inventory boundary is then defined 
by 

T(k) = Z lj2 (k) + a \,2 (k) 

where a is some arbitrary constant. If the enhanced data sample 
. Z(i,j ,k) T (k) 

for at least one channel, then that sample will be considered to locate on or 
close to the inventory boundary. On the other hand, if Z(i,j,k) < T(k) for all 
k, then the sample is not considered on the inventory boundary. 

In the above dxscussion, the 2-dimensional moving average and boundary 
enhancement have been treated as two independent operations. That is, any 
given data will be processed by the moving average first and then use the out-r 
put data to be processed by the boundary enhancement. Actually, these two 
operations may be combined into one. Further, in many applications one needs 
only to use d=l or 2. The algorithm to accomplish these required operations 
simultaneously is given below. 

Z 1 (i,j,k) = W^(i,j,k) + W 2 (i,j,k) 

+ W^(i,j,k) + W^(i,j,k) 

and 

Z 2 (i,j,k) = W 1 (i,j,k)j + W 2 (i,j,k) 

+ W 3 (i,j,k) + W 4 (i,j,k) 

Here, the expressions in the right hand side of the equations for d=l and d=2, 
respectively, are defined as follows. 
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For d=l : 


(i> 3 »k) 
W 41 (i,j ,k) 


7 

~j ’=— u x(i— ^ N— l,j+j * ,k) 
N 

I x(i-K,j+j>,k) 

j*=-N 


N 

W (i,j,k) = l [-x(i4-i\j-N-l, k) + x(i+i' ,j+N,k>] 
i'=-N 

N 

Wo(i>j>k) = £ [x(i+i' ,j-R,k) - x(i+i',j+N+l,k)] 

i’=-N 

N 

W (i,j,k) = l x(i+N,j+j ' ,k) 
j’=-N 

N 

w 4 o(i» 3 »k) = -J x(i+N+l,j+j ’,k) 

4 j '--N 


W 3 (x,j,k) = W 31 (i,j,k) +W 32 (i,j,k) 
W 4 (i,j,k) =W 41 (i,j,k) +W 42 (i,j,k) 


For d=2: 

N 

W (i,j,k) = - l x(i-N-2 ,j+j ' ,k) 
j'=-N 

N 

W (i,j,k) = - l x(x-N-l,j+j ' ,k) 
j’=-N 

N 

W /-1 »3 >k) = l x(i-N-l,j+j ,k) 

3 N 

N 

W 42 (i,3» k ) = l x(i-N,j+j ' ,k) 
j'=-N 
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w 3 (i,j,k) =w 31 (i,j,k) +w 32 (i,j,k) 

+ w 33 (i f j,k) + w 34 (i,j,k) 

W 4 (i,j,k) = W 41 (i,j,k) + W 42 (i,j,k) 

+ W 43 (i,j,k) + W 44 (i,j,k) 

Care has been taken in arranging the above equations such that the input 
data x(i,j,k) can be read in from magnatic tape in the most sequential manner 
with respect to data storage in tape for any channel k. 

2.5 STATISTICAL CHARACTERISTICS OF AN INVENTORY AREA 

Two methods for detecting the inventory boundaries in any given target area 
have been presented in subsections 2.3 and 2.4. Thus, the area was divided 
into some inventory areas, each of which represents a homogeneous population. 

How to describe statistically the properties or characteristics associated 
with each inventory area will be discussed in the following paragraphs. 
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In a two-dimensional inventory area, the multispectral data are expressed 
as x(i,j,k). Here, the scan line i may run from I to 1^, and for each 
particular i, the sample number j may run from J (i) to J (i) , and k represents 

S 6 

the channel number (or components) for each sample. 


The mean vector for an inventory area is then defined by 

l e J e (±) 

M(k) = ~ I l x(i,j,k) 


i-f „ j=J (i) 


where N is the total number of samples within the inventory area. 


The correlation matrix of the inventory area is defined by 

I J ^ (i) 

R(k ,k ) = — l l x(i,j,k ) x(i,j,k,) 

x i=I j=J (i) 1 z 

s s 


The covariance matrix of the inventory area is defined by 


I J (i) 
e e N 


c (kl>k 2 ) = ~r I l [xCi.jjk^ - H(k,)] [x(i,j,k 2 ) - M(k 2 )] 


i=3l s j=:J s Ci) 


= R(k l5 k 2 ) - M(k x ) M(k 2 ) 


The normalized correlation matrix is defined by 


vw = 


R(k l ,k 2 } 


yR(k 1 ,k 1 ) R(k 2 ,k 2 ) 


The normalized covariance matrix is defined by 


C n <k l> k 2 ) 


G(k 15 k 2 ) 


/cc^.k^ C(k 2 ,k 2 ) 


Finally, the univariate (first-order) probability density function is 
defined by 
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pt^ < x(k) < x 2 ] 


1 

N 


{ Number of samples whose value In k 
channel lies between x^ and x 2 


More conveniently, if x^ and x^ will be given in exact multiple of prescribed 
increment Ax, then the above probability density function may be written as 


Pi(x(k)) = p[(i- -|)Ax < x(k) < (i+ -|) Ax] 

t 

_ jL Number of samples whose value in k channel lies - 

N * between (i-1/2) Ax and (i+1/2) Ax. 

Clearly, the above quantities defined for the inventory area are just a 
few of the simplest and more significant statistical features that one may 
choose. 


In many applications, the multivariate Gaussian distribution is a quite 
realistic approximation for the multispectral data, ps is also predicted by 
the Central Limit theorem of the probability theory. One of the desirable 
features of this distribution is because the distribution is completely deter- 
mined by the mean vector and covariance matrix only. Since this distribution 
occupies such a central important position in the pattern classification 
problem, it is given below. 

-|C(k r k 2 )| K K 

0 h h 

1 k x =i k 2 =i 


p(x(k)) = 


f2 1 r’) K/2 lc('k..k„') 


exp 


^ I 


|C*(k 1 ,k 2 )|[x(k 1 ) - MCk^Hxd^) - M(k 2 )]j 

where | C(k^,k ,) | is the determinant of the matrix C(k. ,k 2 ) while |c*(k^,k 2 )[ 
is the cofactor of CCk^jlO . 


Another form of the expression is 

1 


p(x(k)) = 


(2u) K/2 |c(k 1 ,k 2 ) 


exp 


~ [x(k) — M(k) ] 1 C 


. [x(k) - M(k)] 


} 
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where the prime represents the transpose and C is the inverse of 

C(k^,k 2 ). Also 

C “V ,k ) = 1 — 2 - |C*(k ,k ) 

|c(k 1 ,k 2 )| 

2.6 SIMILARITY CRITERIA BETWEEN INVENTORY AREAS 

So far, we have formulated methods for determining inventory areas and pre- 
sented a statistical description for each inventory area. The next question that 
arises is how to determine the similarity or dissimilarity among those inventory 
areas that do not have any common inventory boundary. In other words, we want 
to know whether two inventory areas far apart belong to the same population or 
not based on their statistical characteristics described in the last subsection. 


There are many criteria one can choose to establish the similarity between 
two inventory areas that are characterized by the given statistical discriptions , 
such as mean vectors and covariance matrices. One of the most appropriate 
criteria available is the so-called "divergence" measure, which can be considered 
to be a dimensionless "distance" between two arbitrary multivariate probability 
density functions p^ (x(k)) and p^(x(k)), for i^ and inventory areas, 
respectively. Since the mean vector and covariance matrix can be calculated 
from the probability density function, this "divergence" measure does take 
into consideration all the information from the inventory areas of interest. 

This measure is defined as follows (ref. 4) 




P i (x(k)) 

E [ l0g ( P j(x(k)) } 


i 


] 


-E 


p.(x(k)) 
[ l0g ( Pj (x(k)) 



Thus, it is the difference between two expected values of the logarithms of the 
likelihood ratios of a ramdom input sample vector x(k) , (k=l,2, . . .K) . In other 
words, the divergence is a measure of the difficulty for discriminating between 
the two given probability density functions. The divergence possesses all the 
properties of "distance" or "norm" as well as several others: 
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• D(i,i) = 0 

• D(i, j ) 0 for 1 1 j* 

• D(i,j) = D(j,i) 

• If every component of the sample vector is statistically independent, 
then the divergence is equal to the sum of the divergences for each 
component, individually, i.e. 

K 

D(i, j ) = l D (i,j) 
k=l k 

• If one more component is added to the original sample vectors, then the 
divergence for the (K+l) - dimensional sample space will never decrease, 

i.e. 


D(i, j [ x(k)) <_ D(i,j | x(k), x(K+l)). 


We shall now consider serval special cases of the divergence measure. If the 

probability density function are Gaussian with the mean vector ik(k) and (k) , 

and covariance matrices C.(k, ,k n ) and C.(k, ,k„), then the divergence becomes 

x 1 l j 1 l 


D(i,j) 


JL 

2 




tr «j [C^(k^,k 2 ) - C (k^l^)] 
X[ ^j _1(k l’ k 2 ) " S i 


( 2 - 11 ) 


+ | tr |[C i " 1 (k 1 ,k 2 ) + C^ 1 (k l9 k 2 )] 

*[5^) - M_.<k 2 )] tM jL (k 2 ) - M j (k 2 )]'| 

where tr denotes the trace, the exponent -1 denotes the inverse, and the prime 
the transpose. 


Further, if the covariance matrices are the same, “ C_j - C (k^,k 2 ), the 
above expression reduces to 
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D(i,j) = tr 


C 1 (k l9 k 2 ) [M.Ck^ - M (k 1 )] 


x[M (k 9 ) - M. (k ) ] 

i 2 3 2 


■} 


( 2 - 12 ) 


[M i (k 1 ) - M (k 2 ) ] ’ C ~ 1 (k 1 ,k 2 ) [M.(k 2 ) 


M j (k 2 )3 


If the components are statistically independent, then equation (2-11) can 
be reduced to 

D(i,j) = ~ l j[aj(k) - a 2 (k)] 2 + [a?(k) + o?(k) ] 

k=l L 3 3 


x[M.(k) - Mj (k)]H / [a 2 (k) a 2 (k) ] 


(2-13) 


where a^(k) is the covariance vector which is the diagonal of the original 
covariance matrix for i^ inventory area. 


2 

Further, if the covariance vectors are equal, i.e., o^(k) 
then the above expression reduces to 

D(i,j) - f 

k=l „ 

o (k) 


0j (k) = o 2 (k) , 


(2-14) 


On the other hand, if the mean vectors are equal, then 


D(i,j) = 


1 

2 


K 

I 

k=l 


2 2 2 
[°j(k) ~ 

0 2 (k) 0 2 (k) 


(2-15) 


Finally, the ordinary Euclidean distance between two vectors is obtained 

2 

from expression (2-14) by letting a (k) => 1, or from expression (2-12) by 
letting C(k^,k 2 ) equal to the identity matrix. 

D(i,j) = f [M(k)-H(k)] 2 
k=l J 

From the above various special cases, it can be seen that the divergence 
measure is to very general criterion for similarity. 
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Section HI 

COMPUTER PROGRAMS DESCRIPTION 


3.1 INTRODUCTION 

Computer programs based on the algorithms discussed in the last section 
have been developed for processing and analysis of multispectral data. The 
purposes of these programs, their input and output formats, and the program 
listing will be presented in this section. A list of these programs is as 
follows : 

• Tape Conversion Program 

9 Alphanumeric Plot Program 

• Statistical Inventory Boundary Program 

• Boundary Enhancement Program 

• Mean, Correlation, and Covariance Program 

• Piecewise Correlation and Covariance Program 

e Probability Histogram Program 

• Divergence Matrix Program 

• Si- Table Program 


3.2 TAPE CONVERSION PROGRAM 

This program is devised for converting a standard FORTRAN and/or binary 
tape to the standard format of input data tapes to the following programs: 

• Boundary Enhancement Program " 

• Alphanumeric Plot Program 

• Statistical Averages Program 

• Piecewise Correlation and Covariance Program 

• Statistical Inventory Boundary Programs 
0 Mean, Correlation, and Covariance Program 

• Probability Histogram Program. 

The outputs of the program is a FORTRAN binary or BCD tape with a control 
parameter (fixed-point binary or BCD) and data sets of some specified number 
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of channels (floating point binary or BOD) . The data set is a 3 -dimensional 
array, DATA (i,j,k), with i designating the scan line number, j the same 
number along each scan, and k the channel number. The control number is set 
to be 1 at the end of each scan line, 2 at the end of each specified field, 
and 3 at the end of each data set. 

It is optional in this program to read a header identification record at 
the start of each data set and also optional to write this on the output tape 
in either BCD or FORTRAN binary of some specified length. 

Because of the characteristic of most experimental data sets to contain 
erroneous data samples from some undetermined origin, it becomes desirable to 
replace these erroneous data samples with some data samples that best characterize 
the data set. The capability to replace these erroneous data samples is optional 
in the program. 

The control parameters in this program are inputed as NAMELIST parameters. 

The documentation on the usage of this system package is in IBM Manual C28-6390-2. 
There are two major sets of NAMELIST calls in this program NAMDAT and NAMFIL. 

The following is the definition of the input parameters used in NAMDAT and 
NAMFIL. 

NAMDAT 

o NFILE - Number of files of data to process from the input tape. 

• NWHICH - Number of the selected files to be processed. 

• LTN - Logical tape number of the input tape to be reformated. 

« LTNl ~ Logical tape number of data reformated by the program. 

9 LTN2 - Logical tape number of output header identification record. 

9 NSCANS - Number of scan lines to be processed in this field. 

9 IDWORD - Number of words in the input tapes header identification 
record . 

9 LLL - Data sample to start with from input data record in generating 
the output tape. 

9 NBEAMS - Number of channels of data. 

9 IHEAD - Option to read or write a header identification record. IHEAD 
not equal zero executed. 
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• IBCDH - Option to read header identification record in BCD. IBCDH ^ 0 

header record read as BCD information. IBCDH = 0 read as Binary. 

• IHMRT - Option to write header identification record on output tape. 

IHMRT / 0, record written. 

9 ISKBIN - Option to determine if any files are to be skipped on input 

tape to get to data file. ISKBIN = 0, no files to be skipped. 

9 IBCDW - Option to write output data tape in BCD. IBCDW ^ 0, tape written 
in BCD. IBCDW = 0, tape written in Binary. 

• IDROP - Option to replace selected data samples in data set. IDR0P ^ 

0, data samples replaced. 

NAMFIL 

9 XMEAN - Variables used to replace erroneous samples in first region of 
field. 

9 XMEANl - Variables used to replace erroneous samples in second region 
of field. 

9 XMEAN2 - Variables used to replace dropouts in third region of field. 

9 ISAM - Specifies location of regions in the data set to be replaced 

by means. 

9 XMIN - Minimum value of each region for a field. 

9 XMAX - Maximum value of each region for a field- 


Figure 3-1 is a flowchart of the program and the following is a listing 
of the program. 


3.3 ALPHANUMERIC PLOT PROGRAM 

The program is devised for printing a two-dimensional alphanumeric plot 
for grey-level display or classification results. The input tape can be BCD 
or binary. The other input parameters are the upper and lower class limits 
and desired alphanumeric symbols for each class . 

Definition of Input Data: 

CARD 1 (Al) 

IALPHA - Alpha-numeric characters for printout 
LLP - Border designation. 
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C 


TAPE 

DIMENS I ON 
DIMENS I ON 
DIMENSION 
DIMENS I ON 


CONVERSION PROGRAM 

ISAM ( 6 ) 

XMEAN ( 12 ) , XMEANI ( 12) » XMEAN2 ( 12 ) 
DATA ( 12 ) 

DAT (5000) ,NWHICH(4) 


DIMENSION RCDID(20) 

DIMENSION DUM (12) 

NAMELIST /NAMDAT/ NFI LE » NWH ICH, LTN ,LTNt , LTN2 , N$C ANS , I DWeRD, 
lLLL T KKK t NCT,NBEAMS, I HE AD» I BCD , I FLD , I SUB ,NA , N A8 . N B » I BCDH , 
2IMRD» IMRT, IHMRT,ISK8IN, IBCOR,IBCPW, I DROP 
NAMELIST /NAMFIL /XMEAN , XMEANI, XMEAN2_» I SAM, X MI N » XMAX 

I ZER0=O 


I ONE- 1 
I TW0=2 


ITHREE=3 


NDUM-0 
DO A J “ 1 » I 2 

4 DUM=0 . 0 

READ (5, NAMDAT) 

WRITE! 6, NAMDAT ) 

REWIND LTN 
REWIND LTNi 
D0 998 IAA*1,NFILE 
READ(5, NAMFIL ) 

WRIT£(6, NAMFIL) 

C SEARCH FOR FILE T0 PROCESS 

IF ( ISKBIN.EQ.O) GO TO II 
REWIND LTN 
I AB=NWH ICH ( IAA)“I 
IF{ IAB) 11, ll, 12 
12 CONTINUE 

DO LO I AC= I , I AB 
CALL SKFB I N ( 9 » 1 , RD ) 

10 CONTINUE 

11 CONTINUE 

READ RECORD ONE 

READ HEADER IDENTIFICATION RECORD 

IF(IHEAD) 14,9,14 
14 CONTINUE 

I F { IBCDH.NE .0 ) GO TO 5 

READ (LTN) (RCDID(I)»I=1» I DWORD ) 

I F ( IFLD.EQ.O) GO TO 6 
NCT=0 

CALL FLC (NCT,NA,NAB,NB,RCDID( ISUB) ) 

GO TO 6 

5 CONTINUE 

READ (LTN, 1006) (RCDIDl I ) , 1=1, I DWORD) 

6 CONTINUE 

I F ( IFLD.S'E.O) GO TO 21 
RCO I D ( ISUB )=0.0 
RCDIDI ISUB )=NCT 

21 CONTINUE 

IF (IHMRT) 22,23,22 

22 CONTINUE 

WRITE (LTN2) ( RCD ID ( MM) , MM= l , I DWORD ) 
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23 CONTINUE 

WRITE (6, 1003) {RCDID(K)»K=1, 10 WORD ) 

9 C0NTINUE 

C READ INPUT DATA 

NWW=NCT*N8EAMS 
DO l I=1,NSCANS 
IF{ IRCDR -NE.O) G0 T0 8 
CALL REOTPCtLTN, IMRD, I ERR » NW,NWW , DAT ) 

GO T 0 {50, 51, 50,50,51 ), IERR 
51 WR I TE ( 6 » 1004 ) I 
GO T 0 999 
50 CONTINUE 
0 CONTINUE 

LL=LLL 
KK=KKK 

13 CONTINUE 

C REFORMAT INPUT DATA 

DO 2 J=1»NCT 

I F ( I8CDR • EQ. 0 ) GO TO 15 

READ { LtNi , 1005 ) I U ( DATA { 1 1 1) » 1 1 1 =1 f NBEAMS ) 

NNN = 1 

DO 31 NN=LL,KK 
OAT( NN ) =OATA( NNN ) 

NNN=NNN+l 
31 CONTINUE 

15 CONTINUE 

C REPLACE ERRONEOUS DATA SAMPLES 

IFIIDROP.EQ.O ) GO TO 160 

IF (DAT (LL ).LT. XMAX ,AND«DAT(LL).GT.XMIN) GO TO 160 
IF {J.GE. ISAMU). AND.J.LE.ISAM(2) ) GO TO 15C 
IF (J.GE.ISAM(3 )*AND*J.LE. I S AM( 4 ) ) GO TO 152 
IF {J.GF.tSAM(5).AND.J.LE.ISAM(6)) GO T 2 154 
GO TO 160 

150 I Z=LL 

DO 151 K = 1 » KKK 
DAT ( I Z ) =XMEAN {K ) 

151 IZMZ+1 
GO TO 160 

152 I Z=LL 

DO 153 K=l , KKK 
DAT ( IZ)=XMEANl(K) 

153 IZ=IZ+1 
GO TO 160 

154 I Z-LL 

DO 155 K=1 i KKK 
DAT ( I Z ) =XMEAN2 { K ) 

155 IZ=IZ+1 

160 CONTINUE 

I F ( IBCDW • EQ. 0 ) GO TO 61 

WRITE(LTN2, 1000) I ZERO, l DATA{K ) ,K = LL, KK ) 

GO TO 64 
61 CONTINUE 

WRITE (LTN2) I ZERO, ( DATA { K ) *K = LL t KK ) 
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64 CONTINUE ' ' 

LL-KK+1 

KK=LL+N BEAMS- 1 
2 C0NTINUE 

I F < I.EQ.NSCANS) G0 T0 3 
I Ft IBCDW.EQ.O) G0 TO 65 

W R I T E ( L TN2 » 1 0 CTO ) 'I 0NE , ( fJOMMYTR ) ,R = 1 7NBE AMS ) 

G0 T0 1 

65 CONTINUE 

WRITE!LTN2 ) I0NE , ( DUMMY! K ) ,K=1 , NBEAMS ) 

1 CONTINUE 

3 CONTINUE 

I F { IBCDW.EQ.O) GO T0 63 

WRITE(LTN2 f 1000) IT WO ♦ ( DUMMY ( K ) f K=l » NBEAMS ) 
GO T 0 998 
63 CONTINUE 

WRITE (LTN2) IT W0 ,( DUMMY ( K ) ,K=1» NBEAMS) 

998 CONTINUE 
IFtlBCDW . EQ. 0 ) G0 T0 67 

WRltE(LTN2, 1000) I tHREE » ( DUMMY! K ) »K=I »NBEAM 
G0 T0 999 
67 CONTINUE 

WRITE(LTN2 ) I THREE , (DUMMY! K ) ,K=1 , NBEAMS ) 

999 CONTINUE 

END FILE LTN2 
REWIND LTN 
REWIND LTN 1 
REWIND LTN2 
200 FORMAT ( 2014) 

1000 F0RMAT( I it 8X» 12F6.3 ) 

1001 FORMAT ( ID 

1002 F0RMAT{ IH1 //2X, 24HINPUT RECORDS BAD - ST0P) 

1003 FORMAT (lX»12A6»I6 t 5X»2F8.4) 

1004 FORMAT ( 2X» 14) 

1005 FORMAT! I1»8X, 12F5. 3) 

1006 FORMAT ( 20A6 ) 

1007 FORMAT! Ii t 8X^12F8.3) 

1008 FORMAT! lX~, 2110) 

STOP 

END 

SNAMDAT 

NPILE=4 

NWHICH=2*8fl9 f 12 
L T N=9 
LTN L = 9 
ISKBIN=1 
I BCDR=0 
I BCDW= 1 
I DR0P= 1 
LTN2=1 1 
NS CANS =2 7 
I DWORD-15 
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LlL = l ' 

KKK = 12 
NCT = 167 
N8EAMS=12 
IHEAD= 1 
. IBCD=0 
IFLD=0 
I SUB = 1 3 

NA= 1 9 ' ' 

NAB= 17 
NB = 1 
IBCDH=0 
1 MRD = 1 
I MRT= 1 
I HMRT -0 
$END 
iNAMFIL 

XME AN = 1. 10 35 8, 1.6092, 1.1 33 19, 1.02380, 1.04938, 0.99 5 Alt 0.77604, L. 0 35 85, 
0.50232,0.63149,0.81489,0.63113 

XME AN l -0.92628, 1 . 36367 , l .00306 ,0 . 95745 , 0 . 90845 , 0 . 7667 3 , 0. 636 C9 , 1 . 039 73 

0.59239,0.59394,0.39618,0.30986 _ ___ 

XMEAN2=1.27359, 1. 70424, 1.143 78,1. 0214 i;o. 94322, 0.82919*, 0.621 34, 0.83949 
0.40804, 0.46054,0.58812,0.43592 
ISAM=1,63, 63, 118,118, 167 
XMIN=0.0 

XMAX=3.0 

$END 

iNAMFIL ' * 

XMEAN= 1.50995, 2. 15190,1.52394,1.34540, 1.29378,1.15872,0.92575, 1.24317, 
0.63664,0.68725,0.76504,0.60184 

XMEAN1 =0.97784, 1 , 47237, 1 .08407 , 1 . 01881 , 0.976C9 , 0. 84300 , 0 . 708 77 , 1 . 13348 
0.63643,0.64560,0.44608,0.36001 

XME AN 2=1. 13874, 1 . 861 39, l .27976 , 1 . 13856 , 1. 03394 ,0 . 86392 ,0 . 68685 ,0 .968 13 
0.49192; 0.52221,0/49081 ,'0.38598 ' ~ 

I SAM = 1,42, 42, 111,111,167 

XMIN=0.0 

XMAX-3.0 

$END 

SNAMFIL 

X ME A N= 1.31144, 1.85 966,1. 28 023', Ti 13480 ,'T 71056 r,T.~CO 866\ 6.76919,1.02518, 
0.49790,0.56153,0. 72571,0.54945 

XME AN1 =0.98687, 1.46594,0.97921,0.85527,0.84487,0.84443,0.60006,0.75333 
0.33964,0.57555,1.18249,0.89077 

XMEAN2= 1.29660, 1 . 74870, 1 . 17140 , 1 . 0542 i , 0 .97544 , 0. 85978 , 0 . 642 8 i , 0 .89028 
0.42091, 0.49407,0.61681 ,0.46276 

I SAM* 1,39, 39, 10 8, 108,167 

XMI N=0.0 
XMAX=3 .0 
SEND 
iNAMFIL 

XME AN= 1.00572, 1.47489, 1.03701,0.96775,0.91567,0. 60819 ,0. 63198,0.99724, 
0.53699,0. 56149, 0/48'457, 0 /S6704 

XME AN 1=1. 16048, 1.58077,1.05821,0.94602,0.87191 ,0.76703,0.54934,0.78006 
0.37171,0.42442,0.53810,0.39657 
I SAM* l, 111, 111, 167,167,167 
XMI N =0*0 

XMAX=3 .0 ; 

iEND 3 -8 

I 
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CARD 2 (2F10.0) 

XSKSC - Number of scan lines to skip (i.e. every other, 2, 4, etc.) 

XSKFL - Number of samples to skip in scan line 

CARD 3 (714) 

ITERM - Set to 1 for new case. 

JCNT - Skip numbering across scan line. 

NFLSKC - Option to skip data files on input tape. 

NFLSFL - Option to skip data sample in scan line. 

NFLSKN - Option to skip scan lines. 

IAP — = 0 to initialize numbering of data samples across scan line. 

IPRNT - = 0 for data from boundary enhancement program. 

= 1 for raw data tapes. 

CARD 4 (414) 

N - Number of classes or intervals used, 

NPTS - Number of data samples in scan line. 

NCR - Number of channels used. 

MSTOP - Total data samples used in scan line. 

CARD 5 (F6.0) 

XLOW^ - Lower limit for each channel, i = 1, NCH 
CARD 6 (F6.0) 

XUPP - Upper limit for each channel. 

CARD 7 - Card NCH (Al) 

IWAVE - Wave length in microns for each channel. 

CARD NCH+1 (15) 

ICU - Channels used. 

CARD NCH+2 - CARD N (F10.0) 

XINCR - Increment size/channel for each N(card 4) 

The output format is shown in Figure 3-2 . 

The flowchart is given in Figure 3-3 which is followed by a program 
listing. 
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MULTI SPECTRAL SCANNER DATA 2X2 


CHANNEL I 

WAVE LENGTH' (MICR0NS) 0.40 - 0.4*' 


SYMB0L 


1 

2 

3 

4* 

5 

6 
7 
6 
9 
0 


INTERVAL 
0.50 - 0.65 

0.65 - 0.80 

6. 80 - 0.95 

0._95 - 1. 10 

l.IO’-' "1.25 
1.25 - 1.40 

1.40 - 1.55 

1.55 - 1.70 

1.70 - 1.85 

1.85 - 2.00 

2.15 - 


W0P"S C5- 1 7,5- 22*5-5*5-8* ALtlTUDE 2000 FT. DATE RECORDED 6/30/66 


11111111112222222222333333333344444444445555555555666666666677777777778888 

12345678901234567890123456789012345678901234567890123456789012345678901234567890123 

**** *****♦******♦***♦**»»»»****»*****»»***»♦*»**»*»»*»»»»»♦»»*»»»»»♦♦»♦« »»♦»»«♦»»♦♦»♦ 

1 *555555 55666555555655554555444444444443334433 3333334333443344455^5665665666766677767 

2 *55545555565555^55555544555544443333433434333333333333333334345655556556666666666666 

3 *55544555566555455555554555543434443433333333333333343333333345555555556666666666666 
* *55545555666555555455544455544343443333333333333234343333333345556555565566666666666 
5 *55555555555555555554554455534443334434333334343333333333333345555556655665666666766 

6 *55545565555555555555554555544444444333343434343434334333433345556556566666666567666 

7 *55545559655445455544554444543343444334333343333333344333333345555565655565666666666 

8 *555555556555455555 5554 4445534444433333333434343343434333333345556555666666666666667 

9 *555 ^55 5 56553565355554444555443444 3434343 3433433 33 344443 334435 556555556666 66 76666777 

10 *55555555555445 455554 554455544444444443333344433333444343343345566555555566666666666 

11 *55 55 5555555 545455555 5 454455 343344 334443 3 3334433 3333 3333 333334555656656666 6666666 666 
,,12 *55545555555454445445444445534344444443333333433333443333333345555556556565666666766 

13 *555445 56555555<f5 5554454445544444444444433344433333343333333335555545555566666666666 

14 *77888888888888877777784455445444444444443443333434344444577677777777777777766656555 

15 *778877778877788/7 776784444444434443333344333334443433'334566666767676776777766'655555 

16 *77787787878888877767784454454444444444434433334334 4344 44566776666676676667666655555 

17 *7778788 BST7T7TT887777844455444444443344434333434'44344343467766667666777676766555555 

18 *77778878778777877776784444434444434444344333333433334333466777666766666767666555545 

19 *677 7 78 888 88 8 788 77 777 794455444444444444444443 3334343 3344346667666 7677/67 77766655555T 

20 *77877877777777777777684444444344343444343334333343333343566666666666776766766555545 

21 *77777777777B77m'6776T44444443lf44344W3“4T3*3TJ442«3TJW66666667776676666666555544 

22 *77878778877787887777784445444444333444443333333333333444466666766666666666666554545 

23 *78888S78777777877T7 67744444444444444444 34333*343 3433 3443446666666 767777766 666665 5 545 

24 *78888767677778787676674444444443444444343343334333433333466666666666667666666555555 

25 *7 788B8 8 87 78 8 787 87 7777 84445444444844444344343 33333443444446666766776767766 6666545545 

26 *77778878778778777776774444444344333544433344331334334334455666666667666666656545544 

27 *76677677666666S6&6658655444'4444'44443 43434434T43444333 344555566656666656566676666666 

28 *66676667666566556554564444344433333434444433434434343444556665665556556566666666566 

29 *666 5 76 66666 655666~S6566'44S4444343344444444443 3133334 34*43455 66565 665 565665 5677665665 

30 *67666777766655557554664444444444444434343343443444454434555665566566666666777666665 

31 *667677 66566565666565654444454444444433444343334444344433556666665656656565666655566 

32 *66667675666566565665664445334444443444444443344434434343566656565666656656566665655 

33 *66666776566556666665664444444444343444334433333334334344467566566666566557666665666 

34 *66666776776666665765654444454445433444444333344443444444556666665666665656676667566 

35 *567? ?66T7 766565 556655 54445444444444 34 3344433 43434444343 3556565566665566666676666656 

36 *66655766666666566654664445444443334443444343333443333344465665556665666666677677655 

37 *666667646666665655656544444344443344344333434434444444455W6'7555i5656566566666?66665 

38 *67656776666666666655665445444444344343444443344443444334556665556666665666676665655 

39 *666666 77665665666665654454444434433434434334 1343344 344434566666656655666*6776676666 

40 *44445444333433333232222332333333322336545554566666666656666765664565567444344444344 

41 *444 5 5?T313"5H3H!3Ty32 322233212233233355455515665565 66651T566766655S 65556444444444433 

42 *44445543343333333323233333333333233346445455666666666666666666665556667454344444444 

4 3 *44444533333323222232322223232233233335444545466655666655667765555565557443444433333 
94 *44445434333333323233 32 3222332323333336444545555656665666677666555566566434444444434 

45 *44555543433333333233333232333333333335544555666656666556666666555566668444444444444 

46 *454 55 54333232 333 3333 3 33232333232 2 33334544455565656655666667766665566558444444444*43 

47 *44545543343333333323322222333233333335545554556655666655666666655565657444444444433 

48 *45455544333333333333333333333332333334445545555566666665666665655565668454444444443 

49 *44555544333333433323334223233333333334545545565565666656666766655556667444444444444 

50 *4554554 43 3333332233232332333333333333454554SS66S555 566655666765555656 6744 4444333444 

51 *45555644334333333333433333233333333334545555566655666666666666656566567554444444444 

52 *55555643333 33333332343433233333333333554555566555566665566666555555 566754 5444444443 

Figure 3-2. MULTI SPECTRAL SCANNER DATA 2 X 2 
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ALPHANUMERIC PLOT PROGRAM 

61 mensiWTalph a i lT^oTTxtOH iilT7xuPPli3 jVxxlowT : 13 j » data f if o j 

DIMENSI ON IARAY t!70 ] »XDUM !13 ] , I ALP C 170 1 , JA RAY [108# 60 ) 

DIMENSION IDATAti5)#XJNCRtl3#2QJ# ICH [133 # XDWN 13,3] # I WAVE (13# 131 

IN =5 _ __ 

|N1*9 

LOUT = 6 

'DO 299 L-l* 13 

READ! IN , 99) 1 1 ALPHA [ L# H # I g U2QJ # LL P 

299 CONTINUE 

300 CONTINUE 

READ! IN, 114) XSKSC,XSKFL 

READ II N, 1D0 I I TERM* JCNT# NFLSKC# NFLSF L# NF LSKN# I AP# IPRNT 
iF[NFLSKC,EQ,0] GO TO 1001 

NFLSkTc GREAiER” THAN 0 "SKIP FILES 

DoToFo I si, NFLSKC 

CALL SKF BINEINIjIjRDI 

1000 CONTINUE 

IJOl fiONTLNUE 

READ EIN, 100 3 N , NPTS, NCH# MSTOP 

BEADJJ.Nj.il, 61, I X L OW [13 1 1 .T UNCH J 

READ! JN#ll6J IXUPP! I 3, I=1#NCHJ 

DO 1002 L sliNCH ^ _ 

READ! IN, 99) (IWAVEIL/U# I*1#13J 

1Q02 CONT INUE ' „ __ 

READ II N, 115) rrCHUJ#Ii*l»NCH) 

NUMB=N 

DO 600 L*1,NCH 

_ READ! IN, 114) I X I NCR !L# I) # I»1#NUM8 1 . 

DO 600 J=1 , NUMB 

X I NCR I L, J)=XUPP!L3/XINCR[L # J ] 

XXLOWtUsXLOWIU 

600 CONTINUE . 

XN = N 

IFLSKN = 0. ... 

ISCAN«0 

L=0 . 

NCALL=1 

INDEXES 

NSTART=1 

NSTOP= MSTOP 

JPASS=1 

WRITE I IOUT #110 ) . _ 

WRITE! IOUT#102) 

WRITEJL0U T,1 Q6 ) 

WRITE! IOUT/111) 

MSKFL = X_SKtl 

MPTS=NPTS/MSKFL+2 

DO 1 1=1# MPTS . 

I ALP E I ) = LLP 

1 CONTINUE _ . 
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c 

KKsN+1 


c 

c 

CALCULATE 

AND PRINT INTERVALS FOR EACH CHANNEL 


DO 309 L=1,NCH 
IFMCHILJ .EQ.03 GO 

TO 309 


WR ITE 1 1 OUT #103} JCHEU 

DO 2 1=1 i_KK 

IFtl.EQ.lJ GO TO 9 

I F [ U EQ.23 GO TO 3 _ 

4 WR I T£ 1 1 OUT .1053 I ALPHA EL, 1 3 , XDWNCL J /XUPP E L J 


GO TO 31 


9 

WRITE C 10UT.U21 
GO TO 32 

UWAVEEL.Kl.Ksl.13 3, 1 ALPHA tL> I 1.XLOWEL1 

3 

31 

WRITE! IOUT, 105} 
CONTINUE 

I ALPHA CL, n,XLOWELl#XUPP!U 

32 

XDWN (L 3 = XLOW EL J+XINCR EL* 1*13 
CONTINUE 

2 _ 

XUPPCll=XLOWELJ+XINCREL,Il 

CONTINUE 

ZUPP=XL0WIL3+XINCRCL*NUM8J 
WRITE! I OUT .105 3 I ALPHA EL‘, N + 2 3 » ZUPP 


WR TTE 1 I OUT #106 3 
CONTINUE 



IF [ IPRNT , GT. 0 3 

GO TO 315 


WRITE? IOUT, 102] 

writeiiout.hu 



WRITE I IOUT ,1181 
DO 31_0_.Lsl.iNCH'.. 



IFEICHCLJ.EG.OJ 
WR I TE E I OUT', 119 3 

GO TO 310 

ICHtU, i IWAVEEL.K3.Ksl, 13} 

310 

CONTINUE 

LZaKK+1 


1111 ^^ ^ ** 4 - - ■ 


WR ITE 1 1 OUT ,1201 

E ! ALPHA ELL. 1 1 , 1°1,LZ1 

315 

CONTINUE 
N$TRT1=1 


NSTOPlsNPTS 

GALL, I LABEL [MSTQP, NCALL, NFTS, JCNT, I OUT 3 . 


WRITE! IOUT. 106] 
JiR.LT.EJ I QUTjlLOM 

LLALPllLtl^LiMPTlL __ . ____ 


C 

C INUM COUNTS SCAN LINES 


C 


INUMaQ 


5 
C . 

CONTINUE 


c 

c 

SKIP SCAN 

LINES 

401 

IFLSKNsIFLSKN+ 1 
I NU.M a 1 NUM+1 



IF [NFLSKN ,EQ. 0 ] GO 

TO 205 
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SSs I FLSKN 

SMgS S/X SKSC 

MS=SS/XSKSC 

SS = MS_ 

SNsSH-SS 

I F [ SM . EQ « 0 . 0 ) GO TO 205 

DO 400 I = NST ART # NSTOP 
R E A D ( IN 1 3 I C HEK« CDATACJ1 ,J*1, N C H J 
iniCHEK.EQ.3J GO TO 210 

400 CON TINU E 

READ [ INI J ICHEK. C DATA [ J j , J*l, NCH,J 
_ IF [ ICHEK .EQ.31 60 TO 210 
GO TO 401 

205 CONT INUE __ 

KKsNUMB+1 

M -0 


START PROCESSING ONE SCAN LINE 

DO 90 1 gNSTART .NSTOP 

READ 1 1 N1 J IGHEK# t DATA C j J , J«l, NCH J 

7 IF! I CHEK. EQ.31 GO TO 210 

IFCNFL5FL.eQ.0j GO TO 403 

SS* I 

SNsSS/XSKFL- 

. ms* ssy x skel . . 

SS = MS 

$m=sm-ss 

IFCSM.NE.O.OJ GO TO 90 

40 3 CONTINUE 

M=M + 1 

NNXs|L 


—CLASSIFY DATA S AMPLE S 


D-Q— S-CLO-Jsl t.KK 

NNX*NNX*1 


DS 3Q1.Jji1jlM.CE . 

IFCICHCL J ,EG,0 J GO TO 501 

IFIPATA ILL. G X^XXL OJLC_L.Lf XI N g RJLj MU M 8J.3. GO TO 5 04 
IFCJ.GT.l] GO TO 502 


GO TO 505 


r.Q. 5JUL„_. 


502 I.FJ D A T A W C.L_Lt.XJ N C.R Cl , J_- 1 1 1 GO TO 500 

505 I ARAY CM J = 1 ALPHA CL. NNX J 

501 IF LL.E QjJ^LH J GO TO 90 


500 CONTINUE 


„ GQ TO 90 

504 IARAYCMJ3IALPHAIL#NUMB+2J 
90 CONTINUE . 

MSKFL^XSKFL 

MPTS = NPTSy.HSKFL _ . 


FOR E ACH SCA N L INE 
MAXIMUM OF 110 PTS' /LINE 
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C 

96 WR ITE II OUT * 104 3 INUM# U A_RAY t K ]_# Ksi, MPTS 1 , 1 ALP 1 1 J 
READ [ !N1 ] IDUH, [DATA! J]/j=l,NCHJ 

I FI I DJJM iGE j_3 3 GO TO 211 

IFCNPTS.GE.MSTOPl GO TO 5 

2jp.4_ U»0 _ 

I NDEX S I NDEX+1 

NSTOP2 = NPTS±l _ _ 

DO 206 KKsNST0P2,MSTQP 

. j.«U + t . 

JARAYt INDEX#U = IARAY£KK] 

206 CONTI NUE 

GO TO 5 

210 CONTINUE _ 

211 NCAL1_ = 2 

IHMPsMPTS+2_ _ _____ 

WRITE! I Out* 108] £ I A UP tU » 1 ■ 1 # I MHP 1 

IF INPTS.GE.MSTOP3 GO TO 216 

212 CONTINUE 

WR I TE ( I OUT *110 3 

CALL I LABEL t M STOP > NC ALL* NPTS# JCNT # I OUT ) 

NST0P1= MSTQP«NPTS 

DO 215 L=l» INDEX 

W R LT.E I IOUT #-1-13.1 I JA EAIi L> LL 3 i_LL.= 1 >. N S TOP 1 3 

215 CONTINUE 

REW-I-M.D_.IN1 

IF £ ITERM , GT » 0 3 GO TO 300 

jttsiTEJ-iauifciio i __ 

99 FORMAT C 25 A1 3 

„ JULO FORM A T.U 6„I 4J 

101 FORMAT [ II, I4*4X»12F5,3J 

102 F0RMATI14X» 46HINVENTQRY BOUN DARY MAP BY. BOUNDARY ENHANCEMENT// 1 

103 "FORMAT t2X,8HCHANNEL 1 2, 36X# 22HSYMB0L INTERVAL /3 

1.0 4 FORMAT C6X« I5* 1X*1H»*112A13 _ _ . 

105 FORMAT C50X»A1>8X#F7 *4# 3H ■» ,F7,4J 

_JLQ6 JEMMAIIim 

108 FORMATtlOX»ll5A13 

. .10.9 FQRMATt l2F10.0,12AH 

110 FORMAT tlHl/ 3 

U1 F 0 R M A T I 2X j 83HCRQPS C5»1 7,5 -22 ,5.-5,5 - _ A_L T I T U D E 2 0 0 0 _F T , 

1 DATE RECORDED 6/30/66 //] 

■ ■112 FORMAT. [ 2X « 22HWAVE LENGTH tMICRQNSl ,13A1#13X,A 1 eMjJSjOj- 

13H - / F5 , 2 3 

113 FORMAT C1 1X/112A1) 

114 FORMAT C8F10 ,03 

115 FORMAT I13 I53 ... 

116 FORMAT C13F6.Q 3 

_U_S FORMAT [20X*7HCHANNEL*15X*11HWAVE LENGTH/ 3 _ _ 

119 F0RMAT(22X,I2»14X#13A1] 

120 FORMAT [2X,7HSYMB0L Si5X« Al»18H = LESS THAN MEAf\|, 5X, A1,16H s MEA.N-M 
1EAN+1SD,5X, A1,20H = M6 AN+£sDpMEA'n+2'sD// 

218X* A1 # 2QH = MEANt2SD-MEAN*3SD«5X>Al, . 

325H = GREATER THAN MEAN+3SD 3 

END - " ’ ' 
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SUBR OUTINE I LABEL [ N, I CALL* NST OP * NCNT» KQUT I 

C 

C „ SUBROUTINE NU MBERS ACRO SS SCAN LINE 

C FOR LABEL'IND AND IDENTIFICATION 

C INUM HUNDRE DS 

C JNUM = TENS 

C KNUM s UNITS „ ... 

C 

DIMENSION INUM I ISO Ii JNUMf 15 0 3 « KNUM 115 0 ] 

IF C ICAIL «GT t lJ GO TO 500 

M=NCNT ... 

KsNCNT 


.L^Q 

DO 1 I sNCNT > N> NCNT 


. 

INUM (LI S I/100 
. __ IF IM.EQ.lOQl M s 0 _ 

IFfM.GT *1001 MsM-100 

JRU.Hg.Lr..t1/lP 

M=M+NCNT 

. _ IFIK. 66.10I KPO 

IFtK.GT.10J KsK-10 

KNU_MJ1J«J^ 

KsK+NCNT 


_1 CONTINUE 


ISTOP«NSTOP/NCNT 

WRITE CKOUT,li0J [ JNUM 1 1 1 * in, IsfoPj 


110 FORMATtilX, 18011 1 
.60 , T.Q.-2JLH 


500 K=NSTOP+l 
. „ K sK/NCNT 


M=N/NCNT 


.. mu EjK O U T / U .Q. 1 I I NIjjlLLLLlMjLWi 

WRIT6IK0UT,110 J t JNUM III # I »K, M 3 
m I TE [ KML tll li I K MilLU.t.I P .Ki.iU 

200 RETURN 


END _ 
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3.4 STATISTICAL INVENTORY BOUNDARY PROGRAM 

The inventory boundary program has been developed to automatically 
construct a map of inventory boundaries for sequential multispectral scanner 
data. Inventory boundaries are calculated simultaneously in both the scan line 
(horizontal) and sample number (vertical) directions. The mathematical formula- 
tions are given in subsection 2.3. 


DEFINITION OF INPUT DATA: 


CARD 1 

DOT - Symbol used where no boundaries exist. Usually a period (.). 

DASH - Symbol used for boundary in the flight line (vertical) direction. 

EYE - Symbol used for boundary in the scan line (horizontal) direction. 

PLUS - Symbol used where a horizontal and a vertical boundary coincide. 

CARD 2 

2 

LI - Lower percentage value (a ) for first x table. 

k 2 

U1 - Upper percentage value (a^) for first x table. 

L2 - Lower percentage value (a ) for second x" table. 

L 2 

U2 - Upper percentage value (a^) for second x table. 

IAVG - Number of samples to skip before comparing ERCK to confidence 
limits. 

MACH - Indicates machine and type of input to be used. MACH = 0 indicates 
card input on the IBM 7044 or IBM 7094. MACH = 32 indicates tape 
input on unit 3 (logical unit 54) of the CDC 3200. MACH = 44 
indicates type input on unit 5 (logical unit 1) of the IBM 7044. 
MACH = 94 indicates tape input on Unit B-5 (logical unit 9) on the 
IBM 7094. 

CARD 3 

2 

XSTOP - Number of values in first x table. 


CARD SET 1 

CHILI (m) - 



values for first set of tables (13 values per table). 


CARD 4 

2 

XSTOP — Number of values in second x table. 

CARD SET 2 

x<* u 

CHIUl(m) values for first set of tables (13 values per table). 

v'm 
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CARD 5 

2 

XSTOP - Number of values in third x table. 

CARD SET 3 
X «L 

CHIL2(m) - — values for second set of tables (13 values per table). 

CARD 6 

2 

XSTOP - Number of values in fourth x table. 

CARD SET 4 

X ctu 2 

CHlU2(m) values for second x table (13 values per table). 

CARD 7 

SCALE - Scale value for input data. 

CARD 8 

2 

CONF - Indicates which set of x tables to use. 

CONFl - Percent of confidence within confidence intervals. 

CARD 9 

ALPHA1 — ID card (Normally used for field discription) . 

CARD 10 

ALPHA2 - ID card (Normally used for test discription) . 

CARD 11 

CHAN(I) - Channels to be averaged over and number of channgjLs. For each 
channel N to be used, N must be punched in the n position on 
the card. The number "13" must be punched in all other positions 
except for the thirteenth position which must contain the number 
of channels used . 

DATA 12 

NOTE - If data is on cards they can be followed by the "CHAN (I)" card and 
reloaded for as many channel combinations as required. If data is 
on tape then the "CHAN (I)" cards can be stacked together for 
different channel combinations. 

OUTPUT : 

The program output is given in Figure 3-4 and is self explanatory. 

The input card deck is shown in Figure 3-5. The flowchart of the program 
is given in Figure 3-6, which is followed by a complete listing of the program. 
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. HE L 0 OL SC A I P£1 5 1- _ C5-17 Cj-22 f-S-S C5- 8 

ALTITJOF- 2000 FT. DATE REC0ROED- 6/33/66 



91 PEPCEVT CI'NFUJF'ICE LIMITS 


SYMI^LS* V - VERTICAL I'UJ’IOAftY H - HTRIZvI'ITAL B0UHOARY 
X - 90 TH VERTICAL AND 9<IR I ZO*. TAL BOUNDARY 


S AMPL E NUMOCK 



Figure 3-4. MAP OF INVENTORY BOUNDARIES 
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Figure 3-5. INPUT CARDS FOR INVENTORY BOUNDARY PROGRAM 
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INTEGER Ul“ " 

INTEGER U2 

DIMENSION SYMt90]« ISYM£90]#SYMi(903# JSYM1901 
DIMENSION CHI11E175),CHIU1E175),CHIL2C175 ],CHIU2 1 175 3 
DIMENSION CHAN [ 13 J , I CHAN E 13 3 » DATA U2 ] , V8AR t 12 3 # VBAR1 t 12 1 , VEACUHC12 
1],X ( 3ARE32*90 ]ji_X 8A Rltl 2,903,D EJJK E12,9Q3,XE12 3 # XI C 12J_,X E A C U M { 121 
DIMENSION SUM1 1 90~3 , SUMS £ 90 3 
DIMENSION NV £90 3 
DIMENSION ALPHA1£14}, ALPHA2I143 
DIMENSION IDAT £12] 

NRD-5 

N W T = 6 

READ '['vPD# 6000 ] DOT, DASH, EYE.PLUS, BLANK 
C SET MACH TO 44 FOR IBM 7044 

C SET MACH TO 32 FO'R CDC 3200 

READ C R D * 5 1 } Ll>Ul,i.2#U2 i IAVG/MACH. 

READ £ . «RD# 20 03 3 XSTOP 

MSTOPl = XSTOP _ _____ 

READ URD,2flij £CH I LI £ NU] , NU*1 , MSTOP1 ] 

READ ENPD,20033 XSTOP 
MSTOP2=XSTOP 

READ E '\ R D , 2 0 1 3 £ CHIU1 1 NUi » NUPl«MSTQP2 1 
READ [\RD,2DU3J XSTOP 

M S TO P 3 * X S T 0 P__ __ 

READ £ivRD# 201I t CHIL2 t NU J , NU*1 , MSTOP3 3 

READ ENRD, 2003] XSTOP . 

MST0P4=XSTQP 

READ [ nRD, 201 J ECHIU2CNU3 < NU»1<MS T0P4 3 
READ [MRU, 2003] SCALE 

READ JNRB, 20Q3J CONF/CONFl _ „ . 

fCONFeCONF 
IC0NF1*C0MF1 
WRITE [NWT #1001 ] 

READ I'\RD,2Q04] E AL.PHA1J 1 3 1 1=1^14 J _ . 

RFAD £.\RD* 2004 3 £ ALPH A2 £1 3 , I»l#14 3 

52 CONTINUE 

READ IKRD, 20033 E CHAN £ I J , J *1, 13 3 
DO 6U 1=1,13 

6-j CONTINUE 

i «t CHANs CHAN E 13 3 

N Fjr_A G s u 

NST0P = 0 

IPRINT = f] 

IF EMACH-323 62,61,62 

61 NUN I T = 3 4 _ _ _ 

REWIND NUf-'IT 

J _ GO TO 65 _ __ 
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66 NUN LT s 9„ 

K £ K 1 N D K > U I T 
GO TO 65 


67 NUMT*5 

6b CONTINUE _ * 

DO 6 J s 1 , 9 0 

SUM1J J) = 0 , 0 __ 

SUM2tJ)=U,0 

NVlJlaU 

6 CONTINUE 

101 CONTINUE . . . __ 

DO 5003 J#l,90 

. ...J5.YH LJ I = C 

5000 SYhl[JJ=DOT 

. VSUM1 = 0,0 . 

VSUM2=U ,0 

K” 0 . . 

N*0 

DO 20 ,3=1,90 _ 

IF EMACP] 68,69,68 

69 READ Ii'.UNIT#10J01.MJ5TQPiLCJBAIA.m t I«l A JL2J. . .. 

GO TO 70 

C68 READ MUMTJ M S T O P5 j. I DAIAJU La. l~.lt 1 21 

68 READ I.NUMTJ MST0P5, NDUMMY# NDUKHY* El DAT III , I *1, 12 J 

DjP _ II. .1 =.l ,_± 2 

71 DATA f IJ s I DAT 1 1 3 

c . no 71 1=1,12 . . _ _ . 

C71 DATA £ 1 1 = IDAT [ 1 3 

7 0 CONTINUE 

N = N+1 

. . .jiuiiimji+i 

K = K + 1 

DO 12 1=1,12 

DATAin=DATA£.Ll^SCAL£ 

IF [MST0P5J 9,9,21 

, 9. _ , 

10 XPARU>K]=DATAm 

XSARlt I,K3=DATAII]._ . _ _ 

DELXU, «] = !). G 

11 CONTINUE . ... . 

12 CONTINUE 

CALL g* RORJ y B a U Ni,XC H A jj^ J 
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IF I >1-1] 14,14,13 

13 C 0 sv* T I * ' G _ 

CALL SlOPE [VEACtlM, I CHAM, M * 1, ERCK, VSUM1, VSUM2 J 
IF [ERCK] 113,14^113 

113 CONTINUE 

_ . CALL n u b E LM J g Hi UIulOHJ U1,CH IL 2, CHI U2 , H'STQPl , M S T0P2,H ST OP 3 ,M T 0 P± L N 

1,1 AVG,tRCK, NBOUNO, I COMF1 

114 IF [a; ROUND) 376,14.,376 . 

376 SYrtltNlIsEYE 

JSYHIwiJa3 . 

N = 0 

. . . vsu^i=u A .a. ______ 

V S U M 2 = 0 , 0 

14* CONTINUE ” 

DO 15 1=1,12 
Xm=XoARU,Kl 
XIII 3=X8AR1[I,K] 

15 XEACuMt n=DELXCI,K) 

CAlL ERROR [X, XI, DATA, XEACUM, I CHAM, MVJK H 
IF MVtKl-lJ 'il‘,18,16 " 

16 CONTINUE 

CALL SLOPE [XEACUM, j CHAN, KIV [K ) , NV [ K 3 , ERCK, SUM1 1 K 1 , SUM2 E X I 3 
IF [ ERCK 3 374,18,374 

374 CONTINUE 

CALL COMFLm[CHI L1,CHIU1,C H1L2,CHI U2,MS TOP1,MSTQP2,HSTOP3,MSTOP4, 

IN V [ K 3 , T A V G , ERCK , NBOUND, ICONFf 

NVlsNVtKI 

375 I F [ N b 0 u N D 3 377,18,377 

377 IF [ ISY W IK } - 3 J 379,378,379 

376 SY^i [K ]sPLUS 

N V [ K 3 = C __ 

GO TO 380 

379 SYM £K ]=DASH 
38 0 M V [ K ] = 0 

SUhltK 3=0,0 
SUM2 [K j =0 , 0 

18 CONTINUE _ _ _ 

OO 17 I =1,12 
X9AR[J,*3=xm 
X K A a l n,<3=xin ) 

C C L> il,* )*XEACJ*»f! I 
L 7 COnTI' 

20 CONTI v lE 

21 K = K-1 * ' " 

IF HFlAS] 31,30,31 

JO CALL OuTl [ I CHAN, ICONF ICONFl, K, ALPHA!, ALPHA2, DASH, EYE, PLUS] 

K1 = K 
NFLAG=1 

DO 25 J=1,K_ 
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SYMt J]=SYM1 EvU 

?5 ISYM[ J]SJSYM[ J] _ 

GO TO 101 
-'1 CONTINUE 
K = K1 

IPRIMT = IPRINT+1 

IF T I PR I NT -211 500,501,502 
'00 WRITE [NwT,6GGl3 IPR I NT# £ SYM l JJ , J*i, K J 
GO TO 520 

•01 WRITE [-VWT, 6002 3 I PR I NT, [SYM [ Jj., J«1 , K L 
GO TO 520 

02 IF IIPKINT-23I 503,5 04,5 05 

~U3 "WRITE" [NUT, 60 031 IPR I NT, t SYM I J J # Jsl, K 1 
GO TO 520 

04 WRITE (NWT, 60043 IPR I NT, t SYM £ JI , j*i , K J 
GO TO 520 

05 IF UPRINT-253 506,507,508 

06 WRITE t NwT, 60 05 J IPRIN T, C SYM [ J3 , J-l, K 3 
GO TO 520 

J7’ GO TO 50 0 

38 IF [ IPnINT-273 509,510,511 

39 WRITE [ ,V WT, 60 06 ] IPR I NT , I SYM [ J 3 , J = 1 , K 3 
GO TO 520 

10 WRITE Lf'WT, 6007 3 IPRINT, ISYMIJ3 , 

GO TO '520___ „ 

'll IF II PRINT-29] 512, "513, 514 

12 WRITE [NWT, 60083 IPRJNT A tSYM U 3 , J=l, K J 
GO TO 520 

13 WRITE ENWT, 60 09 3 IPR I NT, [ SYM [ J3 , J=l, K3. 

GO TO 520 

14 GO TO 50 0_ 

20 DO 521 J=1,K 

SYM JI=SYM1EJ3 __ 

21 "iSYM JJaJSYMl J] 

K1 = K 

IF CMSTOP5-23 101,101,53 

3_ WRITE [NWT, 10011 

GO TO 52 

1 F0RMATI8I3J _ __ 

01 FORMAT L13F6, 0,2X3 

"000 FORMATII1,8X,12F5»OJ 

001 FPRNAT13H1) 

J.Q08 FOR V AT[13F 5,OI 

>004 FORMAT [13A6,A2] 

5000 F0RMAT14(9X#A13,5X,A5JL _ 

50 01 FORMATIlUX, I4,1X,90A13 

50 0 2 FORMATE 9X, IMS, I 4, IX, 90 A1 3 

5003 F0RMATC9X,1HC, I4,1X,90A13 

6004 F 0 R M A T I 9_X i 1H A , I 4 , IX , 90A1 3 , 

6005 FORMAT [9X,1HN, I4,1X,90A13 
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6006 F0RMATC9X.1HL, I4 ,1X,9QA 11 

6007 F0RMAT[9X,1HI, I4,1X,90A1 J 

6008 FORNAT [9X,1HN, I4,JJC,90A1J _ 

6009 FORMAT I9X,1HE, I4#1X#90A1) 
END 



SUBROUTINE OUT1 1 ICHAN# JCONF", ICO Mffl,N,At> h A 1 , A L PH A 2 ; D A S H E Yc, PLUS J 

DIMENSION ITEM [90 3 * IONEC90J 

DIMENSION f CHAiv [13,1 CNAN1 1 12 J 

COMMON/C ATA/WlNGTHt 13 3 

DIMENSION WLNGTH [133 

DIMENSION ALPHAltl 3 ,ALPHA2(1J 

DATA It MlNGTH fl 3 , I" SJ , 13 3 = . 40 , .~44 , . 46 , , 487. 5 0 7. 527 . 55 , 38, 62 766 , 
1 72. .60,1,003 

DATA PLNGTH/,4# .44# .46# .48# .50# ,52#.55# .58# .62# .66, . 72# . 8, 1, / 

NWT = 6 

DO 9 1*1,12' ' 

9 ICHAN.lt 1 3 = 0 
"1 1*1 

Jl*ICH»v[l3J 
DO 12 u=l,Jl 
DO 10 1=11,12 
12=1 

IF II chan UJ. " 13 3 11 , 10,11 

10 CONTINUE 

11 ICHAnI t J3 = ICHANtJ23 
11=12+1 

12 CONTINUE 

WRlfE INHT# 10013 

WRITE INwT, 10103 _ ; 

WRITE I'HT# 1007 3 I ALPHA1 1 1 3 # 1*1,14} 

WRITE [MT, 1007 3 I ALPHA2 Cl 3 #J*1 # 14 3 _ 

WRITE [MWT# 10103 
WRITE INWT,1Q06J 
DO 13 1=1, J1 

.M=ICHA^im 

13 WRITE t hkt i i'0 0 8" 3 M , WLNGTH C M j”WLNGTH Hm+3. 3 ' 

WRITE i'wT, 10103 

aRite r-KT.iensj jcg'nfi 

« F I ~E i 1 T » 1 &1 D I 

aRITE C'wT. 1009 1 EYE, DASH, PLUS 

aRITE t^T, 10103 

WRITE £^Hf#lQ02j " ” . - - ' 

WRITE [ ^ wT, 1G10 3 
CO 20 1*1# M 
ITE\>{ n = I /10 
I one t i 3 = i - [ i /i 0 1 *1 o~ 

20 CONTINUE 

WRITE lNwfil6b3T"nTEN[Il# 1 = 10, Nj" ‘ “ 

WRITE ENWT, 1004 3 [ I ONE [13 # I *1, N J 
1001 FORM ATI 1HG,28X#27H MAP OF INVENTORY BOUNDARIES/] 
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1002 FORMAT [36X#13HSAMPLE NUMBER/ 1 

1003 FORMATE24X,90il) " 

1 004 FORMAT U 5X*90J1/1 

1005 FORK AT 1 2 8 X * 1 2 , 2 6 H PERCENT CONFIDENCE LIMITS/'] ~~ 

1006 FORMAT [26X,32HCHANNEL WAVE LENGTH (MICRONS J / ) 

1007 FORMAT U4X13A6,A2/1 

1006 FORMATE2SX, I2,12X,F5.2,3H - ,F5,2,/J 

1009 FORMAT114X,8HSYM80LS=»3X,A1,20H - VERTICAL BOUNDARY, 3X, A1,2?H - HO 

3 R 1 20 a-T A t BO'JNDaRY/? 2X#A1,40H - BOTH VERTICAL AhP HORIZONTAL 80UNDA 
2RY/3 ' ' ' ‘ 

1010 FORMAT I1H0] 

RF7URN 

END 


1 

SO 6'R o U'T Tn r~E'R RWl x m OTTx mCTRT# xTiaTAVxEaCum TTCPa n,' »-'] “* 

DIMENSION XMEANE123 , X^EANIJ 12] , XEACUM E 12 3 , XDATA E 1? 3 , I CHAM 11 

XK=M-1 ~ 

DO 30 3=1,12 

IF IICHANI 1 1-131' 20,30,30 
?0 IF CM-i] 21,21,22 

21 XHEAnI ffaYDATAfn 
XMEANlt I JsXMEAMI 3 
XE ACUM f 1 3 — 0 * U 

GO TO 30 

22 CONTINUE 

XMEAMI ]9[XM/tXM*1.0n»XHEANin»[i- t 0/(XM*l,Dn»XDATAI I 3 
XE AC UM 1 1 3 = I E XR*lTT3 /[ XM + 1 « G7 3 *XEA5 uM ( T 1 1 •0/1 XM + 1 ,0 3 3 * E XMEAN1 113“ 
IXMEANEII 3**2+El,0/[XM*tXM+l. 0 11 1 *£XDATA UJ -XMEAN C I] 3**2 
XMEANl E I 3 =XMEAN [ 1 3 
30 CONTINUE 
RETURN 

END _ _ 


SUBROUTINE SLOPEIEA CUM, I CMAN,M,M0UM,E RCK,SU M l , S U M 2 3 

D I MEN S 1 6Tv“ EACUN [ 1 3 # I CHAN [11 

CHAN13- I CHAN { 13 ] 

XM=M 

HACUH1=0,0 

DO 20 1=1,12 

IF Cl CHANJ I } ~13 3 10,20, 1 0 „ 

10 H ACUM1 sHACUMl + E ACUM [ 1 3 

20 CONTINUE _ _ __ 

HACUMl=$uRT [HACUH1/CHAN13 J 

SUM1 =SUM1 ♦EXM»HACUMll»*2/rxH«l,03 

SUM 2 = SUM2 ' ♦CXM*HACUMli/SORT(XM-i,0J 

JT tSJJMU 30,4 0,30 
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30 IF CSUM23 35,40,35 
40 EPCK=0.0 
GO TO 45 
35 CONTIrluF 

A = S U m 1 /SUM 2 

ERCK = S:>'RT t XM ] *H ACUM1/A 
45 " CONTI HUE " 

RETURN 

END 


" ' , 5UBffijimw _ TnnJFT^^ pstd ps; fstdpt7 

1HST0P4,M, IAVG,ERCK,NBOUND, IC0NF3 
~D f HENS I ON CHILI [ lT , CH I'Ul [ IT', CH I CTO J , CH I U2 [ i J 
IAVG1=IAVG-1 

NSOUND=0 

IF [M~MSTCP3'I 342,342,341 . 

342 CHi0l2="cHir2THl " 

GO TO 343 

341 CMI0L2S0HIL2EMST0P3I 

343 IF [M -s'* ST0P4J345, 345,344 

345 CHIOlJ2 = CHIU2 [MJ‘ . 

GO TO 346 _____ 

344' CH 1 0U2 = CH I U2 fMST&P4 1 ' “ ' " 

346 CONTINUE 

IF [M-NST0P2I 31,31,30 

30 CHl0ul*CHIUlfMSTQP2] 

GO TO 32* " * ’ 

31 CHIOUlsCHIUl [HI 

32 IF £ M-NSTOPlT 34,34,33 

33 CHIOLlaCHILlCMSTOPU 
GO TO 35 

34 CHIOLl*ChILltMI 

35 IFIM-IAVG1J 20,37,37 ' 

37 CONTINUE 

> "IF£IC0**F-1J 372,371,372 ~ ~ 

371 IF EtRCK-CHIOLl J 376,376,36 
36 IF I&RCK-CHIOU11 20,*376,376 

372 IF[ICO\F-2J 376,373,376 

373 IF [ERCK-CHIOL23 376,376,374 

374 I£ 1ERCK-CHI0U2J 20,376,376 

376 NBOUwD#l'“ " " 

20 CONTINUE 

7000 FORMAT C60X,2F10 ,51 
RETURN 

END ' 


I 
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3.5 BOUNDARY ENHANCEMENT PROGRAM 

This program is devised for determination of inventory boundaries through 
the enhancement of the difference in the multispectral data between homogeneous 
inventory areas. The mathematical formulation has been given in subsection 2.4. 

This program consists of two links. Link 1 performs the moving or jumping 
averages of the input data. The output is in tape and may be input to Link 2. 
Link 2 performs the calculation of the finite difference of the averaged data 
and then calculates the absolute value or square of these differences. The 
outputs of Link 2 are (a) means and standard deviations for every channel of 
enhanced data and (b) a digital tape of the enhanced data. These means and 
standard deviations can be used for setting the class intervals for inputting 
to the Alphanumeric Plot Program to process the digital tape. 

The input data to this program is on tape, while the control parameters 
for the program are inputed by NAMELIST statements which are documented in IBM 
FORTRAN IV manual, Form C28-6390-2. 

The definitions of the input parameters to each LINK are as follows: 

LINKO INPUT PARAMETERS : 

i 

LINKNO - Specifies which LINKS to execute. LINKNO equal one LINK1 is executed 
and LINKNO equal two LINK2 is executed. 

LINK1 INPUT PARAMETERS : 

e NLOOP - Number of scan lines to be read in. 

• NCT - Number of samples per scan line. 

• NCHAN - Number of channels of data to be used. 

• LTN - Logical tape number of input data tape. 

• LTNl - Logical tape number of output data tape. 

a NLOOP1 - Number of scan lines to sort, dependent on the number of scan 
lines read (NLOOPl = NLOOP-NSKIPR) . 

a NSKIPR - Number of scan lines to skip, cannot be greater than the 
number averaged (NAVGR) . 

• NAVGC - Number of samples to average in a scan line. 

• NSKIPC - Number of samples to skip in a scan line, cannot be greater 

than the number averaged (NAVGC) . 
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• NAVGR - Number of scan lines to average. 

• NFILE - Number of data files to process. 

« NBEGIN - Scan line to start processing. 

LINK2 INPUT PARAMETERS : 

• NLOOP - Number of scan lines to be read in, must be greater than or 

equal to (2 X NDXFF + 1). 

o NCT - Number of samples per scan line. 
a NCHAN - Number of channels of data. 

• LTN1 - Logical tape number of input data tape. 

• LTN2 - Logical tape number of output data tape. 

« NLOOP 1 - Equals NLOOP minus one. 

• NDIFF - Number of adjacent samples to skip in the cross differencing 

process . 

• NFILE - Number of files of data to process. 

• LABS - Calculates the sum of the absolute values of the differences 

about the adjacent data points. If IABS is not equal zero, it 
is calculated J if it equals zero, it is bypassed. 

• ISQ - Calculates the sum of the squared values of the differences 

about the adjacent data points. If ISQ is not equal zero, it 
is calculated; if ISQ is equal zero, it is bypassed. 

• I SUM - Position in the output array (W(I)) that the sum of the squared 

or absolute value of the differences from the channels is stored. 

• NBEGIN - Scan line to start processing. 

An example of the output is shown in Figure 3-7 , which is followed by a 
complete program listing. 

3.6 MEAN, CORRELATION, AND COVARIANCE PROGRAM 

This program is devised for calculating the mean vector, correlation 
matrix, covariance matrix, and normalized covariance matrix of a given multi- 
spectral data sequence. The format of the input data is specified in the Tape 
Conversion Program. The mathematical expression for these quantities is given 
in subsection 2.6. 

In addition to the printout, a punched card output for the mean vector 
and covariance matrix are available. These card outputs form part of the 
input to the Divergence Matrix Program. 
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INVENT0RY 80UNDARY MAP_8Y BOUNDARY ENHANCEMENT 
CROPS C5-17, 5-22, 5-5, 5-8, ALTITUDE 200C FT. 


DATE RECORDED 6/30/66 


SYMBOLS 


CHANNEL 

13 

= LESS THAN HE5TT 


HAVE LENGTH 


'HEAN-MEAK+1SD 


HEAN*IS0-MTEAN+2S0 


MbAN+2SU-HFA'N43SU 


“T~= — CTFSTFR THAN Hfc AN+3S0 


DOOOOOOOOOOOOOOOOOODOOOOCOOOOOOOOOOOOOOCOOOOOODOOllIllllHimilllUlIIIUllllU 
000011 11122222333334466455555666667777768888999990000011111222223333346664555556 
246802 468O246BOZ468O2460O246BOZ46ffO2468O2468O2468O24'&8OZ468O2468OZ468O2468o2468O 


2 • 



4 


:: * 

6 *•• • . 



8 *. . , . 


* 

i a _ , 



TZ~*; • 





14 



16 

- - 


18 


UIHIMimilllilllfl 

20 *«••• 

. , f f 


22 



— — 22?t>23?2 
26 *2322 

28 *. . . . 



.....2332S...-33332...-. .6 6**62«.. 


: n *rm2\ ' ' 


30 *.... 


- * w 

32 *•••. 

34 * 


•" M 




38 • 



40 + 



* * ##il 



44 ***** 

* * .tt - i ; * 


46 ***** 






48 *• •• • 

50 



52 ***** 



54 



56 *♦.*. 



58 *.... 



6? *- 



64 

- - - -- 


66 *- , . m 



68 - 



70 *.... 



rtrtrzr: 

74 ***** 



76 *. * . * 



78 *22.. 
80 *.... 
82 *.... 



_ 7;7T~ Y~~7TVfVymY2-^ 7XZ T333 7277 771^113 377.T.7I 

™’ , .T.M3132 



84 

86 *.... 


• • • • • 6*2* 3* •» 

88 *.... 
90 *.... 
"92 *.... 


• ••••••2* 32 *» •••«•• 

• •• 3« 33 


94 

96 

98 

10D 

102 


.3.33. 

.3.33. 

.3.33. 

,3.33. 

.3.33. 


Figure 3-7. INVENTORY BOUNDARY MAP BY BOUNDARY ENHANCEMENT 
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C BOUNDARY ENCHANCEMENT PROGRAM 

dimension Link no c 6 3 

N A MEL I S T/NAMDAT/LIN K NO 

READ t 5, NAMDAT J ' 

WRJTE 16, NAMDAT 3 
IFtLINKNOIll . EQ , 0 j GO TO 1 
CALL LINKl _ 

1 CONTINUE 

IF[LIN*N0t2 ] .EQ.QI GO TO 2_ 

CALL LINK2 

2 CONTINUE 
STOP 

END 

$OR I G IN ALPHA, REW 


SLBFTC TWO NQDECK 

SUBROUTINE LINK1 


C 


MOVING 

DIMENSION 

DIMENSION 

NAMELIST/S 


OR_ JUMPJLNG AVERAGE PROGRAM 
DUMMY l 13] 

DATA_[6_,20Qj 12], XI 200*12 J 

ZE/NLOQP, NCT , NCHAN»~LTN, LTN1,LTN2, NLOOP1, NSK IPR, NAVGC, 


inskipc,na vgr, n file,nbegi n _ 

NSK I PR MUST BE LESS THAN OP EQUAL TO NAVGR 

NSK I PC MUST BE LESS THAN OP EQUAL TO NAVGC 

READ [5 , S I ZE 3 

WR ITE 1 6, SI 2JE 1 . 

DO 101 I A8=l, NF ILE 

. ,IZER0 = G „ 

IONE = l 


I TW0=;2 
ITHREE-3 

KS = 0 

XNAVGC=N AVGC 


. . .mkVQHzU 1£R 

100 CONTINUE 

LL = 1 . 

KK=NAVGC 

DO 10 I =NBEG_I.N ttiLOOR 
KSsKS+l 

DO 11 J»1. , NQT 

RFADtLTN] ICHEK, IDUM, [DATA! I , J, L J # L s l* NCHAN ] 

C WRITE 16,10031^. JCH.EK,KS, J , l DATA U ,.J, L ) , L*1 , NCHaN ] 

iniCHfcK-3] 11,7,11 

7 inj.Nt.NCT3 GO TQ 13 . ... 

11 CONTINUE 

UQ CONT I.N.U.E . 

DO 12 1=1, NCT 

DO 12 J=l, NCHAN 
XU,J] = 0.0 

12 CONTINUE 
NCT1=NCT-NAVGC 

ao i .ItIolMCii.... 
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DO 3 J=1,NAVGR 

DO 3 K=LL,KK 

DO 4 M=l, NCHAN 
Xt I,M3=X[!, M3+DATA[J ,K>M3 

4 CONTINUE 

3 CONTINUE _ 

DO 2 M=l, NCHAN 

X II » M J =X [ Ij Ml_/i_XNA.V.GC*XNAVQR 3 
2 CONTINUE 

WRITE tLTNU JLZER.Qj.IX I W M 1 > M = l, NQHANJ. . . 

C WRITE [6# 1004 J IZERO, [XU»M],M=1,NCHANJ 

LL=LU*fi§JULEC _ 

KK=LL+nAVGG-1 

1 CONTINUE 

I F £ I CHEK-3 J 9,8,9 

9 WRITE [UJN11_ I ONE, [ DUMMY 113 , I =1 , NCHAN 3 . 

IFIMSKIPR.EO.NLOOP] GO TO 100 

. _ ■NBEGIN = NL0^QP1±1 - -- 

DO 5 I=l,NLOOPl 

DO. 6 J = L,N£I 

DO 6 1=1, NC HAN 

LZsI*NSKIPR . 

DATA £ I » J,LJ=DATAILZ,J,L3 

. 6 . . GGNLLME 

5 CONTINUE 

GO TO l .00. 

13 WRITE [6,10023 KS 

8 WRITE E.LTNU. - I TWO, [DUMMY [ I J , 1=1, N.CHANJ. „ 

101 CONTINUE 

. . WRJ T £ [ L mil JLB.UHB Y 1 1 h I = lx MfiHAM.1 - - - - . - 

END FILE LTN1 

REWIND ..LTN1 

1000 FORMAT! II, 8X,12F6, 33 

10 01. FORMAT [I1 j lQXx1£L6x31 

1002 F0RMAT[1X,11HEND FILE IN, 15, 9HSCAN LINE 3 

JL 0 0.3__£Q.R M.AI11 

1004 FORMAT [IX, I3,12F6.33 

RETURN .. _ _ 

END 

SOR.I G I-N _ ,AL£HAjJ£W 

SIBFTC THREE NODECK 

C* CROSS DIFERSNC? PROGRAM™''’ 

DIMENSION DUMMY £333 

DIMENSION DATAE6,200,123 ,Wtl5J ,Wl[153 

D I MENS 1 ON X8AR C 15 3 , XMSBARi 15 S1G UJJ , RMS t 15 3 , S I GSQ 115 J 
NAMELIST/S I ZE/NLOOP, NCT, NCHAN, L TNI, LTN2, NLOOP1, ND IFF ,NF ILE, 

1IABS, ISO, IS. UM , NBEGIN _ 

RE AD [ 5 , S I ZE 3 
WRITE [ 6 , S l Z EJ 
NCHANlsNCHAN+1 

DO 102 IAB = l,NFJki „ . 

IZERO=0 

r " TW0 ~: 2 
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I THREE 8*3" “ 

..MS SO 

KS = 0 

_DO 18 I =1 » NCHAN1 
S IG I IJsO.O 

XBAR [ H sQ» 0 

XMSBAR £ 1 3 = 0 » 0 

RHSmgQtO 

siGsomso.o 

..18 CONT I NUE 

100 I F [K$ , GT . 0 J NBEGINsNLOOP 

DO 1 I g NBEG I N, NLOOP 

KSsKS+1 

DO 2 J »1,NCT _ 

READ ILTNl J I CHEK, IDATAIJ#J,LI ,L*1« NCHANJ 

C WRITE [ 6,10033 ICHEK ,KS, J , IDATAI I, J, LI ,L» 1,NCHAN3 

IF 1 1 CHEK-3 } 2 * 4'* 2 

4 IFtJ.NE.NCT3 GO TO 7 

2 CONTINUE 

_i CONTINUE 

NCTlsNCT-NDIFF-1 

ND IFFl g NDIFF»l 

NL00P2=NL00P-NDIFF 

DO 10 I g ND I FF1 , NL0QP2 - . 

DO 11 U*N0IFF1,NCT1 

SMtfMS 

MSsMS+1 

WJ I SUM 3 sQ , Q 

DO 12 Ksl,NCHAN 

-12.s I "NO IFF 

JZ = I *ND I FF 

KZs J+ND IFF 

LZsJ-NOIFF 

_C. ABSOULTE VALUE OF CROSS DjFFERENC E S 

IFUABS) 9,8,9 

jz iimimuE 

AVSUMsABSEDATAtl, J,K]-DATAf IZ, J, K J I + ABS t DATA 1 1, J, K J-DATA [I ,KZ,K 33* 

1.ABS [DA TA I j-<_ J , K 1 _« P A TAX J Z ,J,K 1 1+ AB S tDATAt l, J,K)-DATA 1 1 ,LZ ,K 1 1 

WrKlsAVSUM 

_8 CONTINUE 

C SQUARED VALUE OF CROSS DIFFERENCES 



14- 


13 

12 


IFCIS Ol 14,13,14 
CONTINUE 


1**2*IDATA[I , J,K)-DA TAU, KZ»K11»«2 
l+IDATAn,J,K}-DATAtJZ,J,Kn**2*IDA'TAn,J,K}-0ATA( WU#K]3**2 
W IKIsAVSLJMl 


CONTINUE 

-Mix s U_MJ_s.M.US.O, M JL±wm. 


CONTINUE 


„„DO .15 L g l, NCHAN1 

XMSBAR IUs ISM/ 1 SM+1,0] J *XMS8AR IL 1+U, 0/ tSM+1 . 0 13 *W (U **2 
XBAR a 3 s I.S.M./ I^M±loJ).n*XBAJ5 U.I+ 1 1 . 0. / .1 S Htl * 0 3 I*WJLJ 


15 CONTINUE 

WRITE.LLXNZJ lZ£&fl.t.iJiIjllj. K« lr-NCHANl 3 . 

C WRITE 16,10011 I ZERO, IWIK1 , K*l> NCHAN J 
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C WRITE “[671ml ~TZER0* lWl[K]7k = l/NCHANJ 

11 CONTINUE _ _ _ _ _ 

IFCICHEK-3J 1 D / 3 / 10 

10 WRITE t UTN2J... imeWBUMMYCKI ,Ksl/NGHANlJ 
DO 5 I=l,NLOGPl 

QD 6 J=ltNJGI - - 

DO 6 L=1,NCHAN ^ 

DATAtI#J,U = QAIAJl+l# J./U . 

6 CONTINUE 

5 CONTINUE . ... 

GO TO 100 

7 . ttRII6I$ x liLB21__-KS 

,3 WR ITEILTN2 3 I TWO/ t DUMMY [II / 1 = 1/ NCHAN1 1 

101 CONTINUE 

DO 16 I*l» NCHAN1 

..5IGlII?JU4SJBjABIII*JtEAR[t.l**2 

SIGSQU ]=SQRTf SIGt I ] I 

..... RMS { IlES5R.LU]l^AH_LIJ_l. . _ . . . 

16 CONTINUE 

WRITE (6^100.41 

DO 17 L s l» NCHAN1 

WRITE (6il 0 0.5.3. _li_XBAR (L 1/ XMSBAR CU «S I G ( L 3 /S IGSQ ( LI 

17 CONTINUE 

.. 1.0,2. . CO.NUNgE _ 

WRITE ELTN2 3 I THREE / [DUMMY [ 1 3 * I *1/ NCHANl ) 

END FILE . LTN2_ 

REWIND UTN2 

1 0 0 0 F0RMATJIli8XjLl2Efit-31 . 

1001 F0RMATtU/8X/12F6,31 

U) 0.2 Fp R M A T. [ 1 X / 1 1ME.MD-F J LE_I.N / i 5 / 9 jd S C AN _ LINE 3 . . 

1003 F0RMAT(1X/3I5/12F6,33 

10 0 A FORMAT (iaX*^.yCJii 57fi.4JdllE.AN/ 6X./11HM.ENN-SOU ARg/6X/8h VARIANCE /6X/ 

17HSTD*DEV J 

10 0 5 FOR M AT. ( 10 X ti 3 jl£AJ&j 7 /5 X/F 1D » 7 / 3 X / F_l_Q.t 7 /. 3X / F1_0 1 7 , 3 X , F 1 0 , 7JL . . 
RETURN 

£HD 

SDATA 

SNAMDAT 

LINKNO-0/2/ 4*0 

SEND . . „ 

SSIZE 

.NL0.OP=3 - 

NCT=164 

N C W A N = 1 2 . „ . 

L T N 1 s 9 

L T N 2 s 1 1 . . 

NU00P1-2 

NDJ FF S 1 . ... _ 

NFILE=1 

nbegin=i 

IABSpO 

ISQ*1 

ISUM=13 

SEND . _ . .. _ . 

SSIZE 
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Ml OOP s 5 
NCT = 1 68 
MCHAN = j 2 
LT^ = 9 
LTM1=10 
LTN2*U 

ML00 D 1=4 

NSKlPRsl 

NAVuC*5 
MSKIPC»1 
MB 8 'JC = 1 

NEN 0 Ca 5 
WAVUR=5 
MBP(aIM = l 
MFILFsl 
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The input card setup is shown in Figure 3-8. 

DEFINITIONS OF INPUT PARAMETERS : 

CARD 1 

• NSCAN - The number of scan lines in a field 

• NCHAN - Total number of channels in a field. 

• NCOMP - Number of channel cross products calculated to generate the 

Covariance Matrix. 

• NTOTAL - Total number of samples in a scan line. 

• NFIELD - Number of fields to process in a data set. 

• ISCAN - The number of scan lines to skip before starting calculations 

in a field. 

• NFLD1 - Number of times to process the data set (N X NFIELD) 

CARD 2 

• An 80 Holerith card used for field identification. 

CARD 3 

• NBEGIN - Sample number to start with in a specific scan line for a 

certain field. 

• NEND - Sample number to stop on in a specific scan line for a certain 

field. 

CARD 4 

• An 80 Holerith card used for field Identification. 

A listing of the output from this program is given in Figure 3-9. A 
flowchart of this program is given in Figure 3-10. 

The program listing is given below. 
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C( 5-17) 

SCAN LINE 

SECTION A 
START 

ST ZP 

1 

t 

56 

2 

1 

56 

3 

1 

54 

A 

1 

54 

5 

1 

54 

6 

1 

54 

7 

1 

54 

8 

1 

54 

9 

1 

54 

10 

1 

54 

11 

1 

54 

l? 

l 

54 

n 

1 

54 

14 

1 

54 

15 

1 

54 

16 

1 

58 

17 

l. 

56 

18 

1 

56 

19 

1 

54 

20 

1 

54 

21 

1 

54 

22 


" ' 54 

28“ — 

1 

5 2 

24 

1 

52 

25 

1 

52 

26 

1 

52 

27 

i 

52 

! Figure 3-9. 

OUTPUTS OF THE PROGRAM 
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C ( 5-1 7 ) 


SFCTl ?N 


1.13450 


MEAN VECT2R 

~~2 V " -V 

1.6477 0 L .14.21 6 1.C210 6 . 

CBRHFI.AT12N MATRIX 


1.06431 


1.29448' 
1.37150 
1.29375 
1.16064 
1.20975 
1.18529 , 
0.91027 
1.14810 
0.5269? 
0.71737 
1.0613* 
0.82111 


0.90739 

0.00219 

0.00297 

0.00225 

0.00229 

0.00198 

0.00147 

0.00207 

0.00092 

0.00066 

-0.00034 

-0.00036 


2 

1.871 50 
2.727 80 
1.88613 ' 
1.63567 
1.76 709 
1.72116 
1.32215" 
1.667 16 

0. 76575 
1.04181 
1.5*035 

1. 19234 


1.29875 

U 3_3 6 1 3 

1.31314 

1.16979 

1.21775 

1.19 354 

0.91642" 

1.15721 

0.53188 

0.72263 

1 .06552 

0.92490 


_C 3 V AR I A'lCE 


2 — 

0.00219 
ftT0l290 
0.00420 
0.00348 
0.00343 
0.00243 
0.00224 
0.00271 
0.00180 
0.00090 
-0.00159 
-0. 00073 


0.00297 

0.0042 0 

0.00°61 

0.0035 8 

0.00213 

0.00214 

0.00148 

0.00344 

0.00213 

0.00113 

-0.00332 

- 0.00212 


4 

1. 16064 
1.68 587 
1.16975 
1.C4972 
1.C3374 
1.C6632 

0. 81945" 

1. C 348 5 
0.4 762'? 
0.64605 
0.95135 
0.7367 7 

MATRI X 

4 

0.C022E 
Q.C0346 
0.C0358 
0.C071 6 
"0.00201" 
0. CO 125 
O.C0156 
0.C0341 
O.C0282 
0.C01Q5 
■(T.C0416 
-0.CQ256 


1.2C975 

1-25709 

1.21775 

1.08874 

I.L4035 

1.1116 9 

0.85564 

1.07963 

0.49579 

0.67382 

0.99366 

0.76941 


0.00229 
C . 00343 
0.00213 
0.00201 
0.00809 
0.00149 
0.00306 
0.00455 
0.00233 
0.00145 
-C. 00233 
-3.00124 


’iORMAUZEO C9VARI ANCE “ATR1X 


2 

0.22442' 
1.00000 
0. 39844 
0. 36209 
"0.'33582 
0.25137 
0726 7 3"3 
0.23444 
0. 19443 
0.12931. 
-0. 10191 
-0.05864 


3 

6.37212" 
0.39E44 
l.OOCOO 
0.45614 
’0.25548 
0.27C56 
0.21556 
0.16411 
0.30744 
0.19963 
-6 • 2 62 4" 7 
-0.20957 


4 

" O'. 30 954 
0.36205 

0. 45614 

1 . COOOC 
'“6'. 26425 

0. 1 730 6 
" "6.24922 
0.39575 
0.4079E 
0.20245 
'-0.36066" 
-0.27787 


0.2961?' 
0.33532 
0.25549 
0.26 4.2 *> 

"l I OOOOO" 
0.19454 
0.460"18 
0.49623 
C. 31687 
0.26472 
-0.19027' 
-0.12606 


Figure 3-9.' OUTPUT 


1 


0 


0.80106 __ 


10 11 

4 _CL. 93581 


.0 ,.72.4.0.8 


6 

1.13529 
1.72 116 
l. 19354 
1.0663? 
1.11169 
1 . 0953 1 
6.83689 
1.C5501 
0.48402 
0.65967 
0.97782 
0. 7561 3 


0.91027 

1. 32215 

0.91642 

0.81949 

0.85564 

0.83689 

0.64716 

0.8 1 229 

0.372B9 

0.50695 

0.74945 

0.58046 


8 

1.14810 

1.66716 

1.15721 

1.03485 

1.07968 

1.05501 

0.81229 

',0.47294 

0.63976 

0.93955 

0.72739 


C. 52692 
C. 76575 
0.53188 
C. 47623 
C. 49579 
C. 48402 
0.37289 
C. 47794 
C. 22164 
0.29306 
0.42837 
0.33228 


1JL .. 

0.71737 
_ 1.04181 
0.72268 
0.64609 
0-67382 
0.65967 
0. 50695 

0.29306 

0.40282 

0.59004 

0.45651 


. LI . 

1.06134 
1.J4035 
1.06552 
0.95135 
0.99366 
.. 0,9.77.82, 
0.74945 
.9,939.55 
0.42837 
0.59004 
0.89433 
.0.68781 


. .12 J 

0.82111 ’ 
1-19234 
0.82490 I 
0.73 677! 
0.75 941 | 
.0-75-6L3, 
0.58046 | 
0,7273<tj 
0.33 228 ! 
0.45 651 ! 
0.68781 ' 
JJ.53618 | 


6 

0.00188 
0.00243 
0.00214 
0.00125 
0.00149 
0.00723 
0.00130 
0. CO 129 
0.00039 
0.C0070 
0.00167 
0.00081 


0.00147 

0.00224 

0.00149 

0.00156 

0.00306 

0.00130 

0.00546 

0.00309 

0.00149 

0.00039 

-0.00018 

0.00042 


8 

0.00207 

0.00271 

0.00344 

0.00341 

0.00455 

0.00129 

0.00309 

"0.01038 

0.00459 

-0.00577 

-0.00405 


9 

C.00C92 

0.00180 

0.00233 

0.0028? 

C. 00233 
0.00039 
0.00149 
C. 00459 
~C70T 6 6 8 
C. 00015 
-C. 00551 
-0.00343 


10 

0.00066 
0.0 0090 
0.00U3 
. .CL 00 108 
0.00145 
0.00070 
0.00089 
0.0016C 
0-00015 
_0,.OO373. 
-0.00115 
-0.0009? 


-0.00034 
. — 0.0 0158 
-0.00332 
-0.-00*16. 
-0.00233" 
0.00167 
-0.00018 
--0,00577 
—0.00551 
. tlO, 0,01 L3 
0.01859 
0.01020 


-0.00 036 | 
-o.OMZji 
- 0 . 00212 ' 
— 0. 002S6 I 
-0.00124 ! 
0.00 OB 3 \ 
0.00042 I 

- 0 , 02403 ; 

-0.00343 , 
r0,0Q092j 
0.01020 
0.01 188 . 


0.23174 

0.26711 

0.21556 

0.2492? 

0.46018 

0.20619 

l.OOOOO 

0.41050 

0.24678 

0.19758 

-0.01817" 

0.05237 


8 

o'. 2 161 2 

0.23444 

0.36411 

0. 39575 

0.49623 

_0.J*922 

6.41050 

1.00000 

0.55093 

0.25 650 

'6.41557 

-0.36464 


9 

C. 13042 
0.19443 
C. 30744 
C. 40798 
C. 31687 
C. 0 5594 
C. 246 78 
C. 55093 
1.00000 
..C -03 C99 
-C. 49461 
-C. 38526 


10 

11 . 

12. 

0.12555 

-0.02887 

-0.03883 

0.12931 

-0.10191 

-0.05864 

0. 19968 

-0.26247 

-0.20957 

0.20245 

-0.36068 

-0.27787 

0.26472 

-0.19027 

-0.12606 

0.13470 

0.14408 

0.06981 

0.19756 

-0.01837 

0.05237 

0.2565C 

-0.41557 

-0.36464 

0.03099 

-0.49461 

-0.38 526 

l.OOOOO 

-0.13767 

-0,12312 

-0. 13 _ 767 

1.00000 

0.68642 

.-0.13817 

0.68642 

1.00000 


OF THE PROGRAM (ConcTuded) 
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Figure 3-10. CORRELATION AND COVARIANCE MATRIX PROGRAM 
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C 


C 

C 

C 

C 

C 


34 

2 

71 

3? 


C 


33 


53 

57 


l 


CORK ELAT I?N AND COVARIANCE MATRIX PROGRAM 
DIMERS IPN X XX XX { 1 2 » 12 I 


DIMENSION XMEAht 12) » DATA-{ 1-2 ) » C OPR ( 12*12 ) 

_DJ_M F NS I HN AC2RR ( 12,12) t ACOVR ( 1 2 , 12 ) t XR A R ( 1 2 ) 
LTff logical tape number cf data tape 
LTM= ) 

LTN5 LOGICAL 
LTN5-3 

LT!J6 LOGICAL 
LTN6=6 


NUMBRR 3F INPUT PARAMETERS 
NUMBER OF OUTPUT 



LTN 7 LOGICAL NUMBER OF PUNCH OUTPUT 

LTN7- 7 

'INPUT PARAMETERS ""iN'ITIALfZED 


KS = 0 

SCALE = l CO. "O' " 

IB iE=l 

I FLD 1 = 1 ' ' 

IFIbLO=i _ 

IF (KSf 71,'?;71 
CTI-TI'Iije 

READ' (LTN5", 100) ' NS CAN ,NCHAN, NCEMP", NTOT A L , NF I ELC , I SC AnTKfUM 

C_' I TINUL ^ 

IF( If : IPLC-NFI6LD)32,32,33 
C'N'jTl jUE 

MS-0 ’ ‘ " “ - ” 

"PTI IN to skip scan likes 

I F ( IS CAN) 57, 57, 5 B 

C*MTl ,IIF 

C' 1 5 3 T = l » I SC AN 

DO 55 J=1»N TOTAL 

READ ( LTN") ’ IDUM' ” . 

CONTINUE 

CBN T I lUE” ' " 

I SK I P- I SCAN+- 1 
WRITE' ( L T N 6 , 110) 

W p I TE ( L T N 6 , 10 5) 

R F A n” { L TN 5 , 106 ) ” ~ ' - 

WRITE (LTN6, 105) ' 

WRITE (LTN 6, 105) ' ' ' 

WRITE ( LT.N6, 105 ) 

“WRITE “ (LTN6, 1 0 3') 

D ' 50 I A P= I SKIP, NSC AN 

I )' 1 J=1,12 " ’ - -- -- - 

XMPi'K J )=0.0 
h ’ i i = i , i 2 
cornij, n = o.o 
c/xrT.up 

READ (LTN 5, 100) MPEG I N , NEND 
WRITE (LTN6,104) ' ‘IABjNBEGIM.n'eKD ‘ 

WRITE ( LT M6, 10 5) 

NS AMP = NEND- WEGIN+f ” “ ** 

SAMP* IS AMP I 

N S T A R T = \ 8 E GPL- 1 ’ ' 


IF(NSTART) 62,61,62 
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62 

CONTINUE 

DO 10 I 1 = 1 » NSTART 


* 



10 

READ (LTN) 

IOUM 




61 

CONTINUE 





C 

INPUT DAT READ 

FROM 

LTN 




00 20 I =NBEGI N,NEND 





READ (LTN) 

ICHEK, 

IDUM, (DATA! J) ,J=1, NCHAN) 




D0 26 L= 1, NC0MP 
DO 21 J=l, NCHAN 


' 




C0RR(L, J)=C0RR(L, 

J)+OATA(t)*DATA( J) 



‘ 21 

CONTINUE 





26 

CONTINUE 




* 

C 

MEAN VECTOR 

CALCULATED 




DO 20 K=l, NCHAN 






XMEANtK )=XMEAN(K)+DATA(K) 




20 CONTINUE 


NTERM=N£ND+1 

D0 25 I=NT£RM,NT0TAL 

25 READ (LTN) ICHEK 

SM=M$ 


MS=MS+ L _______________ 

DO 22 J = 1, NCHAN 

XMEANC J)”XMEAN( J)/(SAMP) 

22 CONTINUE ~ ' 

C C0RR ELA JJJN AND COVARIANCE MATRIX CALCULATED 

00 23 J = 1 f NC0MP 

oe 2 3 1 =1, NCHAN 

C0RR( J, I )=C0RR{ J, I) /{ SAMP) ~ 

23 CONTINUE 

00 51 J = I , NC0MP 

XBAR(J) =(SM/(SM+1.0) )*XBAR( J) + ( 1 . 0/ ( 5 M+1 .0) ) *XMEAN( J ) 

00 51 1=1, NCHAN 

AC0RRtJ,n = (SM/(SM+l.on*AC0Ra<J.II-Hl.C/(Sfr+l.O) )*C0RR(J,I ) 
51 CONTINUE 

50 CONTINUE 

WRITE ( LTN6 , 1 10 J 

WRITE ( LTN6 , 105 ) 

READ <LTN5,106) 

WRI TE ( LTN6, 1 06 ) 

WRITE <LTN6,105) 

WRITE { LTN6 , 107 ) 

WRITE ( L TN 6 » 1 0 5 ) ’ 

WRITE ILTN 6 ,H2) ( K ,K = 1 , NCHAN ) 

WRITE "{ LTN 6, 102) 1 0NE , ( XBAR ( J ) f J = 1 f NCHAN ) 

WRI TE ( LTN7 ,113) ( XBAR ( J ), J=1 , NCHAN ) 

WRITE (LTN6.105) 

WRITE ( LTN6 , 108) 

WRITE ( LTN6 , 1 05 ) 

WRITE ( LTN6 , 112) ( K ,K = 1 , NCHAN ) 

DO 53 J = 1 , NC0MP 

WRITE ( LTN6 ,102) J, (AC0RR(J,I), 1 = 1, NCHAN) 

53 CONTINUE 

00 54 J=1,NC0MP 

D0 54 1=1’, NCHAN 

_ ACQVR( J,I)=AC0RR(J,n--»(XBAR{J)*XBAR(I ) ) 
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54 C'JTI'JU'E 

W^ITE (LTN6,105) 

WkITE ( L T N J 6 1 1 09 ) 

,f_ITE (J. T M6tl05 ) 

WRTTG (LTN6.112) < K ,K = 1 , NCH AN } 

0 » 56 J = 1 , NC0MP _ 

W p l TE ( L T '16 » 102 ) J » ( AC? VR ( J » I j , I = 1 , i\iC HAM 
l),’ 2004 L- 1 * 1 2 

xxx xx ( j,l)=ac iv?(j, upscale 

2004_Cj' .TI'jiiE 

W*I TE ! LTM7 t 113) (XXX XX ( J , I ) , I = 1 , NCH AX ) 

56 C M IT1 -m E __ 

DO 60 J s 1 » NC0MP " ’ Z „ 

DO 60 I , NCM AN 

ABSC ■V-ABS(AC°iVR(J, J)=»AC0VR(I t I )) 

AC : >RR ( J » I ) = AC 0VR ( J , I ) / SORT ( A 8SC P V ) 

60 C IXT INUE 

WRITE (LTM6, 105) ” “ ~ ” 

WRITE ( LTM6 * l i 1 ) 

'WRITE (LTN6 , 1 05) 

WRITE ( L T‘16 , 1 12 ) JK»K«1«NCHAN] 

' D'' 1 59 J = l f NCP*MP "" 

WRITE (L TN6tl02) J , ( AC0RR ( J , I ) , I =1 ,KCHAN ) 

59 ~ CONTINUE 

I FLQ1=I FLD1+1 
IFI ELD= IFIELD+1 * 

1 F( I CHEK-2 ) 24 , 2_» 2 

24 W R I T E ( L T N 6 » 115) 

G * TO 3 5 

100 FORMAT ( 7110) 

101 FORMAT! ll,8X, 12F6^3) 

102 F OR HA f ( 2 X , 1 4, 4X , 1H* ,T 2F 10.5)' 

103 F3RMAK20X, 11HSCAN L INE , IPX , 5H S TART , 1C X , 4F STOP ) 

104 FORMAT! 2 OX , 1 5 » 1 4 X* I 5 , l OX , 15) 

105 FOR MAT! IF ) 

106" "FORMAT ( IX, 80H 

1 ) 

107 F ORMAT ( 3 OX , 14HMEAN VECTOR ) " 

108 F TPM A T ( 3QX , 20HC 0RRE LAT I 0N MATRIX) 

'109 'FORMAT (*3CXf 20HC0VARIANCE MATRIX') 

110 FORMAT ( 1F1 ) 

111 ' F CRM A T ( 20X , 32HN ORtMAL I ZED COVARIANCE tM A T R IX) 

112 F ’RMAT ( 4 X , 2HCH, 4X , 1 H* , 1 2 { I 5 1 5X ) > 

113 FORMAT! IX, i 2 F 6 • 3) 

114 F 'RMAT{ IX, 12F6,4)__ _ 

115 F 1RMAT ( l Y , 4OHN0 END 0 F~ FTtLO'"DES I GN AT I ON 0N THIS TAPE)' 

3 3 C£_N TINUE v 

■~”KS=.<S+1 

WRITE ! LTN6, 1001 ,KS _ 

REWIND LTN 

IF! I FLD1-MFLD1) 34_,34 JL 35 __ 

35" CO NT I WUF 

STO P 

END 
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3.7 PIECEWISE CORRELATION AND COVARIANCE PROGRAM 

This program is devised for calculating the piecewise mean, root-mean- 
square, and cross-correlation of a specified channel with respect to the rest 
of the channels in the multi-channel data set. The program also calculates 
the piecewise covariance, the second derivative of the piecewise correlation. 

The associated statistical errors of the above correlation or covariance are 
also calculated. 

DEFINITION OF INPUT PARAMETERS : 

CARD 1 

IPIECE - Option for piecewise printout ^ 0 for printout. 

LTERM - TERMINATES JOB 

4 - 0 for next case 
= 0 to terminate 

NTOTAL - Total NR. of samples in scan line. 

NARRNG - Rearrange input data for correlations. 

# 0 

CARD 2 

st 

NBEGIN - 1 data point to use in each scan line. 

NEND - Last data point to use in each scan line. 

ICHAN - Reference channel for correlations. 

NSAMP - No. of samples /piece, for piecewise calculations cannot be greater 
than (NEND-NBEGIN) . 

NORM - ^ 0 for normalization. 

JFIELD - No. of fields to process. 

The input deck setup is shown in Figure 3-11. 

The flowchart of the program is shown in Figure 3-12 which is followed by 
a listing of the program. 

3.8 PROBABILITY HISTROGRAM PROGRAM 

This program is devised for calculating the univariate (first order) 
probability density functions of tolie multispectral data set within a specified 
inventory area. The program can handle up to 12 channels of data simultaneously. 
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NBEGIN, NEND, ICHAN, NSAMP, NORM, JFIELD 
FORMAT- 14 


46 


1 


11111111112 
1234567890123456789012 


IPJECE , I AC CUM, LTERM, NTOTAL, NARRNG 
(FORMAT - 14) 
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Figure 3-12. FLOWCHART FOR CORRELATION AND COVARIANCE PROGRAM 
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C*************************************** ******************** 

C CORRELATION AND COVARIANCE PROGRAM W/ERRORS 
C*********************************************************** 

DIMENSION XMEAN (12 ) .COVR( 12 ) .ACOVR ( 1 2 > * BCOVR ( 1 2 ) .ECOVR( 12 ) , 
1FCOVR (12) . FRROR ( 12 ) 

DIMENSION F CORR (12) 

DIMENSION ARMSK13) 

DIMENSION DATA( 150.12) .TEMP ( 150. 12 ) .CORR ( 12 ) .DAT ( 150 . 12 ) 
DIMENSION ACORR (12) .ECORR ( 12 ) *DUM ( 12 ) 

DIMENSION DATS ( 150.12) .ERROR (12) * X DAT C 1 2 ) 

DIMENSION ARMS ( 13) .BRMS(13 ) »BCORR( 12) .CCORRt 13) 

EQUIVALENCE ( TEMP.DAT ) 

LTN=54 
LTN 1*6 1 
L T N 2 = 6 0 
2000 CONTINUE 

READ (LTN2.100) I P I ECE » IACCUM .LTERM.NTOT AL .NARRNG 
I CHEK=0 
I F I ELD=0 
JF I ELD = 1 
ICHEK=0 

200 CONTINUE 
IFIELD=IFIFLD+1 

C PUNCH 116. IFIELD 

C 

C JFIELD IS NR. OF CROPS + 1 

C TO TERMINATE JOB. 

C 

I F ( IFI ELD-JFIELD) 201.201.113 

201 CONTINUE 
I CHEK1=0 

DO 302 1=1.12 
ACOVR { I ) =0.0 
ECOVR ( I ) =0.0 
BCORR ( I ) =0 .0 
BCOVR ( I ) =0. 0 
FCORR ( I ) =0 * 0 
FCOVR ( I ) =0. 0 
ACORR ( I ) =0.0 
ECORR ( I ) =0.0 
302 CONTINUE 
MS5 = 0 
MS = 0 

NCOUNT =0 

WRITE (LTN1.108) 

2 CONTINUE 
C 

C READ PARAMETERS FOR EACH SCAN LINE 

C 

READ (LTN2.100) NBEGI N »NEND . I CHAN • NSAMP .NORM » JF I ELD 
NTERM=NEND+1 
NEND=NFND~NBEGI N+l 
MSS=MSS+1 
SAMP=NSAMP 
JEND=NBEGIN-1 
IF (JEND) 501*501.500 
500 CONTINUE 
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DO 10 11=1# JEND 
10 READ! LTN * 100 ) IDUM 
501 CONTINUE 

DO 20 I*1*NEND 

READ ( LTN *101) I CHEK ♦ < DAT < I . J ) . J 
20 CONTINUE 

DO 23 I=NTERM*NTOTAL 
23 READ( LTN *101 ) ICHEK 
PRINT 100.ICHEK 
IFCNARRNG) 50*35*50 
50 CONTINUE 
KM*MSS/2 
XM-MSS/2 
BM*KM 
XM=XM-BM 
KM»2.0*XM 
IF(KM- 1)35*32*32 

32 JJ=NEND 
KKK = 1 

DO 33 1 1 =1 *NEND 

DO 28 J*1 * 12 

DATS! KKK * J ) =DAT ( J J » J ) 

28 CONTINUE 
KKX=KKK+1 
JJaJJ-1 

33 CONTINUE 

DO 305 1*1* NEND 
DO 305 J*1 * 12 
DAT ( I * J ) “DATS { I * J ) 

305 CONTINUE 

35 CONTINUE 

DO 36 11*1 * NEND 

NC0UNT=NC0UNT+1 

DO 37 J=1 * 12 

DATA ( N COUNT * J ) =DAT ( I I * J ) 

37 CONTINUE 

36 CONTINUE 
KK= 1 

DO 25 I * 1 *NCOUNT 
LL* 1+1 

DO 24 J* 1 * 12 

TEMP(KK*J)=DATA( I » J J-DATA < LL * J ) 

24 CONTINUE 
KK=KK+1 

25 CONTINUE 

IF(NCOUNT-NSAMP) 2002*2003.2003 

2002 CONTINUE 

PRINT 100 * ICHEK 
IFIICHEK-2) 2004*200*200 
2004 CONTINUE 

IFINCOUNT-NSAMP) 2*2003.2003 

2003 CONTINUE 
MCOUNT =0 
I COUNT =0 

26 CONTINUE 
MS*MS+1 

DO 1 J*1 *12 


1 * 12 ) 
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COVR ( J ) -O# 0 
ARMS1 ( J) *0.0 
XME AN ( J) a 0.0 -'t 

ARMS ( J ) =0*0 
CORR ( J )=0*0 
ERROR! J) *0.0 
ERROR { J ) *0.0 
1 CONTINUE 

DO 403 J*1 » NSAMP 
MCOUNT *MC0UNT+1 
DO 403 1*1*12 

XMEAN ( I ) “XMEAN ( I ) +DATA ( MCOUNT » I) /SAMP 
403 CONTINUE 

DO 30 KK=1 *NSAMP 
I COUNT = 1 COUNT+1 
DO 29 J* 1 * 12 

COVR ( J ) *COVR ( J ) + < DATA ( I COUNT * ICHAN > -XMEAN t I CHAN ) ) 

1 *( DATA! I COUNT, J) -XMEAN ( J) ) 

CORR ( J ) *CORR ( J ) +TEMP! I COUNT . ICHAN) *TEMP { ICOUNT » J) 

ARMS ( J ) “ARMS ( J ) +TEMP ! I COUNT * J ) #TEMP( ICOUNT.J) 

ARMSK J)=ARMS1 ! J)+DATA< ICOUNT,J)*DATA( ICOUNT* J) 

29 CONTINUE 

30 CONTINUE 
SM*MS-1 

DO 207 J = 1 * 12 

COVR ( J ) =C0VR ( J ) /SAMP 

ARMS1 ( J ) =SQRT ( ( ARMS1 ( J ) /SAMP J-XMEAN ( J ) **2 ) 

ARMS ( J ) =SQRT { ARMS ( J ) /SAMP ) 

C0RR{ J)«CORR< J) /SAMP 

ACORR ( J ) = ( SM/ ( SM+1 .0 J ) *ACORR ( J ) + < 1 .0/ < SM+1 .0 ) ) *CORR ( J ) 

ACOVR ( J ) = ( SM/ ( SM+1 * J )*ACOVR( J) + ( 1.0/ (SM+1,0) )*COVR! J) 

ECORR ( J) =( SM/ ( SM+1 .0 ) ) *ECORR ( J ) + ( 1 . 0/ < SM+1 .0 )) * ( CORR ( J ) **2 ) 
ECOVR ( J) = ( SM/ ( SM+1. ) )*ECOVR ( J)+ ( 1.0/ (SM+1*0) )*COVR( J)**2 
I F ( ECORR ( J)-ACORR( J)**2) 207,207*206 

206 ERROR ( J ) =5QRT ( ( 1 . /SM ) * ( ECORR ( J ) -ACORR ( J > **2 > > 

207 CONTINUE 

I F ( ECOVR ( J ) -ACOVR { J ) **2 ) 204 *-204,20 5 

205 FRROR ( J J “SORT ( ( 1 • 0/SM )*( ECOVR (,J ) “ACOVR ( J )**2 > ) 

204 CONTINUE 

IF(IPIECE) 3001,3000,3001 

3000 IFUCHEK-2) 3003,3002,3003 
3002 CONTINUE 

3001 CONTINUE 

PIECEWISE PRINTOUT 


3003 


WRITEtLTNl 
WRITE(UTN1 
WRI TE( LTN1 
WRITE( LTN1 
WRITE( LTN1 
WRITE! LTN1 
WR I TE < LTNl 
PUNCH 115, 
PUNCH 115, 
WRITE(LTN1 
CONTINUE 


*105) 

*106) MS , ( 
, 600 ) MS , ( 
,110) MS , ( 
,603) MS , ( 
,112) MS , ( 
♦ 602) ,MS ♦ ( 
MS, (CORR < 
MS, (COVR ( 
*109) 


CORR 

( J) 

*J* 

COVR 

( J) 

, J“ 

ACORR 

( J) 

» J= 

ACOVR 

( J) 

,J* 

ERROR 

< J) 

,J* 

FRROR 

( J) 

» J= 

J) »J* 

1,12) 

J) *J= 

1,12) 


* 12 ) 

,12) 

, 12 ) 

, 12 ) 

* 12 ) 

* 12 ) 
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I F ( NORM ) 321*320.321 

321 CONTINUE 

CALCULATE NORMALIZED CORRELATIONS* COVARIANCE 
DO 323 J=1 * 12 

CORR(J)=CORR(J)/ (ARMS (J)*ARM5< I CHINA) > 

COVR( JJs COVR ( J ) / ( ARMS 1 ( J ) +ARMS1 < I CHINA) ) 

BCORR ( J)=(SM/(SM+1.0) ) *BCORR ( J ) + ( 1 . 0/ < SM+1 .0 ) ) *CORR ( J ) 
BC0VR(J) = (SM/( SM+1 • ) ) *BCOVR ( J )+ ( 1*0/ ( SM+1 *0) ) *C0VR { J ) 

FC0RR ( J ) = ( SM/ ( SM+ 1*0)) *FCORR ( J ) + ( 1 • 0/ ( SM+1 . 0 ) ) * ( CORR ( J ) **2 ) 
FCOVR ( J) =( SM/ ( SM+1 . ) ) *FCOVR { J ) + { 1 * 0/ ( SM+1 * 0 ) ) *COVR ( J ) **2 
I F ( FCORR ( J ) -BCORR ( J ) **2 ) 3 1 8 » 3 19 ♦ 3 1 9 

319 ERROR ( J ) =SQRT (( 1 • /SM )*( FCORR < J ) -BCORR ( J ) **2 ) ) 

GO TO 322 

318 ERROR ( J ) =0 • 0 

322 IF (FCOVR ( J ) -BCOVR ( J ) **2 ) 317 *316*316 
316 FRROR ( J ) =SQRT (( 1 . O/SM )*( FCOVR { J ) -BCOVR ( J ) **2 ) ) 

GO TO 323 
317 FRROR ( J ) =0 • 0 

323 CONTINUE 
IF(IPIECE) 3004*3005*3004 

3005 I F ( I CHEK— 2 ) 320*3004,320 
3004 CONTINUE 

WRITE (LTN1 *325 ) 

325 FORMAT (35X*47H********** NORMALIZED CORRELATIONS ******** 
WR I TE ( LTN1 * 105 ) 

WR I TE ( LTN1 , 114 ) MS • C ARMS ( J ) ,J=1*12) 

WRITE! LTN1, 106) MS, (CORR(J) *J*1*12) 

WR I TE ( LTN1 » 600 ) MS , ( COVR ( J ) »J=1,12> 

WRITE! LTN1.110) MS* (BCORRt J) *J=1*12 ) 

WRITECLTN1.601 )MS* ( BCOVR ( J) *J=1,12 ) 

WR 1 TE ( LTN1 *112) MS, { ERROR ( J) *0=1*12) 

WRITE(LTN1»602)MS*{ FRROR ( J ) * J=1 * 1 2 ) 

WR I TE ( LTN1 *605 ) MS » ( XMEAN ( J) * J=1 * 12 ) 

WR I TE ( LTN1 *604 ) MS » { ARMS1 ( J) »J = 1*12) 

PUNCH 115* MS* (CORR(J) ,J=1*12> 

PUNCH 115* MS* (COVR(J) ,J=1 ,12) 

WR I TE { LTN1 , 109 ) 

320 CONTINUE 
NCHEK=NCOUNT*I COUNT 

CHECK FOR SAMPLES LEFT IN SCAN LINE TO SORT AND GROUP 
WITH BEGINNING OF NEXT SCAN* 

NCHEK = 0 - NO SAMPLES LEFT ON PRESENT SCAN LINE 

I F(NCHEK-NSAMP) 44*26,26 

44 NCOUNT =0 
IF(NCHEK)45* 98*45 

45 CONTINUE 
NSTOP= ICOUNT+NCHEK 
DO 40 1 1 =1 COUNT *NST0P 
NCOUNT=NCOUNT+l 
DO 41 J=1 , 12 

DATA(NCOUNT*J)=DATA( II ,J) 

41 CONTINUE 
40 CONTINUE 


3-51 



TR-860 


NORTHROP 

HUNTSVILLE " 


98 IFUCHEK-2) 2*200*200 ! ; 

100 FORMAT (614) 1 

101 FORMAT ( I 1*6X*12F6.3) 

105 FORMAT <1X//7X*4HCH 1»5X,4HCH 2»5X*4HCH 3-*5X*4HCH 4,5X,4HCH 
15X * 4HCH 6 *5X *4HCH 7*5X»4HCH 8»5X*4HCH 9 * 5X * 4HCH10 *5X *4HCH1 1 * 
24HCH12/) 

106 FORMAT( 1X»I3*12F9.6*8H PC CORR/) 

110 FORMAT ( 1X*I3*12F9«6»8H AC CORR/) 

112 FORMAT {1X»I3»12F9«6*8H ER CORR/) 

114 FORMAT (1X»I3*12F9*6»8H RMS COR/) 

115 FORMAT { I 3 »2X * 1 2F6*4 ) 

116 FORMAT (20X*14H****** FIELD I3*8H ******) 

108 FORMAT { 1H1// ) 

109 FORMAT < IX/ ) 

600 F0RMAT(1X*I3*12F9*6*8H PC COVR/ ) 

601 FORMAT (1X»I3*12F9#6*8H AC COVR/) 

602 FORMAT 1 1 X, I 3 * 1 2F9.6 *8H ER COVR/) 

603 FORMAT (1X*I3*12F9#6*8H AC COVR/) 

604 FORMAT {1X*I3*12F9«6*8H PC STD /) 

605 FORMAT (1 X ♦ I 3 » 1 2F9 .6 *8H PC MEAN/) 

113 REWIND LTN 

IF(LTERM) 2000*2001*2000 
2001 CONTINUE 
STOP 
END 
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The main outputs are univariate probability histrograms for every channel. In 
addition, a normal distribution curve is superimposed on the histrograms. The 
normal curve has the mean equal to the most probable value of the calculated 
probability density function and the same standard deviation as the calculated 
one. 


It is optional to read a header identification record at the beginning of 
each data (IHEAD / 0) or to input some header identification via NAMELIST 
variable IDENT (IHEAD = 0). The one selected can be printed out as a title 
heading for each data set for the plotted probability density functions and 
for the printed mean and variance. Bit manipulation on a single variable on 
the read in header identification may also be done (IFLD f 0) . 

The program uses as input a tape reformated by the tape conversion program 
or a standard FORTRAN and/or Binary tape. The program's input control parameters 
are inputed by NAMELIST and the calling routines are NAMFIL and NAMDAT. 


DEFINITION OF INPUT PARAMETER ; 

NAMFIL 

• NFILE1 - Total number of data sets (files) to be processed cannot 

exceed four. 

NAMDAT 

• NBEAMS - Total number of channels of data to be processed (cannot 

exceed 12) 

9 NWHICH - Specify the numerical value of the files on tape to be processed 
(not to exceed 4) . Does not apply to tapes reformated by the 
tape conversion program. 

« TIMSAM - Specifies the time displacement between samples in the same 
channel. 

• LTN - Logical tape number of the input tape. 

0 NBLK - Total number of variables for each channel in a data record. 

For a tape reformated by the Tape Conversion Program, it specifies 
the number of samples in a scan line. 

0 NBEGIN - The sample number of the data value in a scan line to start 
with that will be used in probability density calculation. 

e NSTOP - The sample number of the data value in a scan line to stop 
with that will be used in probability density calculation. 
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• NTOTAL - The total number of intervals in the range of the values to be 

considered in the probability density function. 

• TINCRE - The increment of data value range considered m the probability 

density calculation. 

• ISKIP - Number of data records to skip at the start of each data set. 

• NSCANS - Total number of scan lines in a data set. 

• LTABLE - Minimum value of the data range to be considered in the 

probability calculation. 

• IBCD - Option to read header identification record at start of each 

data set in either Binary (IBCD = 0) or Alpha-Numerically 
(IBCD # 0). 

• IFLD - Option to do bit manipulation on a variable in the header 

identification record. IFLD + 0 bit manipulation executed. 

• NA - Starting bit position in NBLK where bits will be transferred in 

direction left to right. 

• NAB - Number of bit to be transferred. 

c NB - Starting bit position in RCDID(ISVB) where bit will be transferred 
in direction left to right, 

• ISVB - Array position in RCDID from which bit will be transferred to 

NBLK. 

• IDWORD - Number of header identification words to be read in. 

• IMODE - Mode that system REDTPC is to be read in. REDTPC is documented 

in 7094 Computer Manual of Marshall Space Flight Center. 

• IHEAD - Option to read header identification record. IHEAD ^ 0 header 

record read. IHEAD = 0 identification put on with IDENT 
NAMELIST variable. 

• IRDTPC - Option to either read input tape under system REDTPC (IRDTPC 

/ 0) or under format of output from tape conversion program 
(IRDTPC = 0). 

» IDENT - Labeling used in place of header identification record. 


Figure 3-13 shows the flowchart of the computer program, which is followed 
by a complete listing of the program. Some plot output^ are shown in Figures 
4-31 through 4-40. The symbol A's designate the actual calculated probability 
curves, while the symbol B’s designated an expected normal distribution curve 
with its mean equal to the most probable value of the actual curve and with the 
same standard derivation as the actual curve. Figure 3-14 is a printout result 
of the program. 
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EARTH - EFN SOURCE STATEMENT - IFN($) 


C MULTI -CHANNEL PROBAAB I L I TY DENSITY PROGRAM 

DIMENSION NHEG I M( 27 ) , NST0P ( 27) 

DIMENSION SAM P ( 12) 

DIMENSION R TABLE! 12), STABLE < 12) 

DIMENSION ANAX( I 2 ) , P ( 2 0 1 ) 

DIMENSION <» CO ID ('15") , DATA! 5000) * D A T (~50 00 ) ♦ f A R L E (201 ) » PUM 20 1 » 12) 
DIME *SI3N AVGCL ( I?) , AAVGCL ( 12) , RMS ( 12 ) , ARMS { 1 2 ) , S IG ( 12) , AS I G ( 1 2 ) 
DIMENSION XARRAYl L2 > , XLAM( 12 ) ,N WHICH! 5 ) , I DEIST (5) 

1007 FORMAT ( 12A6 ) 

EQUIVALENCE ( RC DID ( IT ) , NC0UNT ) 

CALL CAMRA V ( 9 ) _ _ 

NAMELIST /NAMDAT/ N3E AM S » NWH ICH » Tl MSAM , L TN » IS ? LK j \BEG "l N, 

INSTOP, N TOTAL, TINCR£,I SK IP , NSCANS ,L TABLE , I BCD , 

2IFLD.NA, .NAB, NR , I SUB , I DWORD , I MODE ,~I HE AD ,‘l ROT PC , IDENT 
READ 15,1007) (XARRAY(I),I = 1,12) 

REAO( 5, 100 7) ( X L AM ( I ) , l = 1,12) 

NAMELIST /NAME I L/ N FI LEI 
READ ( 5, NAMF I L ) 

D 0 798 I AA = L , tF 1 L£ 1 
C INPUT NAMELIST DATA 

READ (5, NAMDAT) 

WRITE (6, NAMDAT) 

C INI TIL I ZE CONSTANTS 

TIME = 0.0 

NCT =NtfLK 

100A FORMAT { IH 14F8.4) 

DELT = TIMSAM 
DO 0 1=1, N TOTAL 

P( I ) = 0.0 _ __ _ 

tablk ( i ) = o.o ' 

T CONTINUE 

DO 106 I = 1, NBEAMS 
GTARL E ( I) =0.0 
BTABL E { I) = 0.0 
S I G ( I ) = 0.0 
AMAX ( I ) = 0.0 

AAVGCL ( I ) = 0.0 
ARMS I I ) = 0.0 
106 CONTINUE 

DO 71 II = 1,12 
DO 71 I = 1,201 
71 PLAN ( I, II) = '0-0 

C SET DATA RANGE TO BE C3NSIDIERED IN PROBABILITY DENSITY 
TABLE ( I )=LTABLE 
D0 1011 I = 2 , NTOTAL 
TASLE(I) = TABLE CI-D+TINCRE 
1011 CONTINUE 

C POSITION TAPE FILES AND NUMBER 0F FILES 

I F ( IRDTPC.EO.O) GO TO l l 
REWIND LT:N 
IAB = NWHICHt I A A )-l 
IF ( I A 3 ) 11,11,12 


heproducibie 
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12 CONTINUE 

DO 10 IAC=1,TAB 
CALL 5K.FBIN19, l,RD) 

EARTH - £FN SOURCE STATEMENT - IFM(S) - 

10 csjnti iur- 

11 CONTINUE 

C READ HEADER I DENT I F I C ATI 0N RECORD 

IF< IHEAO.EQ.O) GO TO 1060 
I F ( IHCD.EQ.O) GO TO 1013 

REAO ( LTN, 105 1 > ( R CD I D ( I ) t I =1 , 1 DW3P0 ) 

1051 FORMAT ( 15A6 ) 

RCDIDf I SUB ) =0.0 
RCD I D ( I SUB )=NGLK 

•mct=-nblk 

GO TW 1 0 LA 
10 L3 CONTINUE 

READ (LTM) ( PC D I D { f ) i I =1 , 1 DWORD ) 

I F ( IFLD.EQ.O) GO T 3 1014 
NCT =0 

CALL FLD(NCT,NA,NA8,NB,RCDID{ l SUB) ) 

1014 CONTINUE 

NC0UNT = NCT 
NHLK = NCT 

WRITE (6,1000) ( RCD ID ( I ) , I = 1,15) 

1000 FORMAT ( LX, 12A6, I6,5X,?F8.4) 

1060 CONTINUE 

I F ( I He AD «N E .0 ) GO TO 1070 
WRI TE(6, 1071) 

WRITE16, 105L) ( I DE NT ( I ) , I *1 , 1 0 WORD ) 

WRITE(6,107l) 

1071 FORMAT (1H0) 

DO 1072 I = l » I DW 8 P D 
RCD ID { I ) = I DENT ( I ) 

1072 CONTINUE 
1070 CONTINUE 

WRITE (6, 24) 

24 FORMAT (17V, 3HCH 1 , 6X , 3HCH2 , 6X , 3HCH3 , 6X , 3FCH4 ,6X,3HCH5,6X , 

1 3HCH6 , t>X » LHCH 7 , 6X, 3HC H8 , 6X , 3HCH9 , 6X , 4HCH 1 0 , 5 X ,4HCH 1 L , 5X ,4HCH12) 
NCOUNT = NCT 
NIJP = JCOU'\iT*NBEAMS 
MS = 0 

C OPTION T 1 SKIP SCAN LINES 

UR I N= I SK I P + 1 

I F ( ISKIP.EQ.O) GO TS. 1032 
04 10 T 1 J = 1,ISKIP 
CALL SKRBI N ( 9 , 1 , RD ) 

1031 CONTINUE 

1032 CONTINUE 

DM I I AU=N R IN, N SCANS 

NWDS= ( NSTOP ( I AD ) + L-NBEG IN ( I AO) ) * { NREAMS + 1 ) 

DO 104 I = 1 » GBE AMS 
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AVoCL ( I ) = 0.0 
RMS ( I ) = 0.0 
S AMP { I ) =0.0 
104 CJNTMUE 
MS = MS+i 
IA = 1 

IF< IR0TPC.GE.0) G0 T 2 1061 
D? 60 J= L * NU P , U3EAM S 
JZ=J+NREAMS-l 

EARTH - EFN SOURCE STATEMENT - IFN(S) - 

READ { L T N ) I CHEK , ( DAT (L ) , L =J , J l J 
60 C3MTIwUE 
GO TO 1062 

1061 C 5UTI SUE 

CALL REDTPCILTN, I M0DE , I ERR » NW ,NUP »DAT ) 

GO T M ( 50,51,50,50* Si J, IERR 
51 GO T/ 999 
50 CONTINUE 

1062 CONTINUE 

NSTART = N3EGI.'!( I AD) 

MTERM = NST0P ( I AO ) 

D? 30 ICC=UST ART ,NTESM 
NFRM= ICC*NB5AMS-NBEAMS+ 1 
NT0=NFRM+MBEAMS-1 
J = 0 

D ') 31 I = MFRM, NTO 
J = J + l 

I F ( DAT ( I ).LT. TABLE! 1) ) G0 T3 70 
I F ( DATC I ).GT. TABLE! MTk'TAL) ) GO TO 70 
AVGCL ( J ) = AVGCL ( J ) + D A* ( I ) 

RMS(J) = RMS ( J ) +DAT ( I )**? 

70 l F( DAT ( I ).LT. TARLEt i) ) SAMP (J ) =S AMP ( J )+ L . C 

IF( OATtD.GT. TABLE ( 'ITCTAD ) SAMP ( J ) = SAMP { J ) + l . 0 
31 CONTINUE 

DAT A ( l A ) =T I ME 

T IME=TIME+DELT 

DO 'if I * NFRM,NT3 

IA = Ia+l 

DAT A ( I A ) =DA T ( I ) 

6 3 CONTINUE 
IA = I A + 1 

I F ( I A.GE.NWUS ) fjH Tv) 100 
30 CONTINUE 
100 CU'iTI'lUE 

SAMPU3=riST3P( TAD)-NBEGIN{ IADJ + l 
SM = MS 

DO 101 I = l , NBEAMS 

RMS ( I )=RMS( I ) / ( SAMPNO-S AMP ( I ) ) 

AVGCL ( I ) =AVGC L { I ) / ( SAKPN0-SAMP ( I ) ) 

ARMS ( I ) = ( ( ( SM-1. ) /SM)*ARMS(I ) +RMS(I)/SK) 

A AV jC L ( l ) = ( { (SM-l.)/SM)*AAVGCL(I)+AVGCL(I )/SM) 

SIt»( 1 ) = SORT ( RMS ( I ) -AVGCL ( I )**2) 

RMS ( I ) = SORT (RMS { I ) ) 
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WRITE (6,9000) ( AVGCLd) ,I=i,NBEAMS) 

WRITE (6,9001) ( RMS ( I ) ,1=1,NBEAMS) 

WRITE (6,9002) < SIGH) ,I=1,NBEAMS) 

9000 format ( ix, 8 hpc mean, ix , 12F9.4) _ 

9001 FORMAT ( LX, 8HPC RMS" , IX,12F9.4) 

9002 FORMAT ( IX, 8H D C SO , 1X,12F9.4) 

DO 94 I * 1 , N (TEAMS 
S I G ( I ) = SORT ( AR M S ( I ) — A AVOC L { I )**2) 

RMS ( I ) = SQR T (ARMS ( I ) ) 

94 C.MTI'IUE 

WRITE (6,9003) ( A AVGCL ( I) , I = 1 , N«3E AMS ) 

WRITE (6,9004) ( RMS ( I ) ,I=1,NBEAMS) 

EARTH - EFM S0URCE STATEMENT - IFN(S) 

WRITE (6,9005) ( SIG(I) ,I=1,NREAMS) 

9003 F3RMAT( IX, 8HAC MEAN, IX , 1 2F 9 .4 ) 

9004 FORMAT l IX, 8HAC RM S, 1 X , 1 2 F9 . 4 ) 

9005 F0RM AT { LX, 8HAC SO , 1X,12F9.4) 

II = 0 

MXYZ=MST2P( I AD ) -rMBEGI H( I AD ) + 1 __ 

09 3? T - l, NXYZ 
II = I I +1 

DO 32 III = 1 , N B E AM S 
L = 1 
II » Il+l 

33 CONTINUE 

I F ( DATAC II ) .LT.TACLEt 1) )G3 T0 30 
IF(OATA( II ) .GT.TARLECNT0TAL) ) G0 T3 30 
IF(OATA( II ) .GT.TABLE(L) ) GO T9 34 

G3 r? 30 

34 I F ( DA F A ( II ) .LT.TABLE( L+L) ) G0 TO 31 
G3 T 1 30 

31 PLAM ( L , III ) =' PLAM( L, II I ) + 1.0 
GO T« 32 

30 CONTINUE 

IF(DAT\( II ) ,LT.TABLc( l) ) 8TABLEUII) = BTABLE ( II I ) + 1 . 

I F ( DATA (II) .GT. TABLE I NT0TAL ) ) GTABLE ( I 1 1 )=GTABL£( 1 1 I >*1. 
IF(DATA{ II ) ,LT.TADLE( 1) )G0 10 32 
IF(DATA( II).GT.tABLE(NT0TAL) ) G 8 TO ‘32 
L =L + 1 

I F( L. LT . NT0TAL ) G0 TV 33 
GO TO 31 

32 CONTINUE 

100? F TRMAT ( 1H 1 3F 9 « 4 ) 

1003 FORMAT! 1H 12F9.4) 

1 CONTINUE 
999 CONTINUE 

O'* 90 1=1, NBE AMS 

1U 90 J = untotal 

IF{ PLAM( J, I ) .GT.AMAX( I ) ) G0 TO 92 

GO TO 90 

9? 4MA X ( I ) = PLAM( J , I ) 
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90 C '3 JTI «Ut 

WRITE: (6, 9006) 

9006 F2RMATI 35X, 37HNUMBER 3F VALUES LESS THAN TABLE) 
WRITE (6,1003) ( STABLE ( I ) » I * 1 » N B E A M S ) 

WRITE {6, 900 7) 

9007 FORMAT( 35X, AOHNUMBER OF VALUES GREATER THAN TABLE) 
WRITE (6,1003) { GTABL E l I ) , I = l,N n EAMS) 

WRI TE { 6 , 900 o ) 

9008 F ff R.»|AT ( 35X » 3 1 HM AX I MUM VALUE OF EACH TABLE) 

WRITE (6,1003) (_AMAX(I),I = UNSEAMS) 

NP * iT2TAL 

DO 9 3 I = UNBEAMS 

DB 95 J = 1 , NTSJTAL 

BASE=( TABLE (J ) - AAVGCL ( I ) ) **2/ l 2 .*SIGl I )**2) 

I F ( BASE. GE. 89 . ) Gv5 17, 73 
P(J) = A MAX ( I )/EXP(BASE) 

GO 17 72 
73 P(J)=0.0 


EARTH - tFN SOURCE STATEMENT - IFN(S) - 


72 


95 


93 

998 


CONTINUE 

CONTINUE 

CALL QIJIK3V(-l, LHA, XARR AY , XL AM 
CALL 0UIK3V ( 0 , 1HB , X ARRAY , XLAM 
CALL PRINTV(90,RCDH), 212, 1010) 
CONTINUE 


CONTINUE 


, -NP, TABLE, P LAM (1,1 ) ) 
,-NP, TABLE, P) 


CALL CLEAN 

ST?P 

END 
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CH3 
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CH7 CHB 
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.PC MEAN 
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0.9519 
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0.0506 
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AC MEAN 
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0.9055 
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0.0313 
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0.0518 

PC MEAN 
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0.6051 


PC 

RMS 

1.2618 
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0.8923 
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PC 
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PC 
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0.8966 
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PC 

SD 
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0.0513 
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AC 

mean 
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AC 
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AC 

SD 
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PC 
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PC 
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0.8356 

0.6181 
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0.8128 

0.5758 
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PC 

SD 

0.0790 
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3.9 DIVERGENCE MATRIX PROGRAM 

This program is devised for calculating the "divergence" measure which is 
a similarity criterion between two homogeneous populations whose probability 
density functions are normally distributed. The mathematical formulation of 
this divergence criterion is given in subsection 2.6. The program is coded 
such that up to 12 homogeneous populations can be compared simultaneously. The 
main input parameters are mean vectors, covariance matrices, and selected 
channels for all populations. The main output is a square matrix of divergence 
between all combinations of the given populations. The divergence matrix 
consists of two terms, which are called Trace 1 and Trace 2 matrix representing 
the first and second terms of the right hand side of equation (2-11) in Section 
II. 

Four subroutines are used in this program: Matrix Multiplication (MMPY) , 
Matrix Subtraction (MSUB), Matrix Addition (MADD) , and Matrix Inversion (MINV). 

DEFINITION OF INPUT PARAMETERS : 

CARD 1 (7110) Format 

a NSCLE1 - Scale factor for Mean Vectors, 

a NSCLE2 - Scale factor for Covariance Matrices. 

e NCYCLE - Number of times to cycle through the program comparing the 
matrices. 

a KKK - Maximum row size of Covariance Matrix to be read. 

o III - Maximum column size of Covariance Matrix to be read. 

o NLOOP - Number of times to cycle through the matrices for different 
row and column configurations . 

« NMEANS - N um ber of Mean vectors or Covariance Matrices to be used in 
the calculation. 

CARD 2 (3110) Format 

9 NLOW - Lower range of field identification to be read, 
o NUP - Upper range of field identification to be read. 

e NPRINT - Option to print out intermediate calculation of cross products 
of comparisons of two populations. 

CARD 3 (20A4) Format 

o Z(L) - Array in which field identifications are stored. 
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CARD 4 (1215) Format 

• ICH(L) - Specifies which columns of Covariance Matrices to use in 

calculations. ICH(L) - 0 column not used. 

CARD 5 (1215) Format 

• NSC(I) - Specifies which rows of Covariance Matrices to use in 

calculations. NSC(l) - 0 row not used. 

CARD 6 (5110) Format 

e NSCTS - Specifies number of rows actually used in Divergence 
calculation from NSC (I) selection. 

• ICHAN - Specifies number of columns actually used in Divergence 

calculation from ICH(L) selection. 

• ICHAN1 - Used to specify that the mean is a 1 X NCHAN vector. Nominal 

value is one. 

e NCHAN - Maximum row size of covariance Matrix to be read. 

• NSCANS - Maximum column size of Covariance Matrix to be read. 

CARD 7 (12F6.Q) Format 

e YBAR(KK,X) - Array that Mean vectors for each population is stored. KK 
specifies population and I the dimension of the vector. 

CARD 8 (12F6.0) Format 

e X(KK,J,I) - Array that Covariance Vectors for each population is stored. 

KK specifies population, I the row dimension, and J the 
column dimension. 

Figure 3-15 shows a flowchart of the divergence matrix program, which is 
followed by a complete listing of the program. Some sample outputs are shown 
in Tables 4-1 through 4-5. 


3.10 Si-TABLE PROGRAM 

This program is devised for calculating the so-called Si-value which 

represents a measure of similarity between two inventory areas which, in turn, 

are characterized by the mean vectors and the standard deviations. Let the 
til til 

mean vectors for i and j inventory areas be denoted by M.(k) and M. (k) , 

J 

and the standard deviations be o^(k) and a . (k) , where k stands for the channel 

1 tli til 

number, (k=l,2...K). The Si-value between i and j inventory areas can be 
expressed as 
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Figure 3-15. FLOWCHART OF DIVERGENCE MATRIX PROGRAM 
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C PROGRAM TO CALCULATE DIVERGENCE OF TWO MATRICES 
DIMENSION Z ( 240 ) 

DIMENSION Z8ARC 12*12) 

DIMENSION X1II2, 12 ) 

DIMENSION DIVfl2,12T 

DIMENSION V! 12,12) ,V1<12, 12) ,VT( 12,12) ,VT 1(12, 12) 

DIMENSION VPROD ( 12, 12) , VXPROD! 12, 12 ) , VDIFF ( 12 , 1 2 ) , VTD IFF { 12, 12) , 
AV SUM { 12, 12) ,VTSUM(12,12) 

DIMENSION XBAR( 12, l 2 ) , XTBAR ( 1 2 , 12 ) , XRRODI 1 2 , 12 ) 

DIMENSION X( 12, 12, 12) , I CHI 12) , NSC (12) 

DIMENSION YB-ARI 12, r21 " ~ - • 

DIMENSION TRACE 112,12) , TRACE II 12,12) 

1000 FORMAT (7110) 

1001 FORMAT ( IX, 12F6.0) 

1002 FORMAT ( 1H0) 

1003 FORMAT ( 80H 

1 “ ) 

1004 FORMAT (IX, 5HDETV=F13.8, 6HDETVI=F13.8,6HTRACE=F10.3, 7HTRACE1=F10. 3 ) 

1005 FORMAT (2X, 1 4, 4X , 1H* , 12F 10 . 3) 

1006 FORMAT ( 1 X , SOH#*#***^#* 1 *#*************##* 1 !' *$$$$*****4##**#*****#*^ 
1 $*$$«***********#*#**##***)!<* ) 

1007 FORMAT ( 4X, 2HCH , 4X, 1H* , 12 ( I5,5X) ) 

1008 FORMAT ( 1215) 

1009 FORMAT I 30 X»4HDIV=F 10.3) 

1010 FORMAT { 1H 1 ) 

1011 FORMAT ! 2X , I4,4X, 1H*, 12F10.5) 

1012 F0RMAT(3OX,20HC0VARI ANCE MATRIX) 

1013 FORMAT ( 30X , 14HMEAN VECTOR) 

1014 FORMAT 1 30X , 29H I INVERSE “ COVARIANCE MATRIX) 

1015 FORMAT! 30X,39H(Vi( J,I )-V2( J, I ) ) ( I NV2 ( J , I ) - I NV1 ( J, I ) ) ) 

1016 FORMAT ( 1 OX , 79H( I'NVl ( J , I ) - I NV2 ( J, I ) ) ( X6ARI ( J , I) -XBAR2 ( J ♦ I) ) (TRANS 
lP0SE(XbARl(J,I)-XBAR2(J,I) )) 

1017 FORMAT I30X,8HFIELD ,15) 

1018 FORMAT ( 30X , 34HPR0DUCT OF MATRIX WITH INVERSE) 

1019 FORMAT ! 30X , 20HD T VERGENCE MATRIX) 

1020 FORMAT { 30X ,16HTRACE1 MATRIX) 

1021 FORMAT! 30X, 16HTRACE2 MATRIX) 

1022 FORMAT! 30 X, 36HCHANNELS USED IN THIS CALCULATION) 

1023 FORMAT ( 20A4 ) 

1024 FORMAT! IX, 5HF I EID,4X, 1H*, 12! 15, 5X) ) 

1025 FORMAT! 10X,5HFIELD, I4,4* f 8HC0MPAHED,4X,5HF I ELD, 14) 

IGNE=1 

READ 1000,NSCLE1,NSCLE2,NCYCLE,KKK, III , NL OOP, NME ANS 

NCOVR=NMEANS 

SCALE 1=NSCLE1 

SCALE2=NSCLE2 

DO 251 IFF=1»NL00P 

READ 1000,NLOW,NUP,NPRI NT 

PRINT 1010 

PRINT 1006 

NBEGIN=NLOW 

NSTQP=NUP 

DO 205 J J= 1 , NMEANS 

IF! IFF-1 ) 2012,2013,2012 
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2013 READ 1023, < Z I L ) , L=NB£GIN ,NSTOP ) 

2012 PRINT 1023,(Z(L),L=N6EGIN,NSTOP) 

NBEGIN=NBEGIN+NUP 

N$T0P=NUP+NBEGIN-1 

205 CONTINUE 

PRINT 1006 “ “ 

DO 206 J = 1 , KKK 
DO 206 1=1, III 
DIVI J, I ) =0.0 
TRACE I J , I)=0.0 
TRACE1 I J , I ) =0.0 
ZBAR (J,I)=0.0 ' 

V( J,I)=0.0 

206 CONTINUE 

READ 1008, ( ICHI I ) , 1=1, KKK) 

READ 1008, (NSC(I), 1=1, III) 

READ 1000 , NSCTS, I CHAN, ICHAN1,NFILE,NCHAN, NS CANS , NSC ALE 

SCALE=N$CALE 

MS= 1 

DO 200 I AB= 1 , NCYCLE 
IFIMS-1) 203,202,203 

202 CONTINUE 

IFIIFF-l) 203,2003,203 

2003 CONTINUE " ' 

DO 21 KK=1 , NCOVR 
PRINT 1010 
PRINT 1002 
READ 1003 
PRINT 1003 

PRINT 1002 ' ' 

PRINT 1013 
PRINT 1002 

READ 1001, CY8ARIKK, I ) , I =1 »NCHAN ) 

DO 201 1=1, NCHAN 
YBARIKK, I )=Y8AR(KK, I) /SCALE 1 
201 CONTINUE 

PRINT 1007, (K,K=1, NCHAN) 

PRINT 1005, ( KK , { YBAR ( KK , I) , I =1 , NCHAN ) ) 

PRINT 1002 
PRINT 1012 
PRINT 1002 

PRINT 1007, (K,K=l,NCHAN) 

DO 21 J= 1 » NSCANS 

READ 1001, XX(KK,J, I ), 1=1, NCHAN) 

DO 5 1=1, NCHAN 
X(KK,J, I)=X(KK,J, I) /SCALE 
5 CONTINUE 

PR INT TO IT, TO,(X(KK,J,I) ,1 = 1, NCHAN ) ) 

21 CONTINUE 

203 CONTINUE 

DO 402 KK= 1, NMEANS 
M = 0 

DO 40 1=1, NCHAN 

IFIICHim 4'2»40*~42' ' 

42 M=H+i 
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ZBAR(KK,M)=YBAR(KK, I ) 

40 CONTINUE 
402 CONTINUE 
M=0 

DO 22 J=1»NCHAN 
IFdCFUJ)) 24,22,24 

24 M=M+1 
N = 0 

DO 23 1=1, NSCANS 
IFINSCC l ) ) 25,23,25 

25 N=N+1 

vnr, m > =x cits, i tj) 

23 CONTINUE 
22 CONTINUE 

C MATRIX INVERSION AND 

CALL MINVl V, VT,DET, ICHAN) 

DETl=DET 

CALTL MMPTt ICHAN,1NSCTS» ICHAN, V, VT ,Xl ) 
IF (NPRINT) 2007,2006,2007 
2007 CONTINUE 

IF (MS-1) 2004,2005,2004 

2005 CONTINUE 

PRINT 1010 
PRINT 1002 "" 

PRINT 1002 
PRINT 1002 
PRINT 1002 
PRINT 1006 
PRINT 1002 
PRINT r02Z 
PRINT 1002 

PRINT 1008, (ICH(I), 1=1, NCHAN) 

PRINT 1002 
PRINT 1006 
2004 CONTINUE 

PRINT lOtO ” 

PRINT 1002 
PRINT 1006 
PRINT 1017, IAB 
PRINT 1002 
PRINT 1006 
PRINT 1002- 
PRINT 1013 
PRINT 1002 

PRINT 1007, (K,K=1, ICHAN) 

PRINT 1005, t IONE, (ZBAR (IAB ,11,1=1, 
PRINT 1002 

PR TNT loir' - " 

PRINT 1002 

PRINT 1007, (K,K=1, ICHAN) 

DO 501 J=1 , NSCTS 

PRINT 1005, I J , ( V (J, I), 1=1, ICHAN) ) 
501 CONTINUE 

PR TNT TO 02" * 

PRINT 1014 


CHECK 


ICHAN) ) 
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PRINT 1002 

PRINT 1007, (K,K = 1, KHAN) 

DO 502 J=l, NSCTS 

PRINT 1005, { J, ( VH J, I) ,1=1, I CHAN) ) 

502 CONTINUE 
PRINT 1002 
PRINT 1018 

DO 503 J=l, NSCTS 

PRINT 1005, ( J, (XI t J, I) ,1=1, ICHAN) ) 

503 CONTINUE 
2006 CONTINUE 

NS=MS+ 1 ‘ " 

DO 20 I AA=NS* NF I LE 
DO 10 J = 1 , NSCTS 
DO 10 1=1, ICHAN 
XBARI J , I ) =0.0 

vkj, i ) = 0 . o 
VT1( J, I) =0.0“ 

VPROD I J , I ) =0. 0 
VXPRODI J, I ) =0 . 0 
VDIFFI J, I ) =0.0 
VTOIFFI J, I ) =0.0 
VSUM ( J , I ) =0.0 

vtsumi j,i)=7r.o ' ■ • 

XTBARI J, I ) =0. 0 
XPRQDI J, I ) =0. 0 
10 CONTINUE 
M=Q 

DO 30 J=1 ,NCHAN 

IF (ICHIJn 31730,31 ' • 

31 M=M+ 1 
N = 0 

DO 33 1=1 , NSCANS 
I F (NSC II)) 32,33,32 

32 N=N + 1 
Vl(N,M)=xrrAA,I,J) 

33 CONTINUE 
30 CONTINUE 

C MATRIX INVERSION AND CHEC< 

CALL MINV(V1,VT1,DET, ICHAN) 

DET2=DET 

CALL MMPVT ICHAN, NSCTS, ICHAN, VI, VT1, XI ) 

IF ( NPRINT ) 2009,2008,2009 
2009 CONTINUE 


PRINT 

1010 

PRINT 

1002 

PRINT 

1006 

PRINT 

1‘CJO 2'" ' * 

PRINT 

1017, IAA 

PRINT 

1002 

PRINT 

1006 

PRINT 

1002 

PRINT 

1013 

PRINT 

1007,1 K,K= 1,1 CHANT 

PRINT 

1005, ( IONE, (ZBAR (IAA , I ) , I =1 , I CHAN ) ) 
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PRINT 1002 
PRINT 1012 
PRINT 1002 

PRINT 1007, ( K »K=i, I CHAN l 
DO 504 J=l, NSCTS 

P R I NT' 1 005 VTO t TV 1 {'J , I ) , r=T,~ICHAN ) ) 

504 CONTINUE 
PRINT 1002 
PRINT 1014 
PRINT 1002 

PRINT 1007, IK,K=i, ICHAN) 

DO 505 J=l, NSCTS 

PRINT 1005, ( J,(VTi( J,I ) » 1=1, ICHAN) ) 

505 CONTINUE 
PRINT 1002 
PRINT 1018 
PRINT 1002 

' PR I N I" rOUTrl'KfK=TnCffSNV " ' 

DO 506 J=1 ,NSCT5 

PRINT 1005, I J, { Xl{ J,IJ, 1=1, ICHAN) ) 

506 CONTINUE 
2008 CONTINUE 

DO 3 J = i, ICHAN1 

DO 5 1 = 1 » TCHAD - ----- - - 

X8AR( J, I )=ZBAR( IAS, I ) -Z3AR ( I AA, I ) 

3 CONTINUE 

DO 8 J=1 , I CHAN1 
DO 8 1=1, ICHAN 
XTSARI l , J ) =XBAR( J , I ) 

8 ~ CONTINUE 

CALL MHPYI ICHAN, ICHAN1, l CHAN, XTBAR » XBAR , XPROD ) 
CALL MSUBI ICHAN, NSCTS, V* VI , VOI FF ) 

CALL MSUBC ICHAN , NSC TS, VT 1 , VT , VTD I FF J 
CALL MADDI ICHAN , NSC TS, VT ,VTi , VTSUM ) 

CALL MMPYI ICHAN, NSCTS, I CHAN, VDI FF, VTD IFF, VP ROD) 
CALL WMPYr rCTTA NTNS CTS , ICHAN f VTSUM, XPROD, VXPROD) 
IF (NPRINT) 2014,2015,2014 
2014 CONTINUE 

PRINT 1002 
PRINT 1015 
PRINT 1002 

PRINT 1002- “ - ” ~ - 

PRINT 1025, I AB# IAA 
PRINT 1002 

PRINT 1007, (K,K=1, ICHAN) 

DO 507 J=l, NSCTS 

PRINT 1005, ( J, (VPROD < J, I ) , 1 = 1, ICHAN) ) 

507 CONTINUE 

PRINT 1002 

PRINT 1016 
PRINT 1002 
PRINT 1002 
PRINT 1025, IAB, IAA 

PRINT "1002 “ ” 

PRINT 1007, (K,K=1, ICHAN) 


3-70 


o rs 


TR-860 


NORTHROP 

HUNTSVILLE ~ 

DO 508 J =1 , NSCTS 

PRINT 1005, t J , ( VXPRODl J, I ) , I =1 » I CHAN) ) 

508 CONTINUE 
2015 CONTINUE 

KK=Q I 

J3=0 " 

00 4 J = l, NSCTS 
KK-KK+ 1 
33=33+1 

TRACE ( IAB,IAA)=TRACE(1AB» IAA)+VPR0D(KK»3J> 
TRACE1IIAB,IAA)=TRACE1< IAB, I AA ) +VXPRQD IRK, JJ ) 

4 ' CONTINUE' " * 

DI V( IAB»IAA)=( TRACE ( I AB» I AA) +TRACE1 ( I A8» I A A ) J/2.0 
PRINT 1009, DIV( IAB* IAA) 

PRINT 1004,DET1,0ET2, TRACE! I AB, I AA ) , TRACE If I AB, IAA) 
20 CONTINUE 
MS=MS+ 1 
200 CONTINUE 

PRINT 1010 
PRINT 1002 
PRINT 1002 
NBEGIN=NLQW 
NSTOP=NUP 

DO 2010 33=T, NMEAMS ‘ ’ 

PRINT 1023, (Z(U , L=NBEGIN»N$TOP) 

NBEGIN=NBEGIN+NUP 
NST0P=NUP+N8EGIN-1 
2010 CONTINUE 

PRINT 1002 

PRINT "1002 ' " 

PRINT 1022 
PRINT 1002 

PRINT 1008, ( ICHU) , 1 = 1, NCHAN) 

PRINT 1002 
PRINT 1019 
PRINT 1002 
PRINT 1002 
PRINT 1006 
PRINT 1002 

PRINT 1024, (K,K=1,NC0VR) 

DO 2000 I AB= 1 , NCQVR 

PRINT 1005, (TAB, ('OTVII AB, IAA) , IAA=1,NC0VR ) ) 

2000 CONTINUE 

PRINT 1002 
PRINT 1006 
PRINT 1002 
PRINT 1002 
PRINT' T002 
PRINT 1010 
PRINT 1002 
PRINT 1020 
PRINT 1002 
PRINT 1002 

PRINT 1006" 

PRINT 1024, (K,K=1,NC0VR) 
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DO 2001 IAB=1,NC0VR 

PRINT 1005, ( IAB, (TRACE < I AB , I AA ) , I AA= 1 , NOOVR ) ) 
2001 CONTINUE 

PRINT 1002 

PRINT 1006 V' * 


PRINT' T0C2 

PRINT 1021 
PRINT 1002 
PRINT 1002 
PRINT 1006 ' 

PRINT 1002 

- - -~-p-RTNT~ 102 ATtKTK =T »NC0VR1 
DO 2002 I AB=1 , NCOVR 

PRINT 1005,( I AB » ( TRACE 1 ( I AB, J AA) , I AA=1 , NCOVR ) ) 
2002 CONTINUE 

PRINT 1002 
PRINT 1006 

251 CONTINUE " 

STOP 

ENO 


SIBFTC TWO 

SUBROUTINE MMPY ( J , K, L , P, Q, R ) 

DIMENSION P( 12,12 ),Q( 12,12 )tR( 12, 12) 
DO 100 JJ=1,J 

DO 100 LL=1»L 

100 — R rJTiTllT=OVO “ 

DO 200 LL=l»L 

DO 200 JJ=1,J 

DO 200 KK=1 , K 

R ( J J , LEI =R"{ JJ , LL ) *'P ( JJ»KK ) *Q( KK, LL ) 
200 CONTINUE 

RETURN ' 

END 


SIBFTC THREE 

SUBROUTINE MADD ( M ,N, A»B» C ) 

DIMENSION A(12,I2),B(12,12),C(12,12) 
DO 2 I = 1 , M 
DO 2 J=1,N 

Cl I,J)=A< I,T)+lTCTi _ " 

2 CONTINUE 
RETURN 
END 


SIBFTC FOUR 

SUBROUTINE MSUB (M,N»A,B,C) 

DIMENSION A( 12, 12), B{12, 12), 0(12,12) 
DO 2 1=1 , M 
DO 2 J=1,N 

C ( I , J ) =7H I , J ) -B'{ I , J ) 

2 CONTINUE 
RETURN 
END 
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$ 1 8 FTC FIVE 

SUBROUTINE MINVI A,AI ,DET,M) 

DIMENSION A< 12, 12), AI< 12,12) , IORD( 12) 

N=M 

DO IOO 1=1, N 

DTI' 100 X=T»W 

100 AI(I,J)=A(I,J) \ 

00 51 1=1, N 
51 IORD( I )=I 

DET=t. 

DO 1 K=1,N 

icpt =icfi" — 

TEST=ABS(AI(K,K) ) 

1 ND=K 

IF ( K-N ) 21 , 22, 21 

21 DO 23 L=KP 1 , N 
TEST1=ABS(AI ( L»K ) ) 

" TFTTE ST-TE STTT24,:2 3 ,‘23 

2A TEST=TEST1 
TND=L 

23 CONTINUE 

IF( IND-K )25,22»25 
25 DO 26 J=1,N 

"rESTT=ATTTND» J) ' ~ 

A I ( IND, J ) — A I ( K,J ) 

26 AI(K,J)=TEST1 
DET=-DET 
LL = IORI)( K ) 

I GRD( K )= I ORD ( IND ) 

TORT)'rTND)^U: 

22 CONTINUE 
IF(AI(KtK) )31, 32*31 

31 AI(K,K)=1./AI (K,K) 

DO 2 J=1»N 
I F ( J-K )3 , 2 , 3 

‘ 3 £JUOr=ATmJT#An“K,KT ~ 

AI(K»J)=QUOT 
IF(QU0T)42,2,42 
42 DO 4 1=1, N 

I F ( I— K )5 , 4, 5 

5 TEST=AI( I*J)-AI( I,K)*QUOT 

I F ( ASST TEST) -ABS ( ATT I , 3 IT* 1 .DE -5 ) 70 ,71,71 

70 A I ( I , J ) = 0. 

GO TO 4 

71 A I ( I , J )=TEST 

4 CONTINUE 

2 CONTINUE 

do n i=trrr - ' 

IF( I“K ) 12, 1 1, 12 
12 AI(l,K)=-AI( I,K)*AI(K,K) 

11 CONTINUE 

DET=DET/AI ( K , K ) 

1 CONTINUE 

55 '"DO 52 I =1 »”N 

I ND= I ORD ( I ) 
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I F { IND-1 ) 53» 52» 53 
53 LL = IDRO( I ) 

IQRDC I )=I0RD( IND) 
IORO< I ND ) =LL 
□0 54 K=1 ,N 
TEST=A I( K » I ND ) 
AI(K.IND)=AI<K,[ ) 
54 A I ( K , I ) =TE ST 
GO TO 55 
52 CONTINUE 
33 RETURN 
32 DET=0. 

RETURN 

END 
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or 



M.Qc) - M .(k) | 

\j A 2 |H i (k) - ML (k) | 



The program is coded such that up to 12 inventory areas can be compared 

simultaneously, i.e., X- . (i»3 = 1, 2, ... 12) 

i » 1 


DEFINITION OF INPUT PARAMETERS : 

CARD 1 

NCHAN = Total number of channels. 

NTAU = Number of channels used. 

ICHAN = Reference channel used in correlations. 

ICHEK = 1 for absolute value calculations of input. 

ITERM = 1 for each additional case or data set combination. 

CARD 2 

ICR0P1 = Identification of data set 1. 

JIDEN = Output identification of ICR0P1. 

ICR0P2 = Identification of data set 2. 

KIDEN = Output identification of ICR0P2. 

CARD 3 

CKN1(I) = Means for data set 1. 

CRN2(I) = Standard deviation for data set 1. 

WHTl(I) = Means for data set 2. 

WHT2(I) = Standard deviation for data set 2. 


The input card setup is shown in Figure 3-16. A listing of the program is 
given below. 
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C *★* + ★★*★*****★*★★★★★★★★★★*★’*★★*★★★*★ 

c calculation of si values for SELECTION or OPTIMUM 

c_ COLOR SPACE PLOTS. 

c *****+*****★****■*★********■*+***+***+ 

DIMENSION CRNl[3,123,CRN2[5,l23,WHTl[3,121 *_W HT2 13,121, 

1 DGVl2b] ,XNUM 125 J, RATI 0125), RZERO 1251 
1 CONTINUE 


READ 80, NCHAN,NTAU, ICHAN, ICHEK, ITERM 
TTaUsNTAU 


XTAUsNTAU 

READ 64, 1CROP1, JIDEN, ICR0P2.KIDEN 

C 

C 

NCHANsNR, OF CHANNEL COMBINATIONS 

C 

C 

NTAU = TOT AL NR, OF CHANNELS 
I CHAN=REF . CHANNEL USED FOR 

CORRELATIONS 

0 

READ 50, ICRN1I1, J] , J?1,NTAUJ 


READ 50, [CRN2I1, J3 , J=1,NTAUJ 
READ 50, [WHTltl, J], J=1,NTAU3 

C 

READ 50, [WHT2E1, J3 , J=1,NTAUI 

PRINT 62 
PRINT 65 


PRINT 69, ICROP1, JIDEN 
PRINT 95 


PRINT 88, CJ,J = 1,NTAU3 

PRINT 98 \ 


PRINT 96, CCRMltl, Jl, Jsi,NTAU3 
PRINT 97, [CRN2I1, J], J»1,NTAU1 



PRINT 61 

PRINT 89, I CR0P2 , K I DEN 



PRINT 95 

PRINT 68, [ J , Jsl , NT AU 3 



PRINT 98 


.. PRINT 96, IWHTICIjJ] ,J-1,NTAU1 


PRINT 97, [WHT2I1, J3, J=1,NTAUI 
IF C ICHEK] 16,19,16 


16 DO 15 J = 1 , N T A U 


CRN2I1, J3 = ABSF[CRN2E1, J3 3 
. . . HH 11 [1, J3=ABSF[HHTiri, J 3 3 



NHT2[l,J]=A8SF[WHT2Il,J33 
.15 CONTINUE 


C 

C 

CALC, bl FOR EACH CHANNFI COMB. 


C 

19. CONTINUE 



NSTOP = 'vTAU 
K-Q 



NSTART=0 

NBEGIN=0 



I COUNT = Q, 
PRINT 62 



PRINT 67 

PRINT 90, ICROP1, JIDEN, ICR0P2,KIDEN 
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60 FORMATl3X,2HCH/3 

61 FORMAT [IX// 3 

62 FORMAT [1H1//3 

64 FORMAT EA1; I2jA1>2I21 

67 FORMAT [ 53X * 8HS I TABLE/] 

69 FORMAT t 9X , 1H* , 12F8 « 4/9X, 1H» 3 

7 0 F OR MAT t IX » 11 9H ************************ ***** ************ *********** 
!**♦*★*** ************************************************ * **** * * * * * 
1*/9X,1H*,98X,12H* *) 

71 FORMAT I18X, 1H*,12F8,3 /18X»1H»3 

80 FORMAT [613] 

82 FQRMAT[1X» I4>4X#lH*»12F8 t 5» 

1 2X/lH*>f8,3i2XilH*/ 

2_9X_dl H*f9 _8. X.i aH* *J 

83 F0RMATI53X,7HCHANNEL/ 

1 9X,1H*,98X,12H» *J 

84 F0RMAT13X,2HCH,4X,1H*.»12CI5,3X3,2X,1H*,11H MEAN */ 

1 9X » 1H* » 98X< 12H* *J 

85 FORMAT [ 51X , 11H I NPUT DATA//] 

88 F0RMAT£9X,1H*. 12115,5X3/ 9X«1H»3 

89 FORMAT £7X,3HC5-, 12, 8X,7HSECTION,5X,A1,10X#15HREF CHANNEL 6//3 

90 FQRMAT[1Xj44HDSRIVATIVE CORRELATIONS -UN-NORMALIZED - C5 - ,I2, 

12X, 4HSEC A1»1X»3HC5«,I2,2X»4HSEC A1,4X,9H30 SAM/PC//] 

91 FORMAT [1X,119H *»*** ***************** ****************************** 

±****** *************************************************** ********* 
2 *1 

92 FORMAT [ 1X///10X, 13HTRACE MEAN = F8.4//J 

_ 93_F0RM.AT(1QX»13HHATRIX MEAjj_.sF8_»_4 3 

94 FORMAT E2X,12F8,3J 

9g _F. gR K AT(55X>7HCHANN^ / 9XU H»] 

96 F0RMATC1X,8HAC MEAN , 1H*# 12F8 . 4/9X , 1H* 3 

97 FORMAT ( IX , 8HSTD DEV , 1H»# 12F8 . 4/9X , 1H* 3 . 

98 FORMAT ClX*108H**************************************************** 
±******************************************1******* *_* *** * / 9X » 1H* 3 

100 CONTINUE ~ ' " ~ ” _ " 

C ********************************************** 

IFt ITERM-13200»1#1 " 

-...20 0 - CONT INUE 
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PRINT 83 ' ' " 

_PRINT 0 4, t' j/jsl f 12] 

PRINT 70 

XSU MaQ . Q 

Y SUM* 0,6 , " 

DO 40 IJsl,NTAU ; 

K = K + 1 

_SUM = 0,U 

DO 38 J = 1 # NTAU 

XN_UiU J}_ s [WHT1[1,K3-CRN1[1,KI 3 **2* t WH T1 [ 1 , J 3 - C R N 1 J 1 , J 1 ] * * 2 _ 

DGVIJ3 s4,0*f[[HHTlCl#KJ*CRNltl»Kn**2 * IWHT2C 1 , K J**2*CRN2 tl> K 3 
1**2 3 *!t [WHTltl, J3-CRNl!l,J3 3**2 * C WHT2 ( 1, J 3 *»2 + CRN'2 U , J 3 **2 3 3 3 3 3 
R AT 1 0 [ J 3 = [2,0*XNUMtJ3 3/SQRTFtDGVCJ3 3 

I F [RAT 10 [J1 >10 0. 0 3 3 0,35, 35 

35 RATIO[J3=C.O 

TTAUaTTAU-1. 0 . _ 

SUM=SUrt+RATIO[ J3 

YSUM = YSUH*RATIQ[ J3 

GO TO 38 

30 SUH=SUM+RATIQ[J3 

YSUP=YSUH+RATIO[ J] 

IC QUNTMCOUNT + 1 

38 CONTINUE 

A V5 = SIT /TTAU 

NREGIK'sMSTART + i ' ~ ' ~ 

X 3 1 G I N = N d E G I N _ 

RZERO Ia'ST ART 3 = 0*0 

„ I F_ IK-1 JJ3 9 ,36,39 

3b CONTINUE 

PRINT o2, K, [RATIO! J3,J=1, NTAU] ,AVG _ 

NSTART=NSTART+1 

X SUN=XSUM+RATIQ[13 _ __ 

GO TO 40 

39 C ONTINUE __ 

PRINT 62, K, [RZEROI J 3 , J*1 , NSTART.] , I R ATI O [ J] , JsNBEG IN , NTAU 1 , A VG 

NSTART=NSTART+1 4 

XBEGINaN8EGIN 

NST OP = *MSTOP-1 

N=N+12 

I F [RAT IQ [NBEGINl-100 .0 332,33,33 

33 RATIO! NREG I N ] =0 , 0 

XTAUsXT Au-1, Q 

XSUM=XSUN+RATTUfNBERIN3 

30 TO 

32 XSjNaXbON+RATIO (N8EGINJ 

4 2 CONTINUE 

PRINT vl " “ ' ' " " ' 

Y TQT= I COUNT 

XSUN=X3UM/XTAU ' ' “ ~ 

_X TAU = N T A U 

TTAUaNTAU ' “ 

YSUMsYSUM/ YTOT 

PRINT 92, XSUM 

PRINT 93, YSUM 

50 FORMAT E 12F6 , 0 3 ' ' 
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Section IV 

APPLICATIONS TO MULTISPECTRAL REMOTE SENSING DATA OF AGRICULTURAL FIELDS 

4.1 INTRODUCTION 

In order to test and demonstrate the capability of the algorithms and 
associated computer programs developed, which have been described in the last 
two sections, we have employed these computer programs to process and analyze 
a set of multispectral data records which were made available to us by Dr. R. R. 
Legault of the University of Michigan. In this section, we shall present the 
results of the analysis of those actual data. 

4.2 MULTISPECTRAL DATA DESCRIPTION 

The data to be presented were obtained by the University of Michigan 
Multispectral Scanner (ref. 11). This particular set of data was obtained over 
an agricultural area near Lafayette, Indiana, from a flight altitude of 2000 
feet at 2:00 p.m., June 30, 1966. In particular, we are dealing with only a 
small portion of the area designated as field C5-17, C5-22, C5-5, and C5-8, 
respectively, as shown in Figure 4-1. C5 represents the square area about 1 
mile square with 5 marked in the figure. The number following the dash 
represents a more specific location of each field. 

The electromechanical scanner onboard this aircraft used a rotating 
mirror to scan the ground normal to the flight path. Successive scan lines 
were separated along the flight path by approximately 6 feet due to the motion 
of the aircraft. All observations were stored on analog tape. Subsequent 
analog to digital conversion used a sampling time interval that corresponded 
to a ground separation of 6 feet along the scan line. The digital data points 
thus represent the radiation intensity that were received from square resolu- 
tion elements of 6 feet by 6 feet area. 

Spectroscopic observations of each resolution element were conducted by 
passing the radiation from the rotating mirror through a prism. This prism 
spectrograph was designed to measure the reflectance between 0.4 and 1.0 p. 

The radiation from the ground refers to reflected sunlight. The spectral range 
was split into 12 different channels which are listed in Table 4-1. 
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Twelve digital time series of the spectral reflectance from the agricultural 
target were obtained. The channel number, X, of these time series is directly 
related to the optical wavelength. The sample number along these time series 
denotes the position of a resolution element along a scan line. Sample and 
scan line numbers thus provide a two-dimensional, Cartesian coordinate for the 
agricultural target. 

The data received from each of the four agricultural fields cover only an 
area of 162 by 1000 feet, i.e. 27 scan lines and 167 samples per scan line. 
Since we want to test our algorithms for inventory boundary detection and other 
compairson studies, we have deliberately patched these four fields together as 
shown in Figure 4-2. 

4.3 RAW DATA PLOTS 

To have a quick and overall view of the magnitude distribution of the data 
record, a plotting program was used to plot the raw digital data record in 
analog form one scan line per agricultural field for selected channels, as 
shown in Figures 4-3 through 4-8. From these plots, one observes that sections 
of distinctly different reflectances exist for each agricultural field. These 
differences are particularly obvious in channels 1, 2, 11, and 12. Although 
only one scan line per field has been shown, it is representative for each 
field. This observation tends to imply that each field may consist of three 
different targets or inventory areas. 

It was also found that some resolution elements have negative value 
(about -1) for all 12 channels, which was obviously caused by human error in 
data handling prior to our data analysis. Since these negative reflectance 
cannot be true, we have replaced them when they occur by the mean reflectance 
of the adjacent resolution elements. 

4.4 GREY-LEVEL PLOTS 

Figures 4-9 through 4-12 show grey-level plots for the constructed 
agricultural field for channel 1, 2, 11, and 12, respectively. The grey level 
plots resemble very crude photographs of the target. In these plots 11 levels 
of grey as indicated by numerical numbers are used. The original data have 
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been averaged over 4 resolution elements (i.e., a area of 12 feet by 12 feet) 
first so that the number of scan lines is reduced to 52 and the number of 

samples per scan is reduced to 83. Human inspection of the plots revealed 

several areas with difference reflectances. Identifiable boundaries between 
these areas are indicated by solid lines. Many of these lines coincide with 
the constructed boundaries that separate the four agricultural fields . This 
agreement gives a good confidence in the qualitative human judgement for 
selection of inventory boundaries. It should be noted that the identifiable 
boundaries for different channels are not exactly the same, as x^ould be 
expected. Superposition of the inventory boundaries from all the channels 
is shown in Figure 4-13. The figure indicates that no boundary exists between 
two middle portions of C5-17 and C5-22, even though they belong to different 
fields. Some of the inventory boundaries are not so definite from the grey- 

level plots and they are so designated In the figure by hashed lines. Those 

boundaries in Figure 4-13 will be regarded as the expected inventory boundaries 
and will be used for comparison with automatic computer results to be presented 
in the next two subsections. 

4.5 BOUNDARY ENHANCEMENT 

Subsection 2.4 discussed the underlying principle for the enhancement of 
inventory boundaries. To illustrate the principle more vividly, the hypothetical 
illustrations in Figure 2—10 are now replaced by actual data. Figure 4-14 
shows the original' raw data for channel 11, of each Scan 13 from the 4 fields. 

It can be seen that fairly large and frequent fluctuations exist in the samples 
from each scan line. Figure 4-15 shows the corresponding smoothed data by apply- 
ing the 2— dimensional moving average over the given raw data. For the results 
shown, a moving average over 25 resolution elements in the form of 5 x 5 array 
was used. The curve in Figure 4-15 are not only smoother than the corresponding 
curves in Figure 4-14, but also reveal more clearly trends of slow variation 
existing in each scan line. Figure 4-16 shows the resulting enhanced data from 
the smoothed data using all 12 channels. The peaks in this Figure indicate the 
location of inventory boundaries in each scan. The larger the peak value, the 
greater difference exists between qny two adjacent inventory areas. 
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The computer program, which is developed according to the above algorithm 
has been applied to the multispectral data of the constructed agricultural field 
containing C5-17, C5-22, C5--5, and C5-8. The computer outputs are shown in 
Figures 4-17 through 4-21. Figures 4-17 through 4-19 employ the second type of 
boundary enhancement, while Figures 4-20 and 4-21 employ the first type of 
boundary enhancement. In Figure 4-17 a 4-level plot is used which prints a 
dot (.) if the enhanced data lies below the mean plus one cf-value, prints 1 if 
it lies between the mean plus one o-value and the mean plus 2 a-value, and 
similar definitions for prints 2 and 3. It can be seen that too many boundaries 
appear in this plot. Figure 4-18 shows a new plot without all those symbol l’s. 
One can now see very distinct inventory boundaries which correspond very well 
to those inventory boundaries determined by the human photo-interpretation 
(Figure 4-13). Figure 4-19 employs only 4 channels instead of all 12 available 
channels. Close resemblance of this plot to that in Figure 4-18 would imply 
that these 4 channels may be sufficient to represent all 12 channels. Figures 
4-20 and 4-21 are counterparts of Figures 4-18 and 4-19. Very close resem- 
blance between these two sets of plots further indicates that the two types of 
boundary enhancements, i.e., either taking the absolute value or taking the 
square of the difference between the 2-dimensional moving averaged data are 
about the same. For economy of computation time, the second type of boundary 
enhancement seems to be better. In order to see the effect of the 2-dimensional 
moving averaging on the enhancement of boundary, we have employed another 
moving average for 9 resolution elements in the form of a 3 x 3 array. The 
resulting inventory boundary map is shown in Figure 4-22, Comparing this with 
Figure 4-18 one can see that the major boundaries are almost the same in the 
two figures, and that Figure 4—22 shows more isolated boundary spots than 
Figure 4-18. This implies that using the 5x5 array of resolution elements 
for moving average yields better and more clear-cut boundaries than using the 
3x3 array. This clearly indicates the importance pf the 2-dimensional 
moving average on the performance of boundary enhancement. 

4.6 INVENTORY BOUNDARY DETECTION BY STATISTICAL VARIANCE TECHNIQUE 

Figures 4-23 and 4-24 show the inventory boundary map obtained by the 
statistical variance technique using channel 1 and channel 11, respectively, 
with 99 percent of confidence limit. The symbol V represents the locations of 
inventory boundaries as detected in the scanwise direction. The symbol H 
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represents the location of inventory boundaries as detected normal to the scan 
line. Finally, the symbol X represents the location of inventory boundaries 
as detected in both directions. Comparing Figure 4-23 with Figure 4-9 (grey- 
level plot) , and Figure 4—24 with Figure 4—11, one can see that the inventory 
boundaries determined automatically by the statistical variance method compare 
reasonably well to those selected by the human photo-interpretation. 

Figure 4-25 shows the inventory boundary map using channels 1 and 11 
together. Comparing this figure with Figures 4-23 and 4-24, one can see that 
the new inventory boundaries fall roughly in between those indicated by 
using the individual channels. Two more maps of inventory boundaries using 4 
channels (1, 2, 11, and 12) and 12 channels, respectively, are shown in Figures 
4-26 and 4-27. Figures 4-25 and 4-26 are quite similar to each other. The 
only major differences between Figures 4-26 and 4-27 appear in field C5-8. 
Finally, comparing Figure 4-27 with Figure 4-13, it is seen that good agreement 
in inventory boundaries are observed for C5-17, C5-22, and C5-5; and that the 
vertical boundaries in C5-8 are completely different in these two figures. This 
poor performance detection in C5-8 might be due to the existance of a vertical 
boundary in C5-8 locating very close to the left end of the field, where the 
scanwise inventory boundary detection is initiated. It is recalled that at 
least 4 samples need to be used for the initial estimate of the statistical 
variations (see subsection 2.5) prior to testing for boundaries. This require- 
ment makes it very difficult to detect the left vertical boundary in C5-8 and 
further cause misdetection later. In this respect, the boundary enhancement 
method (Figure 4-22) is more powerful then*fehe statistical variation method 
for detecting inventory boundaries. 

4.7 STATISTICAL FEATURES OF INVENTORY AREAS < 

Once the inventory boundaries within a given target area (or image) are 
determined as presented in subsections 4.5 and 4.6, the statistical features of 
each individual inventory area, which are the areas bounded by the inventory 
boundaries, may be calculated. For convenience of discussions, we shall 
designate various inventory areas as indicated in Figure 4-13 . 

I 

Figure 4-28 is a computer printout of the mean vector, correlation 
matrix, covariance matrix, and normalized covariance matrix of the inventory 
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area C5-17 Section A, The magnitude of the covariance matrix is about two 
orders of magnitude less than that of the correlation matrix. This implies 
that the standard deviation for each channel is about one-tenth of the corres- 
ponding mean value. The fact that roughly the average value of the element 
of the normalized covariance matrix is about 0.50 implies that the 12 channels 
of the data record are statistically correlated within each inventory area. 

Figure 4-29 shows the mean vectors of reflectance for 6 inventory areas 
in fields C5-17 and C5-22, while Figure 4—30 shows the mean vectors of reflect- 
ance for 6 other inventory areas in fields C5-5 and C5-8. It is noted that the 
mean vectors of C5-17, Section B, and C5-22, Section B, are much closer to each 
other as compared with other mean vectors. This agrees with the computer 
results that no inventory boundary is found between these two inventory areas. 

We shall next examine the univariate probability density functions of the 
multispectral agricultural data. It is recalled that among the 12 available 
channels of data for the four agricultural fields under consideration, we have 
found that only channels 1, 2, 11, and 12 are almost sufficient for determing 
inventory boundaries. Thus, the discussion shall be restricted to the univariate 
probability densities of only these four channels. Figures 4-31 through 4-34 
show the univariate probability histograms for field C5-22 for channels 1, 2, 

11, and 12, respectively. The curves marked by symbol A are obtained from 
actual data, while the curves marked by symbol B are gaussianly distributed 
curves with the same standard deviation and peak value as the actual data. 

The figures clearly indicate that the reflectance over field C5-22 is definitely 
not gaussianly distributed, and that the actual distribution are multiple-mode. 
For channel 1 we can see two distinctly large peaks. For channel 2, there are 
three main peaks. Figures 4-33 and 4-34 show that the distribution is skew to 
the low reflectance value. In ojrder to provide one more check to see whether 
the inventory areas we have established are really composed of statistically 
homogeneous samples, we have also plotted the probability histograms for C5-22, 
Section A, B, and C with channel 1 and 11 (Figures 4-35 through 4-40). These 
figures clearly indicate that the univariate probability histograms are 
practically Gaussian for all the three inventory areas under consideration. 

In other words, field C5-22 is composed of three different homogeneous areas, 
each of which may be approximated by a different gaussian distribution. 
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Therefore, the division of the field C5-22 into three inventory areas as shown 
in Figure 4-13 is reasonable. 

4.8 DIVERGENCE MATRIX AMONG INVENTORY AREAS 

To calculate the divergence matrix among various inventory areas (or 
populations) and thus to compare the statistical dissimilarity (or separability) 
among them, we need to input the mean vectors and covariance matrices of all 
inventory areas into the divergence matrix program. Tables 4-2 and 4-3 show 
two such inputs for C5-17, Sections A and B. Table 4-4 shows the divergence 
matrix for the 12 inventory areas as indicated in field designations. (Note 
that the larger the value of divergence between two inventory areas, the less 
similar they are.) We have observed in an earlier part of the Section that 
C5-17, Section A is more similar to Q5-22 Section B than any others. The table 
gives a divergence of 3.303 between these two inventory areas, which happens 
to be the smallest element of the divergence matrix. Tables 4-5 and 4-6 show 
two contributing parts, called trace 1 and trace 2 matrices, to the divergence 
matrix; they represent the first and second term (without the factor 1/2) of 
the sight-hand side of equation (2-11). If covariance matrices of all inventory 
areas are identical, then trace 1 matrices will vanish. On the other hand, if 
the mean vectors of all inventory areas are the same, then trace 2 matrix will 
vanish. Comparing Tables 4-5 and 4-6 indicate that the elements of trace 1 
matrix are much larger than those of trace 2 matrix. This would imply that 
the dissimilarity among those 12 inventory areas is mainly due to their 
difference in mean vectors rather than their covariance matrices. 
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Figure 4-1. AERIAL PHOTOGRAPH OF THE AGRICULTURE FIELD 
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Figure 4-2. A CONSTRUCTED MAP OF FOUR AGRICULTURAL FIELDS 
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Figure 4-3. RAW DATA PLOT OF CHANNEL 1 FOR ONE SCAN LINE PER FIELD 
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Figure 4-4. RAW DATA PLOT OF CHANNEL 2 FOR ONE SCAN LINE PER FIELD 
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1000 FT 


Figure 4-5. RAW DATA PLOT OF CHANNEL 6 FOR ONE SCAN LINE PER FIELD 
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Figure 4-6. RAW DATA PLOT OF CHANNEL 8 FOR ONE SCAN LINE PER FIELD 


4-14 






TR-860 



Figure 4-7. RAW DATA PLOT OF CHANNEL iT FOR ONE SCAN LINE PER FIELD 
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Figure 4-8. RAW DATA PLOT OF CHANNEL 12 FOR ONE SCAN LINE PER FIELD 
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NORTHROP 


TR-860 


HUNTSVILLE 


CHANNEL 1 

HAVE LENGTH (MICR0NS) 0. 40 - 0.4* 


' SYHB0L 


1 

2 

3 

V 

5 

'6 

7 

3 

9 

0 


interval 
0.50 - 0.65 

0.65 - 0.80 

0.80 - 0.95 

0. 95 - 1.10 

1. io"-' 1.25 

1.25 - 1.60 

1.40 - 1.55 

1.5_5_- 1.70 

' 1.'70 - 1.85' 

1.85 - 2.00 

2.15 - 


"CESPT" C5 :r !’7,5-22*5-5,5-8i SlTTTODE 2000 FT* OATE REC03OEO 6/30/66 


11 1111 11 112 2 2222 222 2 3333 3 3 33 33 4 ******4**555 5 55 55 556 666666666777777 77 T78 888 
12 345678901 2365*67 8901236567890123456789012365678901 234567 890 1234 567890123 556 789 012 3* 

***♦****»***********************$********* ****** 6 * 64*4 ********♦*<:******* *6 ♦♦****■***♦ 


1 4555555556665555556555545554, 

2 *555455555655556555555445555 

3 *555445555665554555555545555 

4 *555455556665555554555444555 

5 *555553555555555555545544555 

Qg iy 6 *555455655555555555555545555 

8 *555555556555455555555444455 

9 *555555556555565555554444555 

10 *555555555554454555545544555 

1 1 *5555 55 5 555 $545455553 5454455 

12 *555455555554544454454444455 

4444444444333443333333343334433444 

4444333343343433333333333333333434 

4343444343335333333333334333333334 

4434344333333333333323434333333334 

3444333443433333434333333333333334 

4444444433334343434343433433343334 

5565665665666766677767 
5655556556666666666666 
5555555556666666666666 
5556555565566666666666 
555 5556655665 666666766" 
5556556566666666567666 

4334344433433334333333334433333334 

3444443333333343434334343433333334 

4*3***3*3*3*33*33*33333****333**35 

4444444444333334443333344434334334 

3*33*433***33333**333333333333333* 

3434444444333333343333364333333334 

5555565655565666666666 

5556555666666666666667 

S5655~$55666666T6666777 

5566555555566666666666 

5556566566666*666666666” 

5555556556565666666766 


ttmmxrawrararai********* 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ■! 

14 *7788888888888887777776 

15 *778877778 aTt 7 887' 7 776 76 

16 *7778778787888887776776 

17 *77787888'ff7777778S7m? 
IB *7777887877877787777676 

44554454444444444434433334343444445 
44444444344433333443333344434333345 
4454454444444444434433 3343344 344445 
44455444444443344434333434443443434 
44444344444344443443333334333343334 

[77677777777777777766656555 
&6666767676776777 76665 5555 
46776666676676667666655555 
i7T66667666777676766555555 " 
667776667666667676665555*5 

19 *677 7/8 880UUH /BB /77777S 
no 20 *7787787777777777777766 

C5 -id. 21 *777777 T777TST7T875775 7 

22 *7787877887778788777776 
"23 *788888T8T7777T877TTS77 
24 *788887-6767777878767667 


b667666T67#76777766655?555 
56666666666776766766555545 
56? 66667776676666 66655 55*4 
566667666666666666665545*5 
566666676777776666666555*5 
56666666666667666666555555 

?5“* ty fi 8 88 3 B 7 ?6 8 7 &7BY f i 77 fi 
26 *7777887077877877777677 

444544444434*4* 

444444434433354 

434434333333*6344464 

443334433333*3343344 

343443434T4445T55445 

4444433434*343434445 

1566676677676776666665*55*5 
556666666676666666565*554* 
555666566666565s 6676666555 
56665665556556566666666566 


55444 

44443 
«54* 

44444 

4443333343 

4*34334***44444*3 3333334344*34 
4444444443434 33434434444544345 

55665656655656655677665665 
5586S566566666 666777666665 


44444544444444334443433344443444335 
4*(S33***4*X3***A4***3 3*4*34* 3*3*35 

56666665656696585666655566 

66656565666656656566665655 


4444444444343"?443 344333333343 343444 
44444544454334444443333444434444445 

S75 665 6666656655 7666665666 
56666665666665656676667566 


444544444444434334443343434^4434335 
44454444433 34443444343 3334433 333444 

5 6 5655 6*6455566066 76666656 
^56655566656666666 77677655 

3/ *6666676666666656536555 
38 *6765677666666666665566 

44444344443 344344 3 5343 44344444444"5? 
54454444443443434444433444434443345 

66674355656566566666766669 
56665556666665 66 6676665655 

1 ■ 1 1 1 1 1 > 1 1 1 1 1 1 1 1 1 1 1 1 II 


56666665665566666776676656 

40 *4444544 


554555*56666666665666676566*365567 

4**3*****34* 

41 *4445553 

42 *4444554 

43 *4444453 
64 *4444543 

333333333333523222332322332333 

334333333332323i23333S33323334i 

333332322223232222323223323333 

433333332323332322233232333333 


44 ** ******53 
45*3*4****** 
4*34***33333 
434**44**43* 

45 *4455554 
C5_8 46 *4545554 

*7 “*445455 4 
48 *4545554 

343333333323333323233333333333 

333232333333333323233323223333 

55**555666656666556666666555566668 
454**55 56565665566666776666556655a 

4444**4**44* 

4»4********3 

334333333332332222233323333333 

433333333333333333333333233333 

554555*556655666655666666655565657 

4*455*5555566666665666665655565668 

***** *4*4*33 
4344**4***43 

49 *4455554 

50 *4554554 

433333343332333422323333333333 

433333332233232332333333333333 

45*55*5565565666656666766655556667 

*5*55*5566555556665566676555565667 

4444444444 

444444333444 

51 *4555564 

52 *5555564 

mmimmwm 

*5*5555566655666666666666656566567 

^5455556655^5866655666665555555667 

>5*44444**4* 


Figure 4-9. GREY-LEVEL PLOT {Channel 1) 
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HUNTSVILLE 


_CHANNEL 2 
HAVE LENGTH (MICRONS) - 


0.44 - 0.46 


CROPS C5-l7»5-22,"$-5,5-8V' 


SYPBBt. 

t 

2 

3 

4 

5 

6 ~ 

7 

8 
9 
0 

ALTITUDE 2000 fT. 


INTERVAL 
1.00 - 1.20 
1.20 - 

1.40 - 
1.60 - 
1.80 - 
2.0 0 - 

“ 2.20 - 

2.40 - 
2.60 - 
2.80 - 
3.20 - 


1.40 
1.60 
1.80 
2.00 
2.20 

2.40 
2.60 
2.80 
3.00 


OATE RECORDED 6/30/66 


C5-17 


C5-22 


C5-5 


C5-8 


111111111 12222222222333333333344444444445555555555666666666677777777778888 
12345678901234567890123456789012345678901234567890123456789012345678901234567800123 

********** ******************************************** ************i««********<.*****40 

— 1 *44^3444 5 5^ 4445 44 444 433 3 443J3 333 332 322 32 3^ 222 222 22 22 222222 2 3 31444444444 4444 4'444445 545 

2 *44444434554444444444444444433333332222222222223222222222223344344444444444444444554 

3 *44444445554344444454434444433333222222222322222222222222222344444344444444444445544 

4 *44444445544344444444333344433233332222222222222222232222222344444434444444544444444 

5 *44444444444344454344433344433332333322222232222222222222222344444444434444444454545 

6 *44444444444343445344433344433333323222222232232222232222222344444444444444444444445 

7 ****4*4445444455441544 333444 3333333233222 3'352222222232222'2223*4444344344544~444445445 

8 *44444545444444444444433334433333333233323323322232232222222344444344444544544445455 

9 *44444555544345344343433334433333333333223322222222332222222244434444444544444445544 

10 *44444444444345444344333444433333333332223323222222322222222244444344434444444445444 

11 *44444545443334444434333333433333333322322223222322322222222243444434434444444445544 

12 *44444*44 3433 344 4443 43333334 3333333233222222 32222 223332222222443443444444444444 4545 4 
_ H ****444 4453 5131454413T317344 33333 3334232 322T3TZ2222 31222 222234 33444444 344*4544445544 

14 *677777776766776656567^344*33333333333342323222333233333^5655555555555555554*4343345 

15 *6777666766667776566677334343433323332332322222222232233: 55555555555555555544*333335 

16 *7767676676766676665587 334343333333333332222222232332233* 55555555555545455554*433335 

17 *6676 777667767777666667 33343343333233323332233223233 3232 3 45655555555S55545544*433334 

18 *677767 666666667656667* 334433333332333233333222222332233; 45555555555555555544* 433335 

' 19 46V6677 777 nT77 766 65667 34443*433 3lY33ll3l3T?223?3333 32 32": 4555555 5555555555544* 333335 " 

20 *6667676666666666656667334334333323333333222232323333233:55555555555555555544*333334 

21 *6666767766766767656667 3333333333323333332222 32323333233 -55555555555555S45544* 333334 

22 *6776667666767776665676 334443433323333322222323322322233* 4555545555545555*544* 333334 

23 *7777767666777666655675344443333332333332322222333223333*55555555555555655544*333335 

24 *6777666766666777666676 3344333323333333233222222223232233 55555555555555545444* 333335 

25 "*7777"76TT66TSi'67 6666561 5T5334323333T333Z5%2322223233'2233 3 46555555555555554443* 333334 
26 *6677777777776776656561 3344433233333333233333222323222333 35555555555555444444*433334 

27 *555665555555545555545* 334333333333332333333322223332333: 54554444444444444565444444* 

28 *556565555555555554*457 3433333333333333333333333333233333 555544444444444445544444444 

29 *5 555 55 55 5 55 5445 55 5445 : 333333332332323323222322222233322: ***544444444444444554443444 

30 *555565655545545554545* 343333333332233333232332323333333: 545544444444444445555444444 
“31 *5 5565155 553 54 5 4 554 5 45 2 54T53l5Tl33332333T11JyS222333Z3233 54*544445444444445544444444 
''32 *5655665555445445455453 3433333333322333333323322233323333 5445444445444444*5544444444 

33 *5555655555545555555453 3543333333323333332233222232323332 545544444544444445544444444 

34 *555566555555545555545* 3 3433 3333 3333 3 3 33 3 3233 333 3 3J22 33 333 545644444444444445554444444 

35 *555 5 55 5 5 55445 5 44454557 33433333333223332 33 333 32 33~333 3333 3 54§544444544544445544444444 

36 *555555555554444455545* 343333333333233333333332333333333: 54*454444544444445544444444 

37 *55555663555 55544445453 5531313333333332251333 33233323333: 544554445444444445554544444 

38 *S55 555 55545554455 5446* 344313333333233333333323333333333 ; 5454564444.444444445^5^4646 
_ 39 *665 55355444454 54445453 3433333333323333313333332333333333 5«46644465 : 4444464555f544444 

*J.» J343^2222222Z212222r22l2222222 223^334334 5544«§545445655"654434454342{232Z2222Z22 , 2 

41 *3 33 3432 32222727772 1272?r2I?Z2IZ22?224f5$33444 44$5555314555S? ‘>54446644372 232'32'22222'21 

42 *3334433 222222222222222212222222222233 5333344445443554445655554*33344363 23 3222222222 

43 *3334442222222222222222222222222222233533434445544555544465655454345445*223222222222 

44 *33344332222212212221222122222222222235334334454445545455555554434444362232222222222 

45 *3334442 32222222 1222222222222222222223 53344'34555S4S5S4445656'5544*44Y4455 232232322222 
j 3222222222222222^2222222222222 543333444545554444 5566 S5544444444S 123232222222 
i 322222222222222222222222222222 5433344445445444445665545433454446 223232322221 

48 *333 4443 222222222222222222222222222223 5434434444455554445566544434444446 232332232222 

49 *33*444: 222222222221222222222222222223 S43333444444455444S566S5S434444446 222233322322 

50 *3344443 322222222222222222222222222222 53344444444445_55544 5665 5 4* 33454446 232 22323 3222 

51 *3334442 2222222222223222122222222222225334444454445454445555554443444446 222222222222 

52 *4 4 4 444 322 222222,2222333222122222222222153 3 4444 555 444555 445565544 44344444612332 23122232 


46 *3334343 

47 *3334433 


Figure 4-10. GREY-LEVEL PLOT (Channel 2) 
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WAVE LENGTH (MICR0NS") 
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- 
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4 

0.80 
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- 
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:‘R0P S C5-17VS -‘22 , 5-TT5=8 , 


ALTITUDE 2 00C FT. 


OATE RECORDED 6/30/66 


C5-17 


C5-5 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

15 

14 

15 

16 
17 
1 


11111111 112222222222333333333344444444445555555555666666666677777777778888 
12345678901 23456TB901234567890123456789C1234 5678901 23456789012345678901234567890123 
************************************************************************************* 

1 *555 566 65556665465555 5566544(33 3332 333323 32233 22 32222 3222 33333 3^4334333 33434344444444' 
*555556 555 565564 55 5555565545 33332332233223232232222232222233233 333333343434434444444 
*55556655556665456555556553523332222232232222222222222233332334433333434334344444443 
*555566555666655555555666545 33322223222222323323223222322323333 433343443333344444444 
*555556 555666555 55 555 5666 54E 333223 32222322 322 3223922 32232233333 333343333434444444444 
*5555565 55665555 555555666545 3332 3222222222222 2 23 3 2222 3222 322233 44 33333 3334 3344 44 3344 

*55555555566655556665556665533333333323322323222222222333233333343434343334434444444 
*555555555665555555555666645 33333323022222222222222323333233233 443333333334444444444 
*55555555566655555565566564533333233332223223233233222222232233332433333344344444443 
*55555555556655665565667665523333332322322323233222332333332333343443343334344443343 
*555565555665555555656677645)33333333223222223232 333333333333 233 33333334333344 444444 3 

?2 5¥33Tf3TJ 3 3232? 5TSy523Z52335532 33 3 3 343 33 3343434444434 3 


C5-22 


27 *444444444444444444454: 

28 *444444444444444444454' 
*444444444444444444454: 
*444444444444444444454: 

'in * 4«47 

32 *444544444444445454454i 

33 *44444444444<444”444454£f 
*444444444544444444454! 
*4444S4444S4444"4 4444546(7 
*444444444444445444454! 

'“ 37 ~* 34 ~ 445444'45 44444 444454 . 

38 *4444444444444454444545) 

39 *444444444544444444454! 


29 

30 


34 

35 

36 


*55455555444444445444433333333332322332333222333333333333433333334333333333333444443 
*45544455444434444445433333323332223222222222323333332333433333333333333332332444433 
*45555455554444444454543333333333332322333223333333333333433333333333333333332344443 
*4555545444444444444544(3333323323223222222222233333333333 433333333333323333323344443 
3333332233222332222222333223333333433333333333333333333334433 
f 9 * 555^555544444444444544 S33223'5333"'32233322?232 323332333332 453433333333333333‘333334443 

20 *4545545544444444444544333233323323233222222333323333333:433333333333333323333333433 

21 *5555545544444444444544333333333323223232233322233333333:433334333333333333333334443 

22 *55555554444444444444443333333232333333232232333333333333433334333333333333223344443 

1333332233222333233232333333333233:433343333333333333223344433 

24 *5444444455444444444544333332332332323322223332333323233:433333333333333233222333343 
"25 *44!T537l445444447t 4445544 33"3332*23333 33235323333 323323322333 4333343*33333333322333344433 
26 *55554444444444444455543333333322223222323333333333332333333333333333333333223343443 


5(7666776665666667666655 56777777666 
66677666567766776666566677777665 
5(776777667567766766556666777767666 
6P6677 7676 6&777677 65S5556777777766 

6(6667777765677777765666 6677 77777676)! 
(7767777765677777766655 6677T767776!ej 
76677777678777776666666777777766^ 
7666677667777777666566777 Tf777676. 
5f7766767666787777766766<r7'7777777766 
7 T7MSSeSh~£TrmTni'(, 6 6 A 7 7 7 77777766) 
776666766677776777666667777767766 
5j77766666S67 87 767766 66 6 677 767? 66766) 


5*7 


C5-8 


40 *322233333333332332323232323322222223333 

41 *3222333333373323232322322232222222 J2Z337 

42 *323233333333323332333232222333323333233 

43 *333323333333332233333332223333222223333! 

44 *332233333 3333 22323233333322323333233333: 

45 *33223333333322232*3322332323232222233333: 

46 *3 222 33 433_3 33 222 32 233333 3333 223 22 2_32 22 33 

47 *332233333233323332323333323223323323233 

48 *323233333322222233322323232322223222233 


9444443444443343344444433323 
9444444444444343443344433333 
3444 44443 4444444444444443 3 33 
6*444444343444343444344333323 
(4544443444 43444*44 3444333232 
(544444444444344344444343333 
444444344444344344444343333 
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Figure 4-11. GREY-LEVEL PLOT (Channel 11) 
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50 *3 323334-33333322223322 3333Z2222322222333 3(4443334-34444333445454333 33 3 343223 3 3 33333 323 
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Figure 4-12. GREY-LEVEL PLOT (Channel 12) 
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Figure 4-13. 


INVENTORY BOUNDARIES OBSERVED FROM GREY-LEVEL PLOTS 
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Figure 4-14. ORIGINAL RAW DATA FOR CHANNEL 11 ( X = 0.72 - 0.80 y) 
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Figure 4-1 7. INVENTORY BOUNDARY MAP BY BOUNDARY ENHANCEMENT, 
CHANNELS 1-12, USING ABSOLUTE VALUE 
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Figure 4-18. 


INVENTORY BOUNDARY MAP BY BOUNDARY ENHANCEMENT, 
CHANNELS 1-12, USING ABSOLUTE VALUE 
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Figure 4-19. INVENTORY BOUNDARY MAP BY BOUNDARY ENHANCEMENT, 
CHANNELS 1, 2, 11, AND 12, USING ABSOLUTE VALUE 
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Figure 4-20. INVENTORY BOUNDARY MAP BY BOUNDARY ENHANCEMENT, 
CHANNELS 1-12, USING SQUARE VALUE 
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Figure 4-21. INVENTORY BOUNDARY MAP BY BOUNDARY ENHANCEMENT, - 
CHANNELS 1, 2, 11, AND 12, USING SQUARE VALUE 
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Figure 4-22. 


INVENTORY BOUNDARY MAP BY BOUNDARY ENHANCEMENT, CHANNELS 1-12, 
USING ABSOLUTE VALUE AND 3 X 3 ELEMENTS FOR MOVING AVERAGE 
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Figure 4-23.' MAP OF INVENTORY BOUNDARIES (Channel 1) 
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Figure 4-24. MAP OF INVENTORY BOUNDARIES (Channel 11) 
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Figure 4-25. MAP OF INVENTORY BOUNDARIES (Channels 1, 11) 
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Figure 4-26“ MAP' OF INVENTORY BOUNDARIES (Channels 1, 2, 11, 12) 
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Figure 4-27. MAP OF INVENTORY BOUNDARIES (Channels 1 through 12) 
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Figure 4-28. MEAN VECTOR CORRELATION MATRIX, COVARIANCE MATRIX, AND NORMALIZED 
COVARIANCE MATRIX OF C5-17, SECTION A 
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Table 4-3. MEAN VECTOR AND COVARIANCE MATRIX OF C5-17 SECTION B 
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Table 4-4. DIVERGENCE MATRIX USING ALL 12 CHANNELS 
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Table 4-5. CONTRIBUTION TO THE DIVERGENCE MATRIX FROM COVARIANCE MATRICES 

TRACE 1 MATRIX USING ALL 12 CHANNELS 
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Table 4-6. CONTRIBUTION TO THE DIVERGENCE MATRIX FROM MEAN VECTORS 


TRACE2 MATRIX USING ALL 12 CHANNELS 
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Figure 4-29. MEAN SPECTRAL VECTORS OF THREE HOMOGENEOUS SECTIONS IN FIELDS 
C5-I7 AND C5-22 
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Figure 4-30. MEAN SPECTRAL VECTORS OF THREE HOMOGENEOUS SECTIONS IN FIELDS 
C5-5 AND C5-8 
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Figure 4-32. UNIVARIATE HISTOGRAM OF CHANNEL 2 OF ENTIRE FIELD C5-22 
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Figure 4-33. UNIVARIATE HISTOGRAM OF CHANNEL 11 OF ENTIRE FIELD C5-22 
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Figure 4-34. UNIVARIATE HISTOGRAM OF CHANNEL 12 OF ENTIRE FIELD C5-22 
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’ Figure 4-35. UNIVARIATE HISTOGRAM OF CHANNEL 1 OF C5-22, SECTION A 
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Figure 4-36. UNIVARIATE HISTOGRAM OF CHANNEL 1 OF C5-22, SECTION B 
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Figure 4-37. UNIVARIATE HISTOGRAM OF CHANNEL 1 OF C5-22, SECTION C 
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Figure 4-38. UNIVARIATE HISTOGRAM OF CHANNEL 11 OF C5-22, SECTION A 


4-51 


NORTHROP 


TR-860 


HUNTSVILLE 



Figure 4-39. UNIVARIATE HISTOGRAM OF CHANNEL 11 OF C5-22, SECTION B 
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Figure 4-40. UNIVARIATE HISTOGRAM OF CHANNEL 11 OF C5-22, SECTION C 
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Section V 

SUMMARY AND CONCLUSIONS 


The most significant results of this study are the development of two 
statistical algorithms for determination of inventory boundaries within a tar- 
get area from multispectral remotely sensing data without reference to any 
ground truth or training data. The samples within each individual area bounded 
by the inventory boundaries, called an inventory area, constitute a homogeneous 
population. The first algorithm is based on a principle of boundary enhance- 
ment which involves area smoothing and taking absolute value of the 2-direc— 
tional difference of the smoothed data. The second algorithm is based on the 
statistical variation of the multispectral data and the confidence limit of the 
statistical estimation. The two algorithms are thus basically opposite in 
their underlying principles, since the first relies more heavily on the mean 
of the data; while the second relies more heavily on the variation of the data 
about the mean. As such, they really compensate each other. These two algorithms 
have been implemented into digital computer programs. 

In addition, several other computer programs have also been developed for 
general statistical analysis of multispectral data. Specifically, these programs 
are for calculating the mean vector, correlation matrix, covariance matrix, 
and univariate probability density distribution from the inventory area, and 
for calculating the divergence matrix among various inventory areas. The 
divergence is a statistical measure for separability or dissimilarity between 
inventory areas. Also developed is a grey level plotting program for 2- 
dimensional digital display of a single-channel record which may be either raw 
data or final classified output data. 

The capability of the above computer programs for automatic determination 
of inventory boundaries has been successfully demonstrated on a set of remotely 
sensed data obtained by the University of Michigan Multispectral Scanner over 
some agricultural fields at a flight altitude of 2000 feet. 

The most important advantage of the data analysis methods is data com- 
pression which is accomplished by replacing voluminous raw data samples by the 
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location of inventory boundaries and associated aerial average over the inventory 
area. Two to three orders of magnitude of reduction in data flow rate and 
total computation time may result. 

One drawback of the two algorithms for inventory boundary detection is that 
the boundaries sometimes do not close. To remedy this drawback, a new algorithm 
for automatic sequential classification of multispectral data into homogeneous 
populations (without reference to any training set) has been proposed (ref. 13). 
The performance of the new algorithm will be reported in the future, for it 
has not yet been implemented. 
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